18 class CeedMassOperator : public mfem::Operator {
20   const mfem::FiniteElementSpace *fes;
31 …static void FESpace2Ceed(const mfem::FiniteElementSpace *fes, const mfem::IntegrationRule &ir, Cee…  in FESpace2Ceed()
33     mfem::Mesh                *mesh  = fes->GetMesh();  in FESpace2Ceed()
34     const mfem::FiniteElement *fe    = fes->GetFE(0);  in FESpace2Ceed()
36     mfem::Array<int>           dof_map;  in FESpace2Ceed()
39         const mfem::H1_SegmentElement *h1_fe = dynamic_cast<const mfem::H1_SegmentElement *>(fe);  in FESpace2Ceed()
45 …const mfem::H1_QuadrilateralElement *h1_fe = dynamic_cast<const mfem::H1_QuadrilateralElement *>(f…  in FESpace2Ceed()
51 …    const mfem::H1_HexahedronElement *h1_fe = dynamic_cast<const mfem::H1_HexahedronElement *>(fe);  in FESpace2Ceed()
57 …const mfem::FiniteElement     *fe1d = fes->FEColl()->FiniteElementForGeometry(mfem::Geometry::SEGM…  in FESpace2Ceed()
58     mfem::DenseMatrix              shape1d(fe1d->GetDof(), ir.GetNPoints());  in FESpace2Ceed()
59     mfem::DenseMatrix              grad_1d(fe1d->GetDof(), ir.GetNPoints());  in FESpace2Ceed()
60     mfem::Vector                   q_ref_1d(ir.GetNPoints()), q_weight_1d(ir.GetNPoints());  in FESpace2Ceed()
61     mfem::Vector                   shape_i(shape1d.Height());  in FESpace2Ceed()
62     mfem::DenseMatrix              grad_i(grad_1d.Height(), 1);  in FESpace2Ceed()
63     const mfem::H1_SegmentElement *h1_fe1d = dynamic_cast<const mfem::H1_SegmentElement *>(fe1d);  in FESpace2Ceed()
65     const mfem::Array<int> &dof_map_1d = h1_fe1d->GetDofMap();  in FESpace2Ceed()
67       const mfem::IntegrationPoint &ip = ir.IntPoint(i);  in FESpace2Ceed()
80     const mfem::Table &el_dof = fes->GetElementToDofTable();  in FESpace2Ceed()
81     mfem::Array<int>   tp_el_dof(el_dof.Size_of_connections());  in FESpace2Ceed()
94 …CeedMassOperator(Ceed ceed, const mfem::FiniteElementSpace *fes) : Operator(fes->GetNDofs()), fes(…  in CeedMassOperator()
95     mfem::Mesh                  *mesh     = fes->GetMesh();  in CeedMassOperator()
98     const mfem::IntegrationRule &ir       = mfem::IntRules.Get(mfem::Geometry::SEGMENT, ir_order);  in CeedMassOperator()
103     const mfem::FiniteElementSpace *mesh_fes = mesh->GetNodalFESpace();  in CeedMassOperator()
173   virtual void Mult(const mfem::Vector &x, mfem::Vector &y) const {  in Mult()