Method DiffJacobi

fem/coefficient.cpp:5105–5129

  }

  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());

    Array<int> resdims;
    resdims += Dimensions();
    resdims += var->Dimensions();

    Array<int> resstride(resdims.Size());
    resstride.Range(0, dims.Size()) = stride;
    resstride.Range(dims.Size(), END) = 1;
    for (int i = resdims.Size() - 1; i >= dims.Size(); i--)
      for (int j = 0; j < i; ++j)
        resstride[j] *= resdims[i];

//    cout << "new stride: " << resstride << endl;
    auto diffc1 = c1->DiffJacobi (var, cache);
    auto res = MakeExtendDimensionCoefficientFunction (diffc1, std::move(resdims), Array<int>(pos), std::move(resstride));
    cache[thisptr] = res;
    return res;
  }

  virtual void NonZeroPattern (const class ProxyUserData & ud,
                               FlatVector<AutoDiffDiff<1,NonZero>> values) const override