On Wed, Jun 4, 2014 at 10:25 AM, Miguel Angel Salazar de Troya < salazardetr...@gmail.com> wrote:
> Thanks for your responses. I've been digging into the code and I've > realized that if I just wanted to use refinement with LAGRANGE_VEC I would > have to modify the function "lagrange_compute_constraints" to accommodate > it for LAGRANGE_VEC and call it from the corresponding function in > FEInterface::compute_contraints. As for the projection, in the class > ProjectionVector, for Lagrange elements the function coarsened_dof_values > is not called. I would have to modify the ProjectionVector::operator() so > it can handle LAGRANGE_VEC the same way that handles the > fe_type.family LAGRANGE. I don't know if I'm leaving something behind. > Honestly, I'm a bit overwhelmed with this. No kidding. You are trying to jump in at the deep end here for sure. > I read this paper: > > Wolfgang Bangerth, Oliver Kayser-Herold > *Data Structures and Requirements for hp Finite Element Software* > ACM Transactions on Mathematical Software, vol. 36 (2009), pp. 4/1-4/31. > > to try to understand the process better. The final conclusion mentions how > long it took them to write the code and it's overwhelming. I guess that's > how much it took you guys to write the libmesh hp-refinement version. I > wanted to understand what I needed to do to achieve my goals of having > adaptive mesh refinement for LAGRANGE_VEC elements. At first, I would like > to have just h-refinement. Maybe later, once I grasp more concepts and > become acquainted with the code, I will extend it to hp-refinement. > 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... -- John ------------------------------------------------------------------------------ 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 Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
