On Wed, May 12, 2021 at 12:50 PM edgar <edgar...@cryptolab.net> wrote:
> On 2021-05-12 16:57, edgar wrote:
> > On 2021-05-10 15:59, John Peterson wrote:
> >> On Thu, May 6, 2021 at 9:48 PM edgar <edgar...@cryptolab.net> wrote:
> >> We renumbered our examples once (many years ago) so it's possible that
> >> references to "ex13" that you see have simply never been updated
> >> properly.
> >> If you can point us to them I'll take a look.
> >
> > I see. I found these:
> >
> > ┌────
> > │ find . -type f -exec grep --color -nH --null -e 'example 13' \{\} +
> > └────
> >
> > ┌────
> > │ ./fem_system/ex1/fem_system_ex1.C\026:// example 13 can be solved
> > using the
> > │ ./systems_of_equations/ex3/systems_of_equations_ex3.C\025:// example
> > 13 can be solved using a scalar Lagrange multiplier
> > │ ./vector_fe/ex2/vector_fe_ex2.C\026:// example 13 can be solved using
> > the
> > └────
>
> I found another one of these in fem_system_ex3.C:
>
> // This is just Systems of Equations Example 6 recast.
>
> In that same example, it reads:
>
> // Declare the system "Navier-Stokes" and its variables.
>
> But the system is named "Linear Elasticity". I also think that these two
> can be merged into 1:
>
> #+begin_src diff
> --- a/examples/fem_system/ex3/fem_system_ex3.C 2021-03-22
> 18:33:18.000000000 -0600
> +++ b/examples/fem_system/ex3/fem_system_ex3.C 2021-05-12
> 12:41:14.847399776 -0500
> @@ -172,7 +172,7 @@
> // Create an equation systems object.
> EquationSystems equation_systems (mesh);
>
> - // Declare the system "Navier-Stokes" and its variables.
> + // Declare the system "Linear Elasticity" and its variables.
> ElasticitySystem & system =
> equation_systems.add_system<ElasticitySystem> ("Linear
> Elasticity");
>
> @@ -195,10 +195,10 @@
> a_system->add_variable("u_accel", FIRST, LAGRANGE);
> a_system->add_variable("v_accel", FIRST, LAGRANGE);
> a_system->add_variable("w_accel", FIRST, LAGRANGE);
> - }
>
> - if (time_solver == std::string("newmark"))
> - system.time_solver = libmesh_make_unique<NewmarkSolver>(system);
> + system.time_solver =
> libmesh_make_unique<NewmarkSolver>(system);
> + }
> +
>
> else if( time_solver == std::string("euler") )
> {
> #+end_src
>
Thanks for pointing those out, I pushed a fix for both.
--
John
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
