Thanks for your response. I am trying to modify the ProjectVector::operator() function to be able to handle the LAGRANGE_VEC. One of the first things this function does is this (line 782 at system_projection.C )
AutoPtr<FEBase> fe (FEBase::build(dim, base_fe_type)) This function does not handle the fe_type.family LAGRANGE_VEC, I would have to call: AutoPtr<FEVectorBase> fe (FEVectorBase::build(dim, base_fe_type)); as John pointed out. How could I handle both cases without a conditional? That's the first thing I tried, but I cannot create a variable inside the conditional. Is there anyway to do this with templates? I don't have much experience with templates (I briefly saw them in a course) and I cannot come up with a solution. Thanks Miguel On Wed, Jun 4, 2014 at 12:25 PM, Roy Stogner <[email protected]> wrote: > > On Wed, 4 Jun 2014, John Peterson wrote: > > p-refinement for Lagrange finite elements is pretty limited in >> libMesh anyway, isn't it? We only go up to p=2 for one thing... >> > > Uniform p-refinement is limited to p=1 or p=2 on 2D/3D Lagrange > elements; Even in 1D I think we only go to p=3. > > Adaptive p-refinement is impossible with Lagrange elements for us > right now. Since the hp constraints code was already looking very > hairy, I simplified slightly by assuming that we only needed it to > work for p-nested bases. > --- > Roy -- *Miguel Angel Salazar de Troya* Graduate Research Assistant Department of Mechanical Science and Engineering University of Illinois at Urbana-Champaign (217) 550-2360 [email protected] ------------------------------------------------------------------------------ Learn Graph Databases - Download FREE O'Reilly Book "Graph Databases" is the definitive new guide to graph databases and their applications. Written by three acclaimed leaders in the field, this first edition is now available. Download your free book today! http://p.sf.net/sfu/NeoTech _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
