Method DiffJacobi

fem/coefficient_matrix.cpp:860–883

    }

    shared_ptr<CoefficientFunction> DiffJacobi (const CoefficientFunction * var, typename BASE::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());

      shared_ptr<CoefficientFunction> res;

      if (this->Dimensions()[0] == 2)
        res = (TraceCF(c1)*IdentityCF(2)-TransposeCF(c1)) -> DiffJacobi (var, cache);
      else if (this->Dimensions()[0] == 3)
        {
          auto trcf = TraceCF(c1);
          auto prodcf = c1*c1;
          res = (0.5*(trcf*trcf - TraceCF(prodcf))*IdentityCF(3) - trcf*TransposeCF(c1) + TransposeCF(prodcf)) -> DiffJacobi (var, cache);
        }
      else
        res = (DeterminantCF(c1) * InverseCF(c1)->Transpose() ) -> DiffJacobi (var, cache);
      cache[thisptr] = res;
      return res;
    }

  
  };