Method DiffJacobi

fem/coefficient.cpp:3167–3193

  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)
      return make_shared<ConstantCoefficientFunction> (1);

    shared_ptr<CoefficientFunction> dv1v1;
    int dimip = c1->Dimension();
    int dimvar = var->Dimension();
    
    auto vc1 = c1->Reshape( dimip );
       
    if (c1.get() == var)
      dv1v1 = c1;
    else
      {
        auto dvc1 = c1->DiffJacobi (var, cache);
        dv1v1 = dvc1 -> Reshape(dimip, dimvar) -> Transpose() * vc1;
        dv1v1 = dv1v1 -> Reshape (var->Dimensions());
      }
    auto res = (1.0 / const_cast<CoefficientFunction*>((CoefficientFunction*)this)->shared_from_this()) * dv1v1;
    cache[thisptr] = res;
    return res;
  }