Method GenerateMatrix

comp/Hddf.cpp:724–736

    template <typename MIP, typename MAT>
    static void GenerateMatrix(const FiniteElement &fel,
                               const MIP &mip,
                               MAT &&mat, LocalHeap &lh)
    {
      Mat<3, 3> jacobian = mip.GetJacobian();
      Cast(fel).CalcShape(mip.IP(), Trans(mat) /* , jacobian*/);
      for (int i = 0; i < fel.GetNDof(); i++)
      {
        FlatMatrix<> shapei(3, 3, &mat(0, i));                                                          // &mat(i, 0)
        Mat<3, 3> trafo_shapei = Cof(mip.GetJacobian()) * Mat<3, 3>(shapei) * Trans(mip.GetJacobian()); // Mat<3, 3> trafo_shapei = (1 / mip.GetJacobiDet()) * mip.GetJacobianInverse() * shapei * mip.GetJacobian(); // check if this is correct
        shapei = trafo_shapei;
      }
    }
  };