Method NonZeroPattern

fem/coefficient.cpp:2265–2274

virtual void NonZeroPattern (const class ProxyUserData & ud, FlatVector nonzero, FlatVector nonzero_deriv, FlatVector nonzero_dderiv) const override { int dim = Dimension(); Vector v1(1), d1(1), dd1(1); Vector v2(dim), d2(dim), dd2(dim); c1->NonZeroPattern (ud, v1, d1, dd1); c2->NonZeroPattern (ud, v2, d2, dd2);

  }
  */
  virtual void NonZeroPattern (const class ProxyUserData & ud,
                               FlatVector<AutoDiffDiff<1,NonZero>> values) const override
  {
    int dim = Dimension();
    Vector<AutoDiffDiff<1,NonZero>> v1(1), v2(dim);
    c1->NonZeroPattern (ud, v1);
    c2->NonZeroPattern (ud, v2);
    for (size_t j = 0; j < dim; j++)
      values(j) = v1(0) * v2(j);
  }