Method SymMatrixDifferentialOperator

fem/diffop.cpp:1241–1263

  SymMatrixDifferentialOperator ::
  SymMatrixDifferentialOperator (shared_ptr<DifferentialOperator> adiffop, 
                                 int avdim)
  : DifferentialOperator(sqr(avdim)*adiffop->Dim(), adiffop->BlockDim(),
                         adiffop->VB(), adiffop->DiffOrder()),
      diffop(adiffop), vdim(avdim)
  {
    if (adiffop->Dimensions().Size() == 0)
      {
        SetDimensions ( { avdim, avdim } );
        Matrix<> vsembedding(sqr(vdim), vdim*(vdim+1)/2);
        vsembedding = 0.0;
        for (int i = 0, ii = 0; i < vdim; i++)
          for (int j = 0; j <= i; j++, ii++)
            {
              vsembedding(i*vdim+j, ii) = 1;
              vsembedding(j*vdim+i, ii) = 1;
            }
        SetVectorSpaceEmbedding(vsembedding);
      }
    else
      throw Exception("no matrix-valued of vector-valued possible");
  }


  SymMatrixDifferentialOperator :: ~SymMatrixDifferentialOperator ()  { ; }