diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index da5825a54e917ee672ef1fd7bec78fe7465cefdf..7a7a56638ba7734559b36640386b61f76e341091 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -114,7 +114,7 @@ void assemble_neumann(GridView const &gridView, FEBasis const &feBasis, LocalVectorType SampleVector(0); // FIXME: random values (time-dependent) SampleVector[0] = - (time <= 0.5) ? sin(time * 2 * M_PI) * 1e5 : (time - 0.5) * 2 * 1e5; + (time <= 0.5) ? sin(time * 2 * M_PI) * 1 : (time - 0.5) * 2 * 1; SampleVector[1] = 0; ConstantFunction<LocalVectorType, LocalVectorType> const fNeumann( SampleVector); diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset index bc48494414d223e5936ade1b94759922c783da3a..0c3082fa4f33e6ff392415b574e1188b0db14391 100644 --- a/src/one-body-sample.parset +++ b/src/one-body-sample.parset @@ -82,7 +82,7 @@ model = Ruina [boundary.friction.state] evolve = true -initial = 10 +initial = 1e-4 [boundary.friction.ruina] # "For rocks, typical values of A and B range from 0.005 to 0.015" @@ -90,7 +90,7 @@ initial = 10 # Continuum mathematical modeling of slip weakening in geological systems a = 0.015 b = 0.005 -L = 1e-4 +L = 1e-6 # The second-order a and b terms in [..] are very important in # determining whether sliding is stable or unstable and (in appli-