Method NonZeroPattern

fem/coefficient.cpp:3452–3465

virtual void NonZeroPattern (const class ProxyUserData & ud, FlatVector nonzero, FlatVector nonzero_deriv, FlatVector nonzero_dderiv) const override { FlatArray hdims = Dimensions(); Vector v1(hdims[0]*inner_dim), v2(inner_dim); Vector d1(hdims[0]*inner_dim), d2(inner_dim); Vector dd1(hdims[0]*inner_dim), dd2

  }
  */
  virtual void NonZeroPattern (const class ProxyUserData & ud,
                               FlatVector<AutoDiffDiff<1,NonZero>> values) const override
  {
    // FlatArray<int> hdims = Dimensions();
    Vector<AutoDiffDiff<1,NonZero>> va(Dimension()*inner_dim), vb(inner_dim);
    c1->NonZeroPattern (ud, va);
    c2->NonZeroPattern (ud, vb);
    
    values = NonZero(false);
    
    for (size_t i = 0; i < Dimension(); i++)
      for (size_t j = 0; j < inner_dim; j++)
        values(i) += va(i*inner_dim+j) * vb(j);
  }
  
  virtual void NonZeroPattern (const class ProxyUserData & ud,
                               FlatArray<FlatVector<AutoDiffDiff<1,NonZero>>> input,