Question regarding the implementation of GridFunction.Set

sblauth · September 9, 2022, 12:15pm

Hello,

I’ve got a question regarding the GridFunction.Set method. Here (1.2 CoefficientFunctions — NGS-Py 6.2.1907 documentation) it is stated how the method works. As I am interested in the exact implementation of this method, I would like to know, whether someone could point me to the corresponding piece of the source code.

Thanks a lot in advance,
Sebastian

joachim · September 9, 2022, 1:40pm

HI Sebastian,

it is done in SetValues:

github.com

NGSolve/ngsolve/blob/8fe25e5e63b0d30628f5b84144e75a0c0fbb5261/comp/postproc.cpp#L384


      
                                        bool applyd,
                                        LocalHeap & lh,
                                        int component);
              
          
          

          

          

          

          
  template <class SCAL>
            void SetValues (shared_ptr<CoefficientFunction> coef,
          		  GridFunction & u,
          		  VorB vb,
                            const Region * reg, 
          		  DifferentialOperator * diffop,
          		  LocalHeap & clh,
                            bool dualdiffop = false, bool use_simd = true, int mdcomp = 0,
                            optional<shared_ptr<BitArray>> definedonelements = nullopt,
                            int bonus_intorder = 0)
            {
              static Timer sv("timer setvalues"); RegionTimer r(sv);

follow the path dualdiffop==false

best, Joachim