On Fri, Aug 25, 2017 at 7:19 PM, John Peterson <jwpeter...@gmail.com> wrote:
>
>
> On Fri, Aug 25, 2017 at 1:02 PM, Renato Poli <rebp...@gmail.com> wrote:
>
>> Hi John,
>>
>> I am still struggling and running out of ideas here.
>>
>
> Are you using FIRST and SECOND, LAGRANGE variables? If so, are you also
> changing the mesh when you switch back and forth between orders?
>
Only second order. The program is straightforward: creates the mesh, solves
and exits. Very similar to introduction_ex3, which I took as example. I
just changed the mesh and added my boundary conditions.
// Create an equation systems object.
EquationSystems equation_systems (lm_mesh);
System & sys = equation_systems.add_system<LinearImplicitSystem>
("Poisson");
sys.add_variable("u", SECOND);
sys.attach_assemble_function (assemble_poisson);
equation_systems.init();
equation_systems.print_info();
sys.solve();
Mesh Information:
elem_dimensions()={2}
spatial_dimension()=2
n_nodes()=216
n_local_nodes()=216
n_elem()=60
n_local_elem()=60
n_active_elem()=60
n_subdomains()=1
n_partitions()=1
n_processors()=1
n_threads()=1
processor_id()=0
EquationSystems
n_systems()=1
System #0, "Poisson"
Type "LinearImplicit"
Variables="u"
Finite Element Types="LAGRANGE"
Approximation Orders="SECOND"
n_dofs()=216
n_local_dofs()=216
n_constrained_dofs()=0
n_local_constrained_dofs()=0
n_vectors()=1
n_matrices()=1
DofMap Sparsity
Average On-Processor Bandwidth <= 8.55556
Average Off-Processor Bandwidth <= 0
Maximum On-Processor Bandwidth <= 33
Maximum Off-Processor Bandwidth <= 0
DofMap Constraints
Number of DoF Constraints = 0
>
> --
> John
>
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
