Re: [Libmesh-users] Example 12 →Example 13

[edgar <edg...@cr...>]

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 <edg...@cr...> 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