]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Revert changes to petsc_ts_03/04
authorDaniel Arndt <arndtd@ornl.gov>
Tue, 13 Jun 2023 12:20:35 +0000 (08:20 -0400)
committerDaniel Arndt <arndtd@ornl.gov>
Tue, 13 Jun 2023 12:20:35 +0000 (08:20 -0400)
tests/petsc/petsc_ts_03.cc
tests/petsc/petsc_ts_03.output
tests/petsc/petsc_ts_04.cc
tests/petsc/petsc_ts_04.output

index 116c940aa67d5b91a0062e9cfdc4171bf7c3eb9a..0f6722101b50588b3f81eb4eb4dc0653680d6431 100644 (file)
@@ -136,6 +136,7 @@ public:
         // solver is used.
         time_stepper.explicit_function =
           [&](const real_type t, const VectorType &y, VectorType &res) -> void {
+          deallog << "Evaluating explicit function at t=" << t << std::endl;
           res(0) = -kappa * y(0);
           res.compress(VectorOperation::insert);
         };
@@ -150,6 +151,7 @@ public:
                                                  const VectorType &y,
                                                  MatrixType &      A,
                                                  MatrixType &      P) -> void {
+              deallog << "Evaluating explicit Jacobian at t=" << t << std::endl;
               P.set(0, 0, -kappa);
               P.compress(VectorOperation::insert);
             };
index e63ac89b614c225884ab6f63e348556855f99f38..d6ba3e381c640ceaeff392108998b6b80485f6c6 100644 (file)
@@ -39,36 +39,84 @@ DEAL::# Test explicit interface (J 0)
 DEAL::Intermediate output:
 DEAL::  t =0.00000
 DEAL::  y =1.00000  (exact: 1.00000)
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit function at t=0.200000
+DEAL::Evaluating explicit function at t=0.200000
+DEAL::Evaluating explicit function at t=0.200000
+DEAL::Evaluating explicit function at t=0.200000
 DEAL::Intermediate output:
 DEAL::  t =0.200000
 DEAL::  y =0.823864  (exact: 0.818731)
+DEAL::Evaluating explicit function at t=0.400000
+DEAL::Evaluating explicit function at t=0.400000
+DEAL::Evaluating explicit function at t=0.400000
+DEAL::Evaluating explicit function at t=0.400000
 DEAL::Intermediate output:
 DEAL::  t =0.400000
 DEAL::  y =0.675134  (exact: 0.670320)
+DEAL::Evaluating explicit function at t=0.600000
+DEAL::Evaluating explicit function at t=0.600000
+DEAL::Evaluating explicit function at t=0.600000
+DEAL::Evaluating explicit function at t=0.600000
 DEAL::Intermediate output:
 DEAL::  t =0.600000
 DEAL::  y =0.551962  (exact: 0.548812)
+DEAL::Evaluating explicit function at t=0.800000
+DEAL::Evaluating explicit function at t=0.800000
+DEAL::Evaluating explicit function at t=0.800000
+DEAL::Evaluating explicit function at t=0.800000
 DEAL::Intermediate output:
 DEAL::  t =0.800000
 DEAL::  y =0.450798  (exact: 0.449329)
+DEAL::Evaluating explicit function at t=1.00000
+DEAL::Evaluating explicit function at t=1.00000
+DEAL::Evaluating explicit function at t=1.00000
+DEAL::Evaluating explicit function at t=1.00000
 DEAL::Intermediate output:
 DEAL::  t =1.00000
 DEAL::  y =0.368009  (exact: 0.367879)
+DEAL::Evaluating explicit function at t=1.20000
+DEAL::Evaluating explicit function at t=1.20000
+DEAL::Evaluating explicit function at t=1.20000
+DEAL::Evaluating explicit function at t=1.20000
 DEAL::Intermediate output:
 DEAL::  t =1.20000
 DEAL::  y =0.300364  (exact: 0.301194)
+DEAL::Evaluating explicit function at t=1.40000
+DEAL::Evaluating explicit function at t=1.40000
+DEAL::Evaluating explicit function at t=1.40000
+DEAL::Evaluating explicit function at t=1.40000
 DEAL::Intermediate output:
 DEAL::  t =1.40000
 DEAL::  y =0.245132  (exact: 0.246597)
+DEAL::Evaluating explicit function at t=1.60000
+DEAL::Evaluating explicit function at t=1.60000
+DEAL::Evaluating explicit function at t=1.60000
+DEAL::Evaluating explicit function at t=1.60000
 DEAL::Intermediate output:
 DEAL::  t =1.60000
 DEAL::  y =0.200048  (exact: 0.201897)
+DEAL::Evaluating explicit function at t=1.80000
+DEAL::Evaluating explicit function at t=1.80000
+DEAL::Evaluating explicit function at t=1.80000
+DEAL::Evaluating explicit function at t=1.80000
 DEAL::Intermediate output:
 DEAL::  t =1.80000
 DEAL::  y =0.163253  (exact: 0.165299)
+DEAL::Evaluating explicit function at t=2.00000
+DEAL::Evaluating explicit function at t=2.00000
+DEAL::Evaluating explicit function at t=2.00000
+DEAL::Evaluating explicit function at t=2.00000
 DEAL::Intermediate output:
 DEAL::  t =2.00000
 DEAL::  y =0.133225  (exact: 0.135335)
+DEAL::Evaluating explicit function at t=2.20000
+DEAL::Evaluating explicit function at t=2.20000
+DEAL::Evaluating explicit function at t=2.20000
+DEAL::Evaluating explicit function at t=2.20000
 DEAL::Intermediate output:
 DEAL::  t =2.20000
 DEAL::  y =0.108719  (exact: 0.110803)
@@ -163,36 +211,72 @@ DEAL::# Test explicit interface (J 1)
 DEAL::Intermediate output:
 DEAL::  t =0.00000
 DEAL::  y =1.00000  (exact: 1.00000)
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit Jacobian at t=0.100000
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit function at t=0.200000
+DEAL::Evaluating explicit Jacobian at t=0.200000
+DEAL::Evaluating explicit function at t=0.200000
 DEAL::Intermediate output:
 DEAL::  t =0.200000
 DEAL::  y =0.823864  (exact: 0.818731)
+DEAL::Evaluating explicit function at t=0.400000
+DEAL::Evaluating explicit Jacobian at t=0.400000
+DEAL::Evaluating explicit function at t=0.400000
 DEAL::Intermediate output:
 DEAL::  t =0.400000
 DEAL::  y =0.675134  (exact: 0.670320)
+DEAL::Evaluating explicit function at t=0.600000
+DEAL::Evaluating explicit Jacobian at t=0.600000
+DEAL::Evaluating explicit function at t=0.600000
 DEAL::Intermediate output:
 DEAL::  t =0.600000
 DEAL::  y =0.551962  (exact: 0.548812)
+DEAL::Evaluating explicit function at t=0.800000
+DEAL::Evaluating explicit Jacobian at t=0.800000
+DEAL::Evaluating explicit function at t=0.800000
 DEAL::Intermediate output:
 DEAL::  t =0.800000
 DEAL::  y =0.450798  (exact: 0.449329)
+DEAL::Evaluating explicit function at t=1.00000
+DEAL::Evaluating explicit Jacobian at t=1.00000
+DEAL::Evaluating explicit function at t=1.00000
 DEAL::Intermediate output:
 DEAL::  t =1.00000
 DEAL::  y =0.368009  (exact: 0.367879)
+DEAL::Evaluating explicit function at t=1.20000
+DEAL::Evaluating explicit Jacobian at t=1.20000
+DEAL::Evaluating explicit function at t=1.20000
 DEAL::Intermediate output:
 DEAL::  t =1.20000
 DEAL::  y =0.300364  (exact: 0.301194)
+DEAL::Evaluating explicit function at t=1.40000
+DEAL::Evaluating explicit Jacobian at t=1.40000
+DEAL::Evaluating explicit function at t=1.40000
 DEAL::Intermediate output:
 DEAL::  t =1.40000
 DEAL::  y =0.245132  (exact: 0.246597)
+DEAL::Evaluating explicit function at t=1.60000
+DEAL::Evaluating explicit Jacobian at t=1.60000
+DEAL::Evaluating explicit function at t=1.60000
 DEAL::Intermediate output:
 DEAL::  t =1.60000
 DEAL::  y =0.200048  (exact: 0.201897)
+DEAL::Evaluating explicit function at t=1.80000
+DEAL::Evaluating explicit Jacobian at t=1.80000
+DEAL::Evaluating explicit function at t=1.80000
 DEAL::Intermediate output:
 DEAL::  t =1.80000
 DEAL::  y =0.163253  (exact: 0.165299)
+DEAL::Evaluating explicit function at t=2.00000
+DEAL::Evaluating explicit Jacobian at t=2.00000
+DEAL::Evaluating explicit function at t=2.00000
 DEAL::Intermediate output:
 DEAL::  t =2.00000
 DEAL::  y =0.133225  (exact: 0.135335)
+DEAL::Evaluating explicit function at t=2.20000
+DEAL::Evaluating explicit Jacobian at t=2.20000
+DEAL::Evaluating explicit function at t=2.20000
 DEAL::Intermediate output:
 DEAL::  t =2.20000
 DEAL::  y =0.108719  (exact: 0.110803)
index 73c0cda0c8d1fc9e89e31d7586a04733f1a4441d..0d1d2ad155f487be5feaf841d98461cd659150a0 100644 (file)
@@ -151,6 +151,8 @@ public:
         // solver is used.
         time_stepper.explicit_function =
           [&](const real_type t, const VectorType &y, VectorType &res) -> void {
+          deallog << "Evaluating explicit function at t=" << t << std::endl;
+
           if (t > last_eval_time + 0.1)
             {
               deallog << "Time step too large: last_eval_time="
@@ -174,6 +176,7 @@ public:
                                                  const VectorType &y,
                                                  MatrixType &      A,
                                                  MatrixType &      P) -> void {
+              deallog << "Evaluating explicit Jacobian at t=" << t << std::endl;
               P.set(0, 0, -kappa);
               P.compress(VectorOperation::insert);
             };
index b731a7e0e0ef4d167ac36a836163a25d0244b45c..26db1c105ab718ee5ca04db836b02bbd515df95a 100644 (file)
@@ -39,9 +39,18 @@ DEAL::# Test explicit interface (J 0)
 DEAL::Intermediate output:
 DEAL::  t =0.00000
 DEAL::  y =1.00000  (exact: 1.00000)
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit function at t=0.200000
+DEAL::Evaluating explicit function at t=0.200000
+DEAL::Evaluating explicit function at t=0.200000
+DEAL::Evaluating explicit function at t=0.200000
 DEAL::Intermediate output:
 DEAL::  t =0.200000
 DEAL::  y =0.823864  (exact: 0.818731)
+DEAL::Evaluating explicit function at t=0.400000
 DEAL::Time step too large: last_eval_time=0.200000, t=0.400000
 DEAL::Time stepper aborted with an exception:
 DEAL::
@@ -96,9 +105,16 @@ DEAL::# Test explicit interface (J 1)
 DEAL::Intermediate output:
 DEAL::  t =0.00000
 DEAL::  y =1.00000  (exact: 1.00000)
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit Jacobian at t=0.100000
+DEAL::Evaluating explicit function at t=0.100000
+DEAL::Evaluating explicit function at t=0.200000
+DEAL::Evaluating explicit Jacobian at t=0.200000
+DEAL::Evaluating explicit function at t=0.200000
 DEAL::Intermediate output:
 DEAL::  t =0.200000
 DEAL::  y =0.823864  (exact: 0.818731)
+DEAL::Evaluating explicit function at t=0.400000
 DEAL::Time step too large: last_eval_time=0.200000, t=0.400000
 DEAL::Time stepper aborted with an exception:
 DEAL::

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.