18 class CeedDiffusionOperator : public mfem::Operator {
20   const mfem::FiniteElementSpace *fes;
32 …static void FESpace2Ceed(const mfem::FiniteElementSpace *fes, const mfem::IntegrationRule &ir, Cee…  in FESpace2Ceed()
34     mfem::Mesh                *mesh  = fes->GetMesh();  in FESpace2Ceed()
35     const mfem::FiniteElement *fe    = fes->GetFE(0);  in FESpace2Ceed()
37     mfem::Array<int>           dof_map;  in FESpace2Ceed()
40         const mfem::H1_SegmentElement *h1_fe = dynamic_cast<const mfem::H1_SegmentElement *>(fe);  in FESpace2Ceed()
46 …const mfem::H1_QuadrilateralElement *h1_fe = dynamic_cast<const mfem::H1_QuadrilateralElement *>(f…  in FESpace2Ceed()
52 …    const mfem::H1_HexahedronElement *h1_fe = dynamic_cast<const mfem::H1_HexahedronElement *>(fe);  in FESpace2Ceed()
58 …const mfem::FiniteElement     *fe1d = fes->FEColl()->FiniteElementForGeometry(mfem::Geometry::SEGM…  in FESpace2Ceed()
59     mfem::DenseMatrix              shape1d(fe1d->GetDof(), ir.GetNPoints());  in FESpace2Ceed()
60     mfem::DenseMatrix              grad_1d(fe1d->GetDof(), ir.GetNPoints());  in FESpace2Ceed()
61     mfem::Vector                   q_ref_1d(ir.GetNPoints()), q_weight_1d(ir.GetNPoints());  in FESpace2Ceed()
62     mfem::Vector                   shape_i(shape1d.Height());  in FESpace2Ceed()
63     mfem::DenseMatrix              grad_i(grad_1d.Height(), 1);  in FESpace2Ceed()
64     const mfem::H1_SegmentElement *h1_fe1d = dynamic_cast<const mfem::H1_SegmentElement *>(fe1d);  in FESpace2Ceed()
66     const mfem::Array<int> &dof_map_1d = h1_fe1d->GetDofMap();  in FESpace2Ceed()
68       const mfem::IntegrationPoint &ip = ir.IntPoint(i);  in FESpace2Ceed()
81     const mfem::Table &el_dof = fes->GetElementToDofTable();  in FESpace2Ceed()
82     mfem::Array<int>   tp_el_dof(el_dof.Size_of_connections());  in FESpace2Ceed()
95 …CeedDiffusionOperator(Ceed ceed, const mfem::FiniteElementSpace *fes) : Operator(fes->GetNDofs()),…  in CeedDiffusionOperator()
96     mfem::Mesh                  *mesh     = fes->GetMesh();  in CeedDiffusionOperator()
99     const mfem::IntegrationRule &ir       = mfem::IntRules.Get(mfem::Geometry::SEGMENT, ir_order);  in CeedDiffusionOperator()
104     const mfem::FiniteElementSpace *mesh_fes = mesh->GetNodalFESpace();  in CeedDiffusionOperator()
175   virtual void Mult(const mfem::Vector &x, mfem::Vector &y) const {  in Mult()