Method DiffJacobi

fem/coefficient.cpp:1939–1954

  }
  
  shared_ptr<CoefficientFunction> DiffJacobi (const CoefficientFunction * var, T_DJC & cache) const override
  {
    auto thisptr = const_pointer_cast<CoefficientFunction>(this->shared_from_this());
    if (cache.find(thisptr) != cache.end())
      return cache[thisptr];
             
    if (this == var)
      {
        if (this -> Dimensions().Size() == 0)
            return make_shared<ConstantCoefficientFunction>(1);
        return IdentityCF(this->Dimensions());
      }
    auto res = scal * c1->DiffJacobi(var, cache);
    cache[thisptr] = res;
    return res;
  }
};