Method DiffJacobi

fem/coefficient.cpp:6262–6274

    }

    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 IdentityCF(cf_then->Dimensions());
      
      auto res = IfPos (cf_if, cf_then->DiffJacobi (var, cache), cf_else->DiffJacobi (var, cache));
      cache[thisptr] = res;
      return res;
    }
  };
  
  extern