Method DiffJacobi

fem/coefficient.cpp:2221–2244

  }

  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];
    int dimvar = var->Dimension();
    int mydim = this->Dimension();    
      
    Array<int> dimres { this -> Dimensions() + var->Dimensions() };

    if (this == var)
      return IdentityCF(this->Dimensions());

    auto diffc1 = c1->DiffJacobi(var, cache);
    auto diffc2 = c2->DiffJacobi(var, cache);
    
    auto prod1_ = c2 -> Reshape(mydim, 1) * diffc1 -> Reshape(1, dimvar);
    auto prod1 = prod1_ -> Reshape(dimres);

    auto prod2 = c1 * diffc2;
    auto res = prod1 + prod2;
    cache[thisptr] = res;
    return res;
  }