Method SymDevMatrixDifferentialOperator

fem/diffop.cpp:1555–1581

  SymDevMatrixDifferentialOperator ::
  SymDevMatrixDifferentialOperator (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-1);
        vsembedding = 0.0;
        for (int i = 0, ii = 0; i < vdim; i++)
          for (int j = 0; j <= i; j++, ii++)
            if (ii < vsembedding.Width())
              {
                vsembedding(i*vdim+j, ii) = 1;
                vsembedding(j*vdim+i, ii) = 1;
              }
        for (int i = 0; i < vdim-1; i++)
          vsembedding.Row(vdim*vdim-1) -= vsembedding.Row(i*(vdim+1));
//        cout << "SymDev vsembedding = " << vsembedding << endl;
        SetVectorSpaceEmbedding(vsembedding);
      }
    else
      throw Exception("no matrix-valued of vector-valued possible");
  }


  SymDevMatrixDifferentialOperator :: ~SymDevMatrixDifferentialOperator ()  { ; }