Method NonZeroPattern

fem/coefficient_matrix.cpp:947–960

    { return Array<shared_ptr<CoefficientFunction>>({ c1 } ); }  

    virtual void NonZeroPattern (const class ProxyUserData & ud,
                                 FlatVector<AutoDiffDiff<1,NonZero>> values) const override
    {
      int hd = Dimensions()[0];
      c1->NonZeroPattern (ud, values);

      for (int i = 0; i < hd; i++)
        for (int j = 0; j < hd; j++)
          {
            int ii = i*hd+j;
            int jj = j*hd+i;
            values(ii) = values(ii)+values(jj);  // logical or
          }
    }
  
    virtual void NonZeroPattern (const class ProxyUserData & ud,
                                 FlatArray<FlatVector<AutoDiffDiff<1,NonZero>>> input,