Method Diff

fem/coefficient_matrix.cpp:171–178

    }

    shared_ptr<CoefficientFunction> Diff (const CoefficientFunction * var,
                                          shared_ptr<CoefficientFunction> dir) const override
    {
      if (this == var) return dir;

      auto thisptr = const_pointer_cast<CoefficientFunction>(this->shared_from_this());
      return (-1) * thisptr * c1->Diff(var,dir) * thisptr;
    }

    shared_ptr<CoefficientFunction> DiffJacobi (const CoefficientFunction * var, T_DJC & cache) const override
    {