Method DiffJacobi

fem/coefficient_matrix.cpp:1024–1037

    }

    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(this->Dimensions());
  
      auto diffc1 = c1->DiffJacobi (var, cache);
      auto res = 0.5*(diffc1 + diffc1 -> TensorTranspose( 0, 1 ));
      cache[thisptr] = res;
      return res;
    }
  };