]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Try fixing symengine tests 17076/head
authorDaniel Arndt <arndtd@ornl.gov>
Fri, 24 May 2024 18:47:05 +0000 (14:47 -0400)
committerDaniel Arndt <arndtd@ornl.gov>
Wed, 5 Jun 2024 12:59:54 +0000 (08:59 -0400)
14 files changed:
tests/symengine/basic_05.with_symengine_with_llvm=true.output
tests/symengine/basic_06.output.1 [new file with mode: 0644]
tests/symengine/basic_07.output.1 [new file with mode: 0644]
tests/symengine/batch_optimizer_01_3a.with_symengine_with_llvm=true.output
tests/symengine/batch_optimizer_01_3b.with_symengine_with_llvm=true.output
tests/symengine/batch_optimizer_02_1a.output.release [new file with mode: 0644]
tests/symengine/batch_optimizer_02_1b.output.debug [new file with mode: 0644]
tests/symengine/batch_optimizer_02_2a.output.release [new file with mode: 0644]
tests/symengine/batch_optimizer_02_2b.output.debug [new file with mode: 0644]
tests/symengine/batch_optimizer_02_3a.with_symengine_with_llvm=true.output
tests/symengine/batch_optimizer_02_3b.with_symengine_with_llvm=true.output
tests/symengine/batch_optimizer_03_3a.with_symengine_with_llvm=true.output
tests/symengine/batch_optimizer_03_3b.with_symengine_with_llvm=true.output
tests/symengine/step-44-sd-quadrature_level_03.with_symengine_with_llvm=true.output

index 14de91cc39e6b3a7abda566072f513d62dfb5b1b..e01b1824ac676d1f821ef6ce022f47d8fddfb605 100644 (file)
@@ -1,4 +1,5 @@
 
+DEAL::2.0*s**2.2*(v_0**2 + v_1**2 + v_2**2)**3*(1.0*t_21*t_12 - 1.0*t_22*t_11)*(-1.0*t_20*(1.0*t_01*t_12 - 1.0*t_02*t_11) + 1.0*t_21*(1.0*t_00*t_12 - 1.0*t_02*t_10) - 1.0*t_22*(1.0*t_00*t_11 - 1.0*t_01*t_10))
 DEAL::Value (subs): 2.91338e+13
 DEAL::Value (lambda): 2.91338e+13
 DEAL::Value (llvm): 2.91338e+13
diff --git a/tests/symengine/basic_06.output.1 b/tests/symengine/basic_06.output.1
new file mode 100644 (file)
index 0000000..f064c42
--- /dev/null
@@ -0,0 +1,13 @@
+
+DEAL::C,  Qi,  Q(C)
+DEAL::f_CQ_symb: 0.5*C**2*Qi
+DEAL::df_CQ_dC_symb: 1.0*C*Qi
+DEAL::df_CQ_dC_symb_subs: 1.0*C*Q(C)
+DEAL::D2f_CQ_DC_dC_symb: 1.0*C*Derivative(Q(C), C) + 1.0*Q(C)
+DEAL::f_CQ_subs: 37.5
+DEAL::f_CQ_subs: eval: 37.5000
+DEAL::df_CQ_dC: subs: 15.0
+DEAL::df_CQ_dC: eval: 15.0000
+DEAL::D2f_CQ_DC_dC: subs 1: 3.0 + 0.6*C
+DEAL::D2f_CQ_DC_dC: subs 2: 6.0
+DEAL::D2f_CQ_DC_dC: eval: 6.00000
diff --git a/tests/symengine/basic_07.output.1 b/tests/symengine/basic_07.output.1
new file mode 100644 (file)
index 0000000..5c54cd6
--- /dev/null
@@ -0,0 +1,13 @@
+
+DEAL::C,  Qi,  Q(C)
+DEAL::f_CQ_symb: 0.5*C**2*Qi
+DEAL::df_CQ_dC_symb: 1.0*C*Qi
+DEAL::df_CQ_dC_symb_subs: 1.0*C*Q(C)
+DEAL::D2f_CQ_DC_dC_symb: 1.0*C*Derivative(Q(C), C) + 1.0*Q(C)
+DEAL::f_CQ_subs: 187.5
+DEAL::f_CQ_subs: eval: 187.500
+DEAL::df_CQ_dC: subs: 75.0
+DEAL::df_CQ_dC: eval: 75.0000
+DEAL::D2f_CQ_DC_dC: subs 1: 6.0*C + 1.0*Qi
+DEAL::D2f_CQ_DC_dC: subs 2: 45.0
+DEAL::D2f_CQ_DC_dC: eval: 45.0000
index a53b5ea67a5fd23530a14fae19ab6f581b66ff67..6efc1bda3a9222e0bb0410609fb1e8255e3647bf 100644 (file)
@@ -1,4 +1,7 @@
 
+DEAL:Float::Optimisation method: llvm
+Optimisation flags:  OptimizationFlags|-O2|
+Skip non-implemented functions: 0
 DEAL:Float:Add:Symbolic computation: Function initialisation::symb_s0: s1
 DEAL:Float:Add:Symbolic computation: Function initialisation::symb_s1: s2
 DEAL:Float:Add:Symbolic computation: Function initialisation::symb_psi: (s1 + s2)**2
@@ -8,6 +11,8 @@ DEAL:Float:Add:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2
 DEAL:Float:Add:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2
 DEAL:Float:Add:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2
 DEAL:Float:Add:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2
+DEAL:Float:Add:Symbolic computation: Optimization::Add: About to optimise... done.
+DEAL:Float:Add:Symbolic computation: Substitution::Add: About to substitute... done.
 DEAL:Float:Add:Symbolic computation: Substitution::subs_psi: 25.00000000
 DEAL:Float:Add:Symbolic computation: Substitution::subs_dpsi_ds0: 10.00000000
 DEAL:Float:Add:Symbolic computation: Substitution::subs_dpsi_ds1: 10.00000000
@@ -24,6 +29,8 @@ DEAL:Float:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2
 DEAL:Float:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -2
 DEAL:Float:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -2
 DEAL:Float:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2
+DEAL:Float:Subtract:Symbolic computation: Optimization::Subtract: About to optimise... done.
+DEAL:Float:Subtract:Symbolic computation: Substitution::Subtract: About to substitute... done.
 DEAL:Float:Subtract:Symbolic computation: Substitution::subs_psi: 1.000000000
 DEAL:Float:Subtract:Symbolic computation: Substitution::subs_dpsi_ds0: 2.000000000
 DEAL:Float:Subtract:Symbolic computation: Substitution::subs_dpsi_ds1: -2.000000000
@@ -40,6 +47,8 @@ DEAL:Float:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2
 DEAL:Float:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 4*s1*s2
 DEAL:Float:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 4*s1*s2
 DEAL:Float:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2*s1**2
+DEAL:Float:Multiply:Symbolic computation: Optimization::Multiply: About to optimise... done.
+DEAL:Float:Multiply:Symbolic computation: Substitution::Multiply: About to substitute... done.
 DEAL:Float:Multiply:Symbolic computation: Substitution::subs_psi: 36.00000000
 DEAL:Float:Multiply:Symbolic computation: Substitution::subs_dpsi_ds0: 24.00000000
 DEAL:Float:Multiply:Symbolic computation: Substitution::subs_dpsi_ds1: 36.00000000
@@ -56,6 +65,8 @@ DEAL:Float:Divide:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2/s
 DEAL:Float:Divide:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -4*s1/s2**3
 DEAL:Float:Divide:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -4*s1/s2**3
 DEAL:Float:Divide:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 6*s1**2/s2**4
+DEAL:Float:Divide:Symbolic computation: Optimization::Divide: About to optimise... done.
+DEAL:Float:Divide:Symbolic computation: Substitution::Divide: About to substitute... done.
 DEAL:Float:Divide:Symbolic computation: Substitution::subs_psi: 2.250000000
 DEAL:Float:Divide:Symbolic computation: Substitution::subs_dpsi_ds0: 1.500000000
 DEAL:Float:Divide:Symbolic computation: Substitution::subs_dpsi_ds1: -2.250000000
@@ -72,6 +83,8 @@ DEAL:Float:Power:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2.0*
 DEAL:Float:Power:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2.0*s1**(-1 + s2)*s2*log(s1) + 2.0*s1**(-1 + s2)
 DEAL:Float:Power:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2.0*s1**(-1 + s2)*s2*log(s1) + 2.0*s1**(-1 + s2)
 DEAL:Float:Power:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2.0*s1**s2*log(s1)**2
+DEAL:Float:Power:Symbolic computation: Optimization::Power: About to optimise... done.
+DEAL:Float:Power:Symbolic computation: Substitution::Power: About to substitute... done.
 DEAL:Float:Power:Symbolic computation: Substitution::subs_psi: 18.00000000
 DEAL:Float:Power:Symbolic computation: Substitution::subs_dpsi_ds0: 12.00000000
 DEAL:Float:Power:Symbolic computation: Substitution::subs_dpsi_ds1: 19.77502060
@@ -88,12 +101,14 @@ DEAL:Float:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -0.5*
 DEAL:Float:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -0.5*s1*s2/(s1*s2)**(3/2) + 1.0*(s1*s2)**(-1/2)
 DEAL:Float:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -0.5*s1*s2/(s1*s2)**(3/2) + 1.0*(s1*s2)**(-1/2)
 DEAL:Float:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -0.5*s1**2/(s1*s2)**(3/2)
+DEAL:Float:Sqrt:Symbolic computation: Optimization::Sqrt: About to optimise... done.
+DEAL:Float:Sqrt:Symbolic computation: Substitution::Sqrt: About to substitute... done.
 DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_psi: 4.898979664
-DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_dpsi_ds0: 0.8164965510
-DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_dpsi_ds1: 1.224744797
+DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_dpsi_ds0: 0.8164966106
+DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_dpsi_ds1: 1.224744916
 DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_d2psi_ds0_ds0: -0.1360827684
-DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_d2psi_ds1_ds0: 0.2041241229
-DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_d2psi_ds0_ds1: 0.2041241229
+DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_d2psi_ds1_ds0: 0.2041241527
+DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_d2psi_ds0_ds1: 0.2041241527
 DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_d2psi_ds1_ds1: -0.3061862290
 DEAL:Float:Exponential:Symbolic computation: Function initialisation::symb_s0: s1
 DEAL:Float:Exponential:Symbolic computation: Function initialisation::symb_s1: s2
@@ -104,6 +119,8 @@ DEAL:Float:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0
 DEAL:Float:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2*exp(s2)*s1
 DEAL:Float:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2*exp(s2)*s1
 DEAL:Float:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: exp(s2)*s1**2
+DEAL:Float:Exponential:Symbolic computation: Optimization::Exponential: About to optimise... done.
+DEAL:Float:Exponential:Symbolic computation: Substitution::Exponential: About to substitute... done.
 DEAL:Float:Exponential:Symbolic computation: Substitution::subs_psi: 66.50150299
 DEAL:Float:Exponential:Symbolic computation: Substitution::subs_dpsi_ds0: 44.33433533
 DEAL:Float:Exponential:Symbolic computation: Substitution::subs_dpsi_ds1: 66.50150299
@@ -120,6 +137,8 @@ DEAL:Float:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0:
 DEAL:Float:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2*s1/s2
 DEAL:Float:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2*s1/s2
 DEAL:Float:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -s1**2/s2**2
+DEAL:Float:Logarithm:Symbolic computation: Optimization::Logarithm: About to optimise... done.
+DEAL:Float:Logarithm:Symbolic computation: Substitution::Logarithm: About to substitute... done.
 DEAL:Float:Logarithm:Symbolic computation: Substitution::subs_psi: 6.238324642
 DEAL:Float:Logarithm:Symbolic computation: Substitution::subs_dpsi_ds0: 4.158883095
 DEAL:Float:Logarithm:Symbolic computation: Substitution::subs_dpsi_ds1: 4.500000000
@@ -136,6 +155,8 @@ DEAL:Float:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi
 DEAL:Float:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1/(s1*s2*log(s1)**2)
 DEAL:Float:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1/(s1*s2*log(s1)**2)
 DEAL:Float:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1/(s2**2*log(s1))
+DEAL:Float:Logarithm with Base:Symbolic computation: Optimization::Logarithm with Base: About to optimise... done.
+DEAL:Float:Logarithm with Base:Symbolic computation: Substitution::Logarithm with Base: About to substitute... done.
 DEAL:Float:Logarithm with Base:Symbolic computation: Substitution::subs_psi: 0.6309297681
 DEAL:Float:Logarithm with Base:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1914323419
 DEAL:Float:Logarithm with Base:Symbolic computation: Substitution::subs_dpsi_ds1: 0.4551196098
@@ -152,6 +173,8 @@ DEAL:Float:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_d
 DEAL:Float:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 0.868588963806504*s1/s2
 DEAL:Float:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 0.868588963806504*s1/s2
 DEAL:Float:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -0.434294481903252*s1**2/s2**2
+DEAL:Float:Logarithm base 10:Symbolic computation: Optimization::Logarithm base 10: About to optimise... done.
+DEAL:Float:Logarithm base 10:Symbolic computation: Substitution::Logarithm base 10: About to substitute... done.
 DEAL:Float:Logarithm base 10:Symbolic computation: Substitution::subs_psi: 2.709270000
 DEAL:Float:Logarithm base 10:Symbolic computation: Substitution::subs_dpsi_ds0: 1.806180120
 DEAL:Float:Logarithm base 10:Symbolic computation: Substitution::subs_dpsi_ds1: 1.954325199
@@ -168,6 +191,8 @@ DEAL:Float:Sin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s2**2
 DEAL:Float:Sin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -cos(s2/s1)/s1**2 + s2*sin(s2/s1)/s1**3
 DEAL:Float:Sin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -cos(s2/s1)/s1**2 + s2*sin(s2/s1)/s1**3
 DEAL:Float:Sin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -sin(s2/s1)/s1**2
+DEAL:Float:Sin:Symbolic computation: Optimization::Sin: About to optimise... done.
+DEAL:Float:Sin:Symbolic computation: Substitution::Sin: About to substitute... done.
 DEAL:Float:Sin:Symbolic computation: Substitution::subs_psi: 0.6183698177
 DEAL:Float:Sin:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1746416092
 DEAL:Float:Sin:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2619624138
@@ -184,6 +209,8 @@ DEAL:Float:Cos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s2**2
 DEAL:Float:Cos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: sin(s2/s1)/s1**2 + s2*cos(s2/s1)/s1**3
 DEAL:Float:Cos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: sin(s2/s1)/s1**2 + s2*cos(s2/s1)/s1**3
 DEAL:Float:Cos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -cos(s2/s1)/s1**2
+DEAL:Float:Cos:Symbolic computation: Optimization::Cos: About to optimise... done.
+DEAL:Float:Cos:Symbolic computation: Substitution::Cos: About to substitute... done.
 DEAL:Float:Cos:Symbolic computation: Substitution::subs_psi: 0.7858872414
 DEAL:Float:Cos:Symbolic computation: Substitution::subs_dpsi_ds0: 0.1374155134
 DEAL:Float:Cos:Symbolic computation: Substitution::subs_dpsi_ds1: -0.2061232775
@@ -200,6 +227,8 @@ DEAL:Float:Tan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*s2*(
 DEAL:Float:Tan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -(1 + tan(s2/s1)**2)/s1**2 - 2*s2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**3
 DEAL:Float:Tan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -(1 + tan(s2/s1)**2)/s1**2 - 2*s2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**3
 DEAL:Float:Tan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**2
+DEAL:Float:Tan:Symbolic computation: Optimization::Tan: About to optimise... done.
+DEAL:Float:Tan:Symbolic computation: Substitution::Tan: About to substitute... done.
 DEAL:Float:Tan:Symbolic computation: Substitution::subs_psi: 0.7868429422
 DEAL:Float:Tan:Symbolic computation: Substitution::subs_dpsi_ds0: -0.3598048389
 DEAL:Float:Tan:Symbolic computation: Substitution::subs_dpsi_ds1: 0.5397073030
@@ -216,6 +245,8 @@ DEAL:Float:ASin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 1.0*1
 DEAL:Float:ASin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ASin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ASin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 1.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 2.0*1/(s1*s2**3*sqrt(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Float:ASin:Symbolic computation: Optimization::ASin: About to optimise... done.
+DEAL:Float:ASin:Symbolic computation: Substitution::ASin: About to substitute... done.
 DEAL:Float:ASin:Symbolic computation: Substitution::subs_psi: 0.1674480885
 DEAL:Float:ASin:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05634361506
 DEAL:Float:ASin:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08451542258
@@ -232,6 +263,8 @@ DEAL:Float:ACos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -1.0*
 DEAL:Float:ACos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ACos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ACos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 2.0*1/(s1*s2**3*sqrt(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Float:ACos:Symbolic computation: Optimization::ACos: About to optimise... done.
+DEAL:Float:ACos:Symbolic computation: Substitution::ACos: About to substitute... done.
 DEAL:Float:ACos:Symbolic computation: Substitution::subs_psi: 1.403348207
 DEAL:Float:ACos:Symbolic computation: Substitution::subs_dpsi_ds0: 0.05634361506
 DEAL:Float:ACos:Symbolic computation: Substitution::subs_dpsi_ds1: 0.08451542258
@@ -248,7 +281,9 @@ DEAL:Float:ATan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -2.0*
 DEAL:Float:ATan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -2.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ATan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -2.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ATan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -2.0*1/(s1**3*s2**5*(1 + 1.0*1/(s1**2*s2**2))**2) + 2.0*1/(s1*s2**3*(1 + 1.0*1/(s1**2*s2**2)))
-DEAL:Float:ATan:Symbolic computation: Substitution::subs_psi: 0.1651486754
+DEAL:Float:ATan:Symbolic computation: Optimization::ATan: About to optimise... done.
+DEAL:Float:ATan:Symbolic computation: Substitution::ATan: About to substitute... done.
+DEAL:Float:ATan:Symbolic computation: Substitution::subs_psi: 0.1651486903
 DEAL:Float:ATan:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05405405536
 DEAL:Float:ATan:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08108108491
 DEAL:Float:ATan:Symbolic computation: Substitution::subs_d2psi_ds0_ds0: 0.03506208956
@@ -264,6 +299,8 @@ DEAL:Float:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 6.0*
 DEAL:Float:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 6.0*s2**2/(s1**2 + s2**2)**2 - 3.0*(s1**2 + s2**2)**(-1)
 DEAL:Float:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -6.0*s1**2/(s1**2 + s2**2)**2 + 3.0*(s1**2 + s2**2)**(-1)
 DEAL:Float:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -6.0*s1*s2/(s1**2 + s2**2)**2
+DEAL:Float:ATan2:Symbolic computation: Optimization::ATan2: About to optimise... done.
+DEAL:Float:ATan2:Symbolic computation: Substitution::ATan2: About to substitute... done.
 DEAL:Float:ATan2:Symbolic computation: Substitution::subs_psi: 1.764007807
 DEAL:Float:ATan2:Symbolic computation: Substitution::subs_dpsi_ds0: -0.4615384936
 DEAL:Float:ATan2:Symbolic computation: Substitution::subs_dpsi_ds1: 0.6923077106
@@ -280,6 +317,8 @@ DEAL:Float:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: s2**2
 DEAL:Float:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -cosh(s2/s1)/s1**2 - s2*sinh(s2/s1)/s1**3
 DEAL:Float:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -cosh(s2/s1)/s1**2 - s2*sinh(s2/s1)/s1**3
 DEAL:Float:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: sinh(s2/s1)/s1**2
+DEAL:Float:Sinh:Symbolic computation: Optimization::Sinh: About to optimise... done.
+DEAL:Float:Sinh:Symbolic computation: Substitution::Sinh: About to substitute... done.
 DEAL:Float:Sinh:Symbolic computation: Substitution::subs_psi: 0.7171584964
 DEAL:Float:Sinh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.2734612525
 DEAL:Float:Sinh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.4101918638
@@ -296,6 +335,8 @@ DEAL:Float:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: s2**2
 DEAL:Float:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -sinh(s2/s1)/s1**2 - s2*cosh(s2/s1)/s1**3
 DEAL:Float:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -sinh(s2/s1)/s1**2 - s2*cosh(s2/s1)/s1**3
 DEAL:Float:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: cosh(s2/s1)/s1**2
+DEAL:Float:Cosh:Symbolic computation: Optimization::Cosh: About to optimise... done.
+DEAL:Float:Cosh:Symbolic computation: Substitution::Cosh: About to substitute... done.
 DEAL:Float:Cosh:Symbolic computation: Substitution::subs_psi: 1.230575562
 DEAL:Float:Cosh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1593685597
 DEAL:Float:Cosh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2390528321
@@ -312,6 +353,8 @@ DEAL:Float:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*s2*
 DEAL:Float:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -(1 - tanh(s2/s1)**2)/s1**2 + 2*s2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**3
 DEAL:Float:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -(1 - tanh(s2/s1)**2)/s1**2 + 2*s2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**3
 DEAL:Float:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**2
+DEAL:Float:Tanh:Symbolic computation: Optimization::Tanh: About to optimise... done.
+DEAL:Float:Tanh:Symbolic computation: Substitution::Tanh: About to substitute... done.
 DEAL:Float:Tanh:Symbolic computation: Substitution::subs_psi: 0.5827829838
 DEAL:Float:Tanh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1467475593
 DEAL:Float:Tanh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2201213539
@@ -328,6 +371,8 @@ DEAL:Float:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -1.0
 DEAL:Float:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1.0*1/(s1**3*s2**5*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 2.0*1/(s1*s2**3*sqrt(1 + 1.0*1/(s1**2*s2**2)))
+DEAL:Float:ASinh:Symbolic computation: Optimization::ASinh: About to optimise... done.
+DEAL:Float:ASinh:Symbolic computation: Substitution::ASinh: About to substitute... done.
 DEAL:Float:ASinh:Symbolic computation: Substitution::subs_psi: 0.1659045517
 DEAL:Float:ASinh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05479966477
 DEAL:Float:ASinh:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08219949901
@@ -344,6 +389,8 @@ DEAL:Float:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s1*
 DEAL:Float:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -s1**2*s2**2/(-1 + s1**2*s2**2)**(3/2) + (-1 + s1**2*s2**2)**(-1/2)
 DEAL:Float:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -s1**2*s2**2/(-1 + s1**2*s2**2)**(3/2) + (-1 + s1**2*s2**2)**(-1/2)
 DEAL:Float:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -s1**3*s2/(-1 + s1**2*s2**2)**(3/2)
+DEAL:Float:ACosh:Symbolic computation: Optimization::ACosh: About to optimise... done.
+DEAL:Float:ACosh:Symbolic computation: Substitution::ACosh: About to substitute... done.
 DEAL:Float:ACosh:Symbolic computation: Substitution::subs_psi: 2.477888823
 DEAL:Float:ACosh:Symbolic computation: Substitution::subs_dpsi_ds0: 0.3380616903
 DEAL:Float:ACosh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.5070925355
@@ -360,6 +407,8 @@ DEAL:Float:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2.0*
 DEAL:Float:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**2) + 2.0*1/(s1*s2**3*(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Float:ATanh:Symbolic computation: Optimization::ATanh: About to optimise... done.
+DEAL:Float:ATanh:Symbolic computation: Substitution::ATanh: About to substitute... done.
 DEAL:Float:ATanh:Symbolic computation: Substitution::subs_psi: 0.1682361215
 DEAL:Float:ATanh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05714286119
 DEAL:Float:ATanh:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08571429551
@@ -376,6 +425,8 @@ DEAL:Float:Erf:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*erf(
 DEAL:Float:Erf:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Float:Erf:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Float:Erf:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -4*exp(-s2**2)*s1**2*s2/sqrt(pi)
+DEAL:Float:Erf:Symbolic computation: Optimization::Erf: About to optimise... done.
+DEAL:Float:Erf:Symbolic computation: Substitution::Erf: About to substitute... done.
 DEAL:Float:Erf:Symbolic computation: Substitution::subs_psi: 8.957901001
 DEAL:Float:Erf:Symbolic computation: Substitution::subs_dpsi_ds0: 5.971933842
 DEAL:Float:Erf:Symbolic computation: Substitution::subs_dpsi_ds1: 0.1860028654
@@ -392,6 +443,8 @@ DEAL:Float:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*erf
 DEAL:Float:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Float:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Float:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 4*exp(-s2**2)*s1**2*s2/sqrt(pi)
+DEAL:Float:Erfc:Symbolic computation: Optimization::Erfc: About to optimise... done.
+DEAL:Float:Erfc:Symbolic computation: Substitution::Erfc: About to substitute... done.
 DEAL:Float:Erfc:Symbolic computation: Substitution::subs_psi: 0.04209961370
 DEAL:Float:Erfc:Symbolic computation: Substitution::subs_dpsi_ds0: 0.02806640789
 DEAL:Float:Erfc:Symbolic computation: Substitution::subs_dpsi_ds1: -0.1860028654
@@ -400,6 +453,9 @@ DEAL:Float:Erfc:Symbolic computation: Substitution::subs_d2psi_ds1_ds0: -0.12400
 DEAL:Float:Erfc:Symbolic computation: Substitution::subs_d2psi_ds0_ds1: -0.1240019053
 DEAL:Float:Erfc:Symbolic computation: Substitution::subs_d2psi_ds1_ds1: 0.7440114617
 DEAL:Float::OK
+DEAL:Double::Optimisation method: llvm
+Optimisation flags:  OptimizationFlags|-O2|
+Skip non-implemented functions: 0
 DEAL:Double:Add:Symbolic computation: Function initialisation::symb_s0: s1
 DEAL:Double:Add:Symbolic computation: Function initialisation::symb_s1: s2
 DEAL:Double:Add:Symbolic computation: Function initialisation::symb_psi: (s1 + s2)**2
@@ -409,6 +465,8 @@ DEAL:Double:Add:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2
 DEAL:Double:Add:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2
 DEAL:Double:Add:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2
 DEAL:Double:Add:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2
+DEAL:Double:Add:Symbolic computation: Optimization::Add: About to optimise... done.
+DEAL:Double:Add:Symbolic computation: Substitution::Add: About to substitute... done.
 DEAL:Double:Add:Symbolic computation: Substitution::subs_psi: 25.00000000
 DEAL:Double:Add:Symbolic computation: Substitution::subs_dpsi_ds0: 10.00000000
 DEAL:Double:Add:Symbolic computation: Substitution::subs_dpsi_ds1: 10.00000000
@@ -425,6 +483,8 @@ DEAL:Double:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0:
 DEAL:Double:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -2
 DEAL:Double:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -2
 DEAL:Double:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2
+DEAL:Double:Subtract:Symbolic computation: Optimization::Subtract: About to optimise... done.
+DEAL:Double:Subtract:Symbolic computation: Substitution::Subtract: About to substitute... done.
 DEAL:Double:Subtract:Symbolic computation: Substitution::subs_psi: 1.000000000
 DEAL:Double:Subtract:Symbolic computation: Substitution::subs_dpsi_ds0: 2.000000000
 DEAL:Double:Subtract:Symbolic computation: Substitution::subs_dpsi_ds1: -2.000000000
@@ -441,6 +501,8 @@ DEAL:Double:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0:
 DEAL:Double:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 4*s1*s2
 DEAL:Double:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 4*s1*s2
 DEAL:Double:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2*s1**2
+DEAL:Double:Multiply:Symbolic computation: Optimization::Multiply: About to optimise... done.
+DEAL:Double:Multiply:Symbolic computation: Substitution::Multiply: About to substitute... done.
 DEAL:Double:Multiply:Symbolic computation: Substitution::subs_psi: 36.00000000
 DEAL:Double:Multiply:Symbolic computation: Substitution::subs_dpsi_ds0: 24.00000000
 DEAL:Double:Multiply:Symbolic computation: Substitution::subs_dpsi_ds1: 36.00000000
@@ -457,6 +519,8 @@ DEAL:Double:Divide:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2/
 DEAL:Double:Divide:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -4*s1/s2**3
 DEAL:Double:Divide:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -4*s1/s2**3
 DEAL:Double:Divide:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 6*s1**2/s2**4
+DEAL:Double:Divide:Symbolic computation: Optimization::Divide: About to optimise... done.
+DEAL:Double:Divide:Symbolic computation: Substitution::Divide: About to substitute... done.
 DEAL:Double:Divide:Symbolic computation: Substitution::subs_psi: 2.250000000
 DEAL:Double:Divide:Symbolic computation: Substitution::subs_dpsi_ds0: 1.500000000
 DEAL:Double:Divide:Symbolic computation: Substitution::subs_dpsi_ds1: -2.250000000
@@ -473,6 +537,8 @@ DEAL:Double:Power:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2.0
 DEAL:Double:Power:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2.0*s1**(-1 + s2)*s2*log(s1) + 2.0*s1**(-1 + s2)
 DEAL:Double:Power:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2.0*s1**(-1 + s2)*s2*log(s1) + 2.0*s1**(-1 + s2)
 DEAL:Double:Power:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2.0*s1**s2*log(s1)**2
+DEAL:Double:Power:Symbolic computation: Optimization::Power: About to optimise... done.
+DEAL:Double:Power:Symbolic computation: Substitution::Power: About to substitute... done.
 DEAL:Double:Power:Symbolic computation: Substitution::subs_psi: 18.00000000
 DEAL:Double:Power:Symbolic computation: Substitution::subs_dpsi_ds0: 12.00000000
 DEAL:Double:Power:Symbolic computation: Substitution::subs_dpsi_ds1: 19.77502120
@@ -489,6 +555,8 @@ DEAL:Double:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -0.5
 DEAL:Double:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -0.5*s1*s2/(s1*s2)**(3/2) + 1.0*(s1*s2)**(-1/2)
 DEAL:Double:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -0.5*s1*s2/(s1*s2)**(3/2) + 1.0*(s1*s2)**(-1/2)
 DEAL:Double:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -0.5*s1**2/(s1*s2)**(3/2)
+DEAL:Double:Sqrt:Symbolic computation: Optimization::Sqrt: About to optimise... done.
+DEAL:Double:Sqrt:Symbolic computation: Substitution::Sqrt: About to substitute... done.
 DEAL:Double:Sqrt:Symbolic computation: Substitution::subs_psi: 4.898979486
 DEAL:Double:Sqrt:Symbolic computation: Substitution::subs_dpsi_ds0: 0.8164965809
 DEAL:Double:Sqrt:Symbolic computation: Substitution::subs_dpsi_ds1: 1.224744871
@@ -505,6 +573,8 @@ DEAL:Double:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds0_ds
 DEAL:Double:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2*exp(s2)*s1
 DEAL:Double:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2*exp(s2)*s1
 DEAL:Double:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: exp(s2)*s1**2
+DEAL:Double:Exponential:Symbolic computation: Optimization::Exponential: About to optimise... done.
+DEAL:Double:Exponential:Symbolic computation: Substitution::Exponential: About to substitute... done.
 DEAL:Double:Exponential:Symbolic computation: Substitution::subs_psi: 66.50150489
 DEAL:Double:Exponential:Symbolic computation: Substitution::subs_dpsi_ds0: 44.33433659
 DEAL:Double:Exponential:Symbolic computation: Substitution::subs_dpsi_ds1: 66.50150489
@@ -521,6 +591,8 @@ DEAL:Double:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0:
 DEAL:Double:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2*s1/s2
 DEAL:Double:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2*s1/s2
 DEAL:Double:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -s1**2/s2**2
+DEAL:Double:Logarithm:Symbolic computation: Optimization::Logarithm: About to optimise... done.
+DEAL:Double:Logarithm:Symbolic computation: Substitution::Logarithm: About to substitute... done.
 DEAL:Double:Logarithm:Symbolic computation: Substitution::subs_psi: 6.238324625
 DEAL:Double:Logarithm:Symbolic computation: Substitution::subs_dpsi_ds0: 4.158883083
 DEAL:Double:Logarithm:Symbolic computation: Substitution::subs_dpsi_ds1: 4.500000000
@@ -537,6 +609,8 @@ DEAL:Double:Logarithm with Base:Symbolic computation: Differentiation::symb_d2ps
 DEAL:Double:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1/(s1*s2*log(s1)**2)
 DEAL:Double:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1/(s1*s2*log(s1)**2)
 DEAL:Double:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1/(s2**2*log(s1))
+DEAL:Double:Logarithm with Base:Symbolic computation: Optimization::Logarithm with Base: About to optimise... done.
+DEAL:Double:Logarithm with Base:Symbolic computation: Substitution::Logarithm with Base: About to substitute... done.
 DEAL:Double:Logarithm with Base:Symbolic computation: Substitution::subs_psi: 0.6309297536
 DEAL:Double:Logarithm with Base:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1914323370
 DEAL:Double:Logarithm with Base:Symbolic computation: Substitution::subs_dpsi_ds1: 0.4551196133
@@ -553,6 +627,8 @@ DEAL:Double:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_
 DEAL:Double:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 0.868588963806504*s1/s2
 DEAL:Double:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 0.868588963806504*s1/s2
 DEAL:Double:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -0.434294481903252*s1**2/s2**2
+DEAL:Double:Logarithm base 10:Symbolic computation: Optimization::Logarithm base 10: About to optimise... done.
+DEAL:Double:Logarithm base 10:Symbolic computation: Substitution::Logarithm base 10: About to substitute... done.
 DEAL:Double:Logarithm base 10:Symbolic computation: Substitution::subs_psi: 2.709269961
 DEAL:Double:Logarithm base 10:Symbolic computation: Substitution::subs_dpsi_ds0: 1.806179974
 DEAL:Double:Logarithm base 10:Symbolic computation: Substitution::subs_dpsi_ds1: 1.954325169
@@ -569,6 +645,8 @@ DEAL:Double:Sin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s2**
 DEAL:Double:Sin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -cos(s2/s1)/s1**2 + s2*sin(s2/s1)/s1**3
 DEAL:Double:Sin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -cos(s2/s1)/s1**2 + s2*sin(s2/s1)/s1**3
 DEAL:Double:Sin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -sin(s2/s1)/s1**2
+DEAL:Double:Sin:Symbolic computation: Optimization::Sin: About to optimise... done.
+DEAL:Double:Sin:Symbolic computation: Substitution::Sin: About to substitute... done.
 DEAL:Double:Sin:Symbolic computation: Substitution::subs_psi: 0.6183698031
 DEAL:Double:Sin:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1746416135
 DEAL:Double:Sin:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2619624203
@@ -585,6 +663,8 @@ DEAL:Double:Cos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s2**
 DEAL:Double:Cos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: sin(s2/s1)/s1**2 + s2*cos(s2/s1)/s1**3
 DEAL:Double:Cos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: sin(s2/s1)/s1**2 + s2*cos(s2/s1)/s1**3
 DEAL:Double:Cos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -cos(s2/s1)/s1**2
+DEAL:Double:Cos:Symbolic computation: Optimization::Cos: About to optimise... done.
+DEAL:Double:Cos:Symbolic computation: Substitution::Cos: About to substitute... done.
 DEAL:Double:Cos:Symbolic computation: Substitution::subs_psi: 0.7858872608
 DEAL:Double:Cos:Symbolic computation: Substitution::subs_dpsi_ds0: 0.1374155118
 DEAL:Double:Cos:Symbolic computation: Substitution::subs_dpsi_ds1: -0.2061232677
@@ -601,6 +681,8 @@ DEAL:Double:Tan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*s2*
 DEAL:Double:Tan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -(1 + tan(s2/s1)**2)/s1**2 - 2*s2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**3
 DEAL:Double:Tan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -(1 + tan(s2/s1)**2)/s1**2 - 2*s2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**3
 DEAL:Double:Tan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**2
+DEAL:Double:Tan:Symbolic computation: Optimization::Tan: About to optimise... done.
+DEAL:Double:Tan:Symbolic computation: Substitution::Tan: About to substitute... done.
 DEAL:Double:Tan:Symbolic computation: Substitution::subs_psi: 0.7868428895
 DEAL:Double:Tan:Symbolic computation: Substitution::subs_dpsi_ds0: -0.3598048295
 DEAL:Double:Tan:Symbolic computation: Substitution::subs_dpsi_ds1: 0.5397072442
@@ -617,6 +699,8 @@ DEAL:Double:ASin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 1.0*
 DEAL:Double:ASin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ASin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ASin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 1.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 2.0*1/(s1*s2**3*sqrt(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Double:ASin:Symbolic computation: Optimization::ASin: About to optimise... done.
+DEAL:Double:ASin:Symbolic computation: Substitution::ASin: About to substitute... done.
 DEAL:Double:ASin:Symbolic computation: Substitution::subs_psi: 0.1674480792
 DEAL:Double:ASin:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05634361698
 DEAL:Double:ASin:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08451542547
@@ -633,6 +717,8 @@ DEAL:Double:ACos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -1.0
 DEAL:Double:ACos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ACos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ACos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 2.0*1/(s1*s2**3*sqrt(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Double:ACos:Symbolic computation: Optimization::ACos: About to optimise... done.
+DEAL:Double:ACos:Symbolic computation: Substitution::ACos: About to substitute... done.
 DEAL:Double:ACos:Symbolic computation: Substitution::subs_psi: 1.403348248
 DEAL:Double:ACos:Symbolic computation: Substitution::subs_dpsi_ds0: 0.05634361698
 DEAL:Double:ACos:Symbolic computation: Substitution::subs_dpsi_ds1: 0.08451542547
@@ -649,6 +735,8 @@ DEAL:Double:ATan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -2.0
 DEAL:Double:ATan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -2.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ATan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -2.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ATan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -2.0*1/(s1**3*s2**5*(1 + 1.0*1/(s1**2*s2**2))**2) + 2.0*1/(s1*s2**3*(1 + 1.0*1/(s1**2*s2**2)))
+DEAL:Double:ATan:Symbolic computation: Optimization::ATan: About to optimise... done.
+DEAL:Double:ATan:Symbolic computation: Substitution::ATan: About to substitute... done.
 DEAL:Double:ATan:Symbolic computation: Substitution::subs_psi: 0.1651486774
 DEAL:Double:ATan:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05405405405
 DEAL:Double:ATan:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08108108108
@@ -665,6 +753,8 @@ DEAL:Double:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 6.0
 DEAL:Double:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 6.0*s2**2/(s1**2 + s2**2)**2 - 3.0*(s1**2 + s2**2)**(-1)
 DEAL:Double:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -6.0*s1**2/(s1**2 + s2**2)**2 + 3.0*(s1**2 + s2**2)**(-1)
 DEAL:Double:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -6.0*s1*s2/(s1**2 + s2**2)**2
+DEAL:Double:ATan2:Symbolic computation: Optimization::ATan2: About to optimise... done.
+DEAL:Double:ATan2:Symbolic computation: Substitution::ATan2: About to substitute... done.
 DEAL:Double:ATan2:Symbolic computation: Substitution::subs_psi: 1.764007811
 DEAL:Double:ATan2:Symbolic computation: Substitution::subs_dpsi_ds0: -0.4615384615
 DEAL:Double:ATan2:Symbolic computation: Substitution::subs_dpsi_ds1: 0.6923076923
@@ -681,6 +771,8 @@ DEAL:Double:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: s2**
 DEAL:Double:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -cosh(s2/s1)/s1**2 - s2*sinh(s2/s1)/s1**3
 DEAL:Double:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -cosh(s2/s1)/s1**2 - s2*sinh(s2/s1)/s1**3
 DEAL:Double:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: sinh(s2/s1)/s1**2
+DEAL:Double:Sinh:Symbolic computation: Optimization::Sinh: About to optimise... done.
+DEAL:Double:Sinh:Symbolic computation: Substitution::Sinh: About to substitute... done.
 DEAL:Double:Sinh:Symbolic computation: Substitution::subs_psi: 0.7171584610
 DEAL:Double:Sinh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.2734612400
 DEAL:Double:Sinh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.4101918600
@@ -697,6 +789,8 @@ DEAL:Double:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: s2**
 DEAL:Double:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -sinh(s2/s1)/s1**2 - s2*cosh(s2/s1)/s1**3
 DEAL:Double:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -sinh(s2/s1)/s1**2 - s2*cosh(s2/s1)/s1**3
 DEAL:Double:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: cosh(s2/s1)/s1**2
+DEAL:Double:Cosh:Symbolic computation: Optimization::Cosh: About to optimise... done.
+DEAL:Double:Cosh:Symbolic computation: Substitution::Cosh: About to substitute... done.
 DEAL:Double:Cosh:Symbolic computation: Substitution::subs_psi: 1.230575580
 DEAL:Double:Cosh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1593685469
 DEAL:Double:Cosh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2390528203
@@ -713,6 +807,8 @@ DEAL:Double:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*s2
 DEAL:Double:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -(1 - tanh(s2/s1)**2)/s1**2 + 2*s2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**3
 DEAL:Double:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -(1 - tanh(s2/s1)**2)/s1**2 + 2*s2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**3
 DEAL:Double:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**2
+DEAL:Double:Tanh:Symbolic computation: Optimization::Tanh: About to optimise... done.
+DEAL:Double:Tanh:Symbolic computation: Substitution::Tanh: About to substitute... done.
 DEAL:Double:Tanh:Symbolic computation: Substitution::subs_psi: 0.5827829453
 DEAL:Double:Tanh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1467475641
 DEAL:Double:Tanh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2201213462
@@ -729,6 +825,8 @@ DEAL:Double:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -1.
 DEAL:Double:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1.0*1/(s1**3*s2**5*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 2.0*1/(s1*s2**3*sqrt(1 + 1.0*1/(s1**2*s2**2)))
+DEAL:Double:ASinh:Symbolic computation: Optimization::ASinh: About to optimise... done.
+DEAL:Double:ASinh:Symbolic computation: Substitution::ASinh: About to substitute... done.
 DEAL:Double:ASinh:Symbolic computation: Substitution::subs_psi: 0.1659045503
 DEAL:Double:ASinh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05479966244
 DEAL:Double:ASinh:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08219949365
@@ -745,6 +843,8 @@ DEAL:Double:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s1
 DEAL:Double:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -s1**2*s2**2/(-1 + s1**2*s2**2)**(3/2) + (-1 + s1**2*s2**2)**(-1/2)
 DEAL:Double:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -s1**2*s2**2/(-1 + s1**2*s2**2)**(3/2) + (-1 + s1**2*s2**2)**(-1/2)
 DEAL:Double:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -s1**3*s2/(-1 + s1**2*s2**2)**(3/2)
+DEAL:Double:ACosh:Symbolic computation: Optimization::ACosh: About to optimise... done.
+DEAL:Double:ACosh:Symbolic computation: Substitution::ACosh: About to substitute... done.
 DEAL:Double:ACosh:Symbolic computation: Substitution::subs_psi: 2.477888730
 DEAL:Double:ACosh:Symbolic computation: Substitution::subs_dpsi_ds0: 0.3380617019
 DEAL:Double:ACosh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.5070925528
@@ -761,6 +861,8 @@ DEAL:Double:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2.0
 DEAL:Double:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**2) + 2.0*1/(s1*s2**3*(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Double:ATanh:Symbolic computation: Optimization::ATanh: About to optimise... done.
+DEAL:Double:ATanh:Symbolic computation: Substitution::ATanh: About to substitute... done.
 DEAL:Double:ATanh:Symbolic computation: Substitution::subs_psi: 0.1682361183
 DEAL:Double:ATanh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05714285714
 DEAL:Double:ATanh:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08571428571
@@ -777,6 +879,8 @@ DEAL:Double:Erf:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*erf
 DEAL:Double:Erf:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Double:Erf:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Double:Erf:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -4*exp(-s2**2)*s1**2*s2/sqrt(pi)
+DEAL:Double:Erf:Symbolic computation: Optimization::Erf: About to optimise... done.
+DEAL:Double:Erf:Symbolic computation: Substitution::Erf: About to substitute... done.
 DEAL:Double:Erf:Symbolic computation: Substitution::subs_psi: 8.957900385
 DEAL:Double:Erf:Symbolic computation: Substitution::subs_dpsi_ds0: 5.971933590
 DEAL:Double:Erf:Symbolic computation: Substitution::subs_dpsi_ds1: 0.1860028682
@@ -793,6 +897,8 @@ DEAL:Double:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*er
 DEAL:Double:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Double:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Double:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 4*exp(-s2**2)*s1**2*s2/sqrt(pi)
+DEAL:Double:Erfc:Symbolic computation: Optimization::Erfc: About to optimise... done.
+DEAL:Double:Erfc:Symbolic computation: Substitution::Erfc: About to substitute... done.
 DEAL:Double:Erfc:Symbolic computation: Substitution::subs_psi: 0.04209961483
 DEAL:Double:Erfc:Symbolic computation: Substitution::subs_dpsi_ds0: 0.02806640989
 DEAL:Double:Erfc:Symbolic computation: Substitution::subs_dpsi_ds1: -0.1860028682
index 936f3726c4d1ba4d660c03d10c8019d8e23d4938..1e71050a8b87d0a62011aa8085276973a98d7748 100644 (file)
@@ -1,4 +1,7 @@
 
+DEAL:Float::Optimisation method: llvm
+Optimisation flags:  OptimizationFlags|cse|-O3|
+Skip non-implemented functions: 0
 DEAL:Float:Add:Symbolic computation: Function initialisation::symb_s0: s1
 DEAL:Float:Add:Symbolic computation: Function initialisation::symb_s1: s2
 DEAL:Float:Add:Symbolic computation: Function initialisation::symb_psi: (s1 + s2)**2
@@ -8,6 +11,8 @@ DEAL:Float:Add:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2
 DEAL:Float:Add:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2
 DEAL:Float:Add:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2
 DEAL:Float:Add:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2
+DEAL:Float:Add:Symbolic computation: Optimization::Add: About to optimise... done.
+DEAL:Float:Add:Symbolic computation: Substitution::Add: About to substitute... done.
 DEAL:Float:Add:Symbolic computation: Substitution::subs_psi: 25.00000000
 DEAL:Float:Add:Symbolic computation: Substitution::subs_dpsi_ds0: 10.00000000
 DEAL:Float:Add:Symbolic computation: Substitution::subs_dpsi_ds1: 10.00000000
@@ -24,6 +29,8 @@ DEAL:Float:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2
 DEAL:Float:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -2
 DEAL:Float:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -2
 DEAL:Float:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2
+DEAL:Float:Subtract:Symbolic computation: Optimization::Subtract: About to optimise... done.
+DEAL:Float:Subtract:Symbolic computation: Substitution::Subtract: About to substitute... done.
 DEAL:Float:Subtract:Symbolic computation: Substitution::subs_psi: 1.000000000
 DEAL:Float:Subtract:Symbolic computation: Substitution::subs_dpsi_ds0: 2.000000000
 DEAL:Float:Subtract:Symbolic computation: Substitution::subs_dpsi_ds1: -2.000000000
@@ -40,6 +47,8 @@ DEAL:Float:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2
 DEAL:Float:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 4*s1*s2
 DEAL:Float:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 4*s1*s2
 DEAL:Float:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2*s1**2
+DEAL:Float:Multiply:Symbolic computation: Optimization::Multiply: About to optimise... done.
+DEAL:Float:Multiply:Symbolic computation: Substitution::Multiply: About to substitute... done.
 DEAL:Float:Multiply:Symbolic computation: Substitution::subs_psi: 36.00000000
 DEAL:Float:Multiply:Symbolic computation: Substitution::subs_dpsi_ds0: 24.00000000
 DEAL:Float:Multiply:Symbolic computation: Substitution::subs_dpsi_ds1: 36.00000000
@@ -56,6 +65,8 @@ DEAL:Float:Divide:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2/s
 DEAL:Float:Divide:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -4*s1/s2**3
 DEAL:Float:Divide:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -4*s1/s2**3
 DEAL:Float:Divide:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 6*s1**2/s2**4
+DEAL:Float:Divide:Symbolic computation: Optimization::Divide: About to optimise... done.
+DEAL:Float:Divide:Symbolic computation: Substitution::Divide: About to substitute... done.
 DEAL:Float:Divide:Symbolic computation: Substitution::subs_psi: 2.250000000
 DEAL:Float:Divide:Symbolic computation: Substitution::subs_dpsi_ds0: 1.500000000
 DEAL:Float:Divide:Symbolic computation: Substitution::subs_dpsi_ds1: -2.250000000
@@ -72,6 +83,8 @@ DEAL:Float:Power:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2.0*
 DEAL:Float:Power:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2.0*s1**(-1 + s2)*s2*log(s1) + 2.0*s1**(-1 + s2)
 DEAL:Float:Power:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2.0*s1**(-1 + s2)*s2*log(s1) + 2.0*s1**(-1 + s2)
 DEAL:Float:Power:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2.0*s1**s2*log(s1)**2
+DEAL:Float:Power:Symbolic computation: Optimization::Power: About to optimise... done.
+DEAL:Float:Power:Symbolic computation: Substitution::Power: About to substitute... done.
 DEAL:Float:Power:Symbolic computation: Substitution::subs_psi: 18.00000000
 DEAL:Float:Power:Symbolic computation: Substitution::subs_dpsi_ds0: 12.00000000
 DEAL:Float:Power:Symbolic computation: Substitution::subs_dpsi_ds1: 19.77502060
@@ -88,6 +101,8 @@ DEAL:Float:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -0.5*
 DEAL:Float:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -0.5*s1*s2/(s1*s2)**(3/2) + 1.0*(s1*s2)**(-1/2)
 DEAL:Float:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -0.5*s1*s2/(s1*s2)**(3/2) + 1.0*(s1*s2)**(-1/2)
 DEAL:Float:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -0.5*s1**2/(s1*s2)**(3/2)
+DEAL:Float:Sqrt:Symbolic computation: Optimization::Sqrt: About to optimise... done.
+DEAL:Float:Sqrt:Symbolic computation: Substitution::Sqrt: About to substitute... done.
 DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_psi: 4.898979664
 DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_dpsi_ds0: 0.8164965510
 DEAL:Float:Sqrt:Symbolic computation: Substitution::subs_dpsi_ds1: 1.224744797
@@ -104,6 +119,8 @@ DEAL:Float:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0
 DEAL:Float:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2*exp(s2)*s1
 DEAL:Float:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2*exp(s2)*s1
 DEAL:Float:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: exp(s2)*s1**2
+DEAL:Float:Exponential:Symbolic computation: Optimization::Exponential: About to optimise... done.
+DEAL:Float:Exponential:Symbolic computation: Substitution::Exponential: About to substitute... done.
 DEAL:Float:Exponential:Symbolic computation: Substitution::subs_psi: 66.50150299
 DEAL:Float:Exponential:Symbolic computation: Substitution::subs_dpsi_ds0: 44.33433533
 DEAL:Float:Exponential:Symbolic computation: Substitution::subs_dpsi_ds1: 66.50150299
@@ -120,6 +137,8 @@ DEAL:Float:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0:
 DEAL:Float:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2*s1/s2
 DEAL:Float:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2*s1/s2
 DEAL:Float:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -s1**2/s2**2
+DEAL:Float:Logarithm:Symbolic computation: Optimization::Logarithm: About to optimise... done.
+DEAL:Float:Logarithm:Symbolic computation: Substitution::Logarithm: About to substitute... done.
 DEAL:Float:Logarithm:Symbolic computation: Substitution::subs_psi: 6.238324642
 DEAL:Float:Logarithm:Symbolic computation: Substitution::subs_dpsi_ds0: 4.158883095
 DEAL:Float:Logarithm:Symbolic computation: Substitution::subs_dpsi_ds1: 4.500000000
@@ -136,6 +155,8 @@ DEAL:Float:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi
 DEAL:Float:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1/(s1*s2*log(s1)**2)
 DEAL:Float:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1/(s1*s2*log(s1)**2)
 DEAL:Float:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1/(s2**2*log(s1))
+DEAL:Float:Logarithm with Base:Symbolic computation: Optimization::Logarithm with Base: About to optimise... done.
+DEAL:Float:Logarithm with Base:Symbolic computation: Substitution::Logarithm with Base: About to substitute... done.
 DEAL:Float:Logarithm with Base:Symbolic computation: Substitution::subs_psi: 0.6309297681
 DEAL:Float:Logarithm with Base:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1914323419
 DEAL:Float:Logarithm with Base:Symbolic computation: Substitution::subs_dpsi_ds1: 0.4551196098
@@ -152,6 +173,8 @@ DEAL:Float:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_d
 DEAL:Float:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 0.868588963806504*s1/s2
 DEAL:Float:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 0.868588963806504*s1/s2
 DEAL:Float:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -0.434294481903252*s1**2/s2**2
+DEAL:Float:Logarithm base 10:Symbolic computation: Optimization::Logarithm base 10: About to optimise... done.
+DEAL:Float:Logarithm base 10:Symbolic computation: Substitution::Logarithm base 10: About to substitute... done.
 DEAL:Float:Logarithm base 10:Symbolic computation: Substitution::subs_psi: 2.709270000
 DEAL:Float:Logarithm base 10:Symbolic computation: Substitution::subs_dpsi_ds0: 1.806180000
 DEAL:Float:Logarithm base 10:Symbolic computation: Substitution::subs_dpsi_ds1: 1.954325199
@@ -168,6 +191,8 @@ DEAL:Float:Sin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s2**2
 DEAL:Float:Sin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -cos(s2/s1)/s1**2 + s2*sin(s2/s1)/s1**3
 DEAL:Float:Sin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -cos(s2/s1)/s1**2 + s2*sin(s2/s1)/s1**3
 DEAL:Float:Sin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -sin(s2/s1)/s1**2
+DEAL:Float:Sin:Symbolic computation: Optimization::Sin: About to optimise... done.
+DEAL:Float:Sin:Symbolic computation: Substitution::Sin: About to substitute... done.
 DEAL:Float:Sin:Symbolic computation: Substitution::subs_psi: 0.6183698177
 DEAL:Float:Sin:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1746416092
 DEAL:Float:Sin:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2619624138
@@ -184,6 +209,8 @@ DEAL:Float:Cos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s2**2
 DEAL:Float:Cos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: sin(s2/s1)/s1**2 + s2*cos(s2/s1)/s1**3
 DEAL:Float:Cos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: sin(s2/s1)/s1**2 + s2*cos(s2/s1)/s1**3
 DEAL:Float:Cos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -cos(s2/s1)/s1**2
+DEAL:Float:Cos:Symbolic computation: Optimization::Cos: About to optimise... done.
+DEAL:Float:Cos:Symbolic computation: Substitution::Cos: About to substitute... done.
 DEAL:Float:Cos:Symbolic computation: Substitution::subs_psi: 0.7858872414
 DEAL:Float:Cos:Symbolic computation: Substitution::subs_dpsi_ds0: 0.1374155134
 DEAL:Float:Cos:Symbolic computation: Substitution::subs_dpsi_ds1: -0.2061232775
@@ -200,6 +227,8 @@ DEAL:Float:Tan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*s2*(
 DEAL:Float:Tan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -(1 + tan(s2/s1)**2)/s1**2 - 2*s2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**3
 DEAL:Float:Tan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -(1 + tan(s2/s1)**2)/s1**2 - 2*s2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**3
 DEAL:Float:Tan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**2
+DEAL:Float:Tan:Symbolic computation: Optimization::Tan: About to optimise... done.
+DEAL:Float:Tan:Symbolic computation: Substitution::Tan: About to substitute... done.
 DEAL:Float:Tan:Symbolic computation: Substitution::subs_psi: 0.7868429422
 DEAL:Float:Tan:Symbolic computation: Substitution::subs_dpsi_ds0: -0.3598048389
 DEAL:Float:Tan:Symbolic computation: Substitution::subs_dpsi_ds1: 0.5397073030
@@ -216,6 +245,8 @@ DEAL:Float:ASin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 1.0*1
 DEAL:Float:ASin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ASin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ASin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 1.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 2.0*1/(s1*s2**3*sqrt(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Float:ASin:Symbolic computation: Optimization::ASin: About to optimise... done.
+DEAL:Float:ASin:Symbolic computation: Substitution::ASin: About to substitute... done.
 DEAL:Float:ASin:Symbolic computation: Substitution::subs_psi: 0.1674480885
 DEAL:Float:ASin:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05634361506
 DEAL:Float:ASin:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08451542258
@@ -232,6 +263,8 @@ DEAL:Float:ACos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -1.0*
 DEAL:Float:ACos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ACos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ACos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 2.0*1/(s1*s2**3*sqrt(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Float:ACos:Symbolic computation: Optimization::ACos: About to optimise... done.
+DEAL:Float:ACos:Symbolic computation: Substitution::ACos: About to substitute... done.
 DEAL:Float:ACos:Symbolic computation: Substitution::subs_psi: 1.403348207
 DEAL:Float:ACos:Symbolic computation: Substitution::subs_dpsi_ds0: 0.05634361506
 DEAL:Float:ACos:Symbolic computation: Substitution::subs_dpsi_ds1: 0.08451542258
@@ -248,7 +281,9 @@ DEAL:Float:ATan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -2.0*
 DEAL:Float:ATan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -2.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ATan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -2.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ATan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -2.0*1/(s1**3*s2**5*(1 + 1.0*1/(s1**2*s2**2))**2) + 2.0*1/(s1*s2**3*(1 + 1.0*1/(s1**2*s2**2)))
-DEAL:Float:ATan:Symbolic computation: Substitution::subs_psi: 0.1651486754
+DEAL:Float:ATan:Symbolic computation: Optimization::ATan: About to optimise... done.
+DEAL:Float:ATan:Symbolic computation: Substitution::ATan: About to substitute... done.
+DEAL:Float:ATan:Symbolic computation: Substitution::subs_psi: 0.1651486903
 DEAL:Float:ATan:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05405405536
 DEAL:Float:ATan:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08108108491
 DEAL:Float:ATan:Symbolic computation: Substitution::subs_d2psi_ds0_ds0: 0.03506208956
@@ -264,6 +299,8 @@ DEAL:Float:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 6.0*
 DEAL:Float:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 6.0*s2**2/(s1**2 + s2**2)**2 - 3.0*(s1**2 + s2**2)**(-1)
 DEAL:Float:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -6.0*s1**2/(s1**2 + s2**2)**2 + 3.0*(s1**2 + s2**2)**(-1)
 DEAL:Float:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -6.0*s1*s2/(s1**2 + s2**2)**2
+DEAL:Float:ATan2:Symbolic computation: Optimization::ATan2: About to optimise... done.
+DEAL:Float:ATan2:Symbolic computation: Substitution::ATan2: About to substitute... done.
 DEAL:Float:ATan2:Symbolic computation: Substitution::subs_psi: 1.764007807
 DEAL:Float:ATan2:Symbolic computation: Substitution::subs_dpsi_ds0: -0.4615384936
 DEAL:Float:ATan2:Symbolic computation: Substitution::subs_dpsi_ds1: 0.6923077106
@@ -280,6 +317,8 @@ DEAL:Float:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: s2**2
 DEAL:Float:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -cosh(s2/s1)/s1**2 - s2*sinh(s2/s1)/s1**3
 DEAL:Float:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -cosh(s2/s1)/s1**2 - s2*sinh(s2/s1)/s1**3
 DEAL:Float:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: sinh(s2/s1)/s1**2
+DEAL:Float:Sinh:Symbolic computation: Optimization::Sinh: About to optimise... done.
+DEAL:Float:Sinh:Symbolic computation: Substitution::Sinh: About to substitute... done.
 DEAL:Float:Sinh:Symbolic computation: Substitution::subs_psi: 0.7171584964
 DEAL:Float:Sinh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.2734612525
 DEAL:Float:Sinh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.4101918638
@@ -296,6 +335,8 @@ DEAL:Float:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: s2**2
 DEAL:Float:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -sinh(s2/s1)/s1**2 - s2*cosh(s2/s1)/s1**3
 DEAL:Float:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -sinh(s2/s1)/s1**2 - s2*cosh(s2/s1)/s1**3
 DEAL:Float:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: cosh(s2/s1)/s1**2
+DEAL:Float:Cosh:Symbolic computation: Optimization::Cosh: About to optimise... done.
+DEAL:Float:Cosh:Symbolic computation: Substitution::Cosh: About to substitute... done.
 DEAL:Float:Cosh:Symbolic computation: Substitution::subs_psi: 1.230575562
 DEAL:Float:Cosh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1593685597
 DEAL:Float:Cosh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2390528321
@@ -312,6 +353,8 @@ DEAL:Float:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*s2*
 DEAL:Float:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -(1 - tanh(s2/s1)**2)/s1**2 + 2*s2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**3
 DEAL:Float:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -(1 - tanh(s2/s1)**2)/s1**2 + 2*s2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**3
 DEAL:Float:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**2
+DEAL:Float:Tanh:Symbolic computation: Optimization::Tanh: About to optimise... done.
+DEAL:Float:Tanh:Symbolic computation: Substitution::Tanh: About to substitute... done.
 DEAL:Float:Tanh:Symbolic computation: Substitution::subs_psi: 0.5827829838
 DEAL:Float:Tanh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1467475593
 DEAL:Float:Tanh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2201213539
@@ -328,6 +371,8 @@ DEAL:Float:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -1.0
 DEAL:Float:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1.0*1/(s1**3*s2**5*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 2.0*1/(s1*s2**3*sqrt(1 + 1.0*1/(s1**2*s2**2)))
+DEAL:Float:ASinh:Symbolic computation: Optimization::ASinh: About to optimise... done.
+DEAL:Float:ASinh:Symbolic computation: Substitution::ASinh: About to substitute... done.
 DEAL:Float:ASinh:Symbolic computation: Substitution::subs_psi: 0.1659045517
 DEAL:Float:ASinh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05479966477
 DEAL:Float:ASinh:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08219949901
@@ -344,6 +389,8 @@ DEAL:Float:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s1*
 DEAL:Float:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -s1**2*s2**2/(-1 + s1**2*s2**2)**(3/2) + (-1 + s1**2*s2**2)**(-1/2)
 DEAL:Float:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -s1**2*s2**2/(-1 + s1**2*s2**2)**(3/2) + (-1 + s1**2*s2**2)**(-1/2)
 DEAL:Float:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -s1**3*s2/(-1 + s1**2*s2**2)**(3/2)
+DEAL:Float:ACosh:Symbolic computation: Optimization::ACosh: About to optimise... done.
+DEAL:Float:ACosh:Symbolic computation: Substitution::ACosh: About to substitute... done.
 DEAL:Float:ACosh:Symbolic computation: Substitution::subs_psi: 2.477888823
 DEAL:Float:ACosh:Symbolic computation: Substitution::subs_dpsi_ds0: 0.3380616903
 DEAL:Float:ACosh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.5070925355
@@ -360,6 +407,8 @@ DEAL:Float:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2.0*
 DEAL:Float:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Float:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**2) + 2.0*1/(s1*s2**3*(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Float:ATanh:Symbolic computation: Optimization::ATanh: About to optimise... done.
+DEAL:Float:ATanh:Symbolic computation: Substitution::ATanh: About to substitute... done.
 DEAL:Float:ATanh:Symbolic computation: Substitution::subs_psi: 0.1682361215
 DEAL:Float:ATanh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05714286119
 DEAL:Float:ATanh:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08571429551
@@ -376,6 +425,8 @@ DEAL:Float:Erf:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*erf(
 DEAL:Float:Erf:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Float:Erf:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Float:Erf:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -4*exp(-s2**2)*s1**2*s2/sqrt(pi)
+DEAL:Float:Erf:Symbolic computation: Optimization::Erf: About to optimise... done.
+DEAL:Float:Erf:Symbolic computation: Substitution::Erf: About to substitute... done.
 DEAL:Float:Erf:Symbolic computation: Substitution::subs_psi: 8.957901001
 DEAL:Float:Erf:Symbolic computation: Substitution::subs_dpsi_ds0: 5.971933842
 DEAL:Float:Erf:Symbolic computation: Substitution::subs_dpsi_ds1: 0.1860028654
@@ -392,6 +443,8 @@ DEAL:Float:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*erf
 DEAL:Float:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Float:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Float:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 4*exp(-s2**2)*s1**2*s2/sqrt(pi)
+DEAL:Float:Erfc:Symbolic computation: Optimization::Erfc: About to optimise... done.
+DEAL:Float:Erfc:Symbolic computation: Substitution::Erfc: About to substitute... done.
 DEAL:Float:Erfc:Symbolic computation: Substitution::subs_psi: 0.04209961370
 DEAL:Float:Erfc:Symbolic computation: Substitution::subs_dpsi_ds0: 0.02806640789
 DEAL:Float:Erfc:Symbolic computation: Substitution::subs_dpsi_ds1: -0.1860028654
@@ -400,6 +453,9 @@ DEAL:Float:Erfc:Symbolic computation: Substitution::subs_d2psi_ds1_ds0: -0.12400
 DEAL:Float:Erfc:Symbolic computation: Substitution::subs_d2psi_ds0_ds1: -0.1240019053
 DEAL:Float:Erfc:Symbolic computation: Substitution::subs_d2psi_ds1_ds1: 0.7440114617
 DEAL:Float::OK
+DEAL:Double::Optimisation method: llvm
+Optimisation flags:  OptimizationFlags|cse|-O3|
+Skip non-implemented functions: 0
 DEAL:Double:Add:Symbolic computation: Function initialisation::symb_s0: s1
 DEAL:Double:Add:Symbolic computation: Function initialisation::symb_s1: s2
 DEAL:Double:Add:Symbolic computation: Function initialisation::symb_psi: (s1 + s2)**2
@@ -409,6 +465,8 @@ DEAL:Double:Add:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2
 DEAL:Double:Add:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2
 DEAL:Double:Add:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2
 DEAL:Double:Add:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2
+DEAL:Double:Add:Symbolic computation: Optimization::Add: About to optimise... done.
+DEAL:Double:Add:Symbolic computation: Substitution::Add: About to substitute... done.
 DEAL:Double:Add:Symbolic computation: Substitution::subs_psi: 25.00000000
 DEAL:Double:Add:Symbolic computation: Substitution::subs_dpsi_ds0: 10.00000000
 DEAL:Double:Add:Symbolic computation: Substitution::subs_dpsi_ds1: 10.00000000
@@ -425,6 +483,8 @@ DEAL:Double:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0:
 DEAL:Double:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -2
 DEAL:Double:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -2
 DEAL:Double:Subtract:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2
+DEAL:Double:Subtract:Symbolic computation: Optimization::Subtract: About to optimise... done.
+DEAL:Double:Subtract:Symbolic computation: Substitution::Subtract: About to substitute... done.
 DEAL:Double:Subtract:Symbolic computation: Substitution::subs_psi: 1.000000000
 DEAL:Double:Subtract:Symbolic computation: Substitution::subs_dpsi_ds0: 2.000000000
 DEAL:Double:Subtract:Symbolic computation: Substitution::subs_dpsi_ds1: -2.000000000
@@ -441,6 +501,8 @@ DEAL:Double:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0:
 DEAL:Double:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 4*s1*s2
 DEAL:Double:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 4*s1*s2
 DEAL:Double:Multiply:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2*s1**2
+DEAL:Double:Multiply:Symbolic computation: Optimization::Multiply: About to optimise... done.
+DEAL:Double:Multiply:Symbolic computation: Substitution::Multiply: About to substitute... done.
 DEAL:Double:Multiply:Symbolic computation: Substitution::subs_psi: 36.00000000
 DEAL:Double:Multiply:Symbolic computation: Substitution::subs_dpsi_ds0: 24.00000000
 DEAL:Double:Multiply:Symbolic computation: Substitution::subs_dpsi_ds1: 36.00000000
@@ -457,6 +519,8 @@ DEAL:Double:Divide:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2/
 DEAL:Double:Divide:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -4*s1/s2**3
 DEAL:Double:Divide:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -4*s1/s2**3
 DEAL:Double:Divide:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 6*s1**2/s2**4
+DEAL:Double:Divide:Symbolic computation: Optimization::Divide: About to optimise... done.
+DEAL:Double:Divide:Symbolic computation: Substitution::Divide: About to substitute... done.
 DEAL:Double:Divide:Symbolic computation: Substitution::subs_psi: 2.250000000
 DEAL:Double:Divide:Symbolic computation: Substitution::subs_dpsi_ds0: 1.500000000
 DEAL:Double:Divide:Symbolic computation: Substitution::subs_dpsi_ds1: -2.250000000
@@ -473,6 +537,8 @@ DEAL:Double:Power:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2.0
 DEAL:Double:Power:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2.0*s1**(-1 + s2)*s2*log(s1) + 2.0*s1**(-1 + s2)
 DEAL:Double:Power:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2.0*s1**(-1 + s2)*s2*log(s1) + 2.0*s1**(-1 + s2)
 DEAL:Double:Power:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2.0*s1**s2*log(s1)**2
+DEAL:Double:Power:Symbolic computation: Optimization::Power: About to optimise... done.
+DEAL:Double:Power:Symbolic computation: Substitution::Power: About to substitute... done.
 DEAL:Double:Power:Symbolic computation: Substitution::subs_psi: 18.00000000
 DEAL:Double:Power:Symbolic computation: Substitution::subs_dpsi_ds0: 12.00000000
 DEAL:Double:Power:Symbolic computation: Substitution::subs_dpsi_ds1: 19.77502120
@@ -489,6 +555,8 @@ DEAL:Double:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -0.5
 DEAL:Double:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -0.5*s1*s2/(s1*s2)**(3/2) + 1.0*(s1*s2)**(-1/2)
 DEAL:Double:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -0.5*s1*s2/(s1*s2)**(3/2) + 1.0*(s1*s2)**(-1/2)
 DEAL:Double:Sqrt:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -0.5*s1**2/(s1*s2)**(3/2)
+DEAL:Double:Sqrt:Symbolic computation: Optimization::Sqrt: About to optimise... done.
+DEAL:Double:Sqrt:Symbolic computation: Substitution::Sqrt: About to substitute... done.
 DEAL:Double:Sqrt:Symbolic computation: Substitution::subs_psi: 4.898979486
 DEAL:Double:Sqrt:Symbolic computation: Substitution::subs_dpsi_ds0: 0.8164965809
 DEAL:Double:Sqrt:Symbolic computation: Substitution::subs_dpsi_ds1: 1.224744871
@@ -505,6 +573,8 @@ DEAL:Double:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds0_ds
 DEAL:Double:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2*exp(s2)*s1
 DEAL:Double:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2*exp(s2)*s1
 DEAL:Double:Exponential:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: exp(s2)*s1**2
+DEAL:Double:Exponential:Symbolic computation: Optimization::Exponential: About to optimise... done.
+DEAL:Double:Exponential:Symbolic computation: Substitution::Exponential: About to substitute... done.
 DEAL:Double:Exponential:Symbolic computation: Substitution::subs_psi: 66.50150489
 DEAL:Double:Exponential:Symbolic computation: Substitution::subs_dpsi_ds0: 44.33433659
 DEAL:Double:Exponential:Symbolic computation: Substitution::subs_dpsi_ds1: 66.50150489
@@ -521,6 +591,8 @@ DEAL:Double:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0:
 DEAL:Double:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2*s1/s2
 DEAL:Double:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2*s1/s2
 DEAL:Double:Logarithm:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -s1**2/s2**2
+DEAL:Double:Logarithm:Symbolic computation: Optimization::Logarithm: About to optimise... done.
+DEAL:Double:Logarithm:Symbolic computation: Substitution::Logarithm: About to substitute... done.
 DEAL:Double:Logarithm:Symbolic computation: Substitution::subs_psi: 6.238324625
 DEAL:Double:Logarithm:Symbolic computation: Substitution::subs_dpsi_ds0: 4.158883083
 DEAL:Double:Logarithm:Symbolic computation: Substitution::subs_dpsi_ds1: 4.500000000
@@ -537,6 +609,8 @@ DEAL:Double:Logarithm with Base:Symbolic computation: Differentiation::symb_d2ps
 DEAL:Double:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1/(s1*s2*log(s1)**2)
 DEAL:Double:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1/(s1*s2*log(s1)**2)
 DEAL:Double:Logarithm with Base:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1/(s2**2*log(s1))
+DEAL:Double:Logarithm with Base:Symbolic computation: Optimization::Logarithm with Base: About to optimise... done.
+DEAL:Double:Logarithm with Base:Symbolic computation: Substitution::Logarithm with Base: About to substitute... done.
 DEAL:Double:Logarithm with Base:Symbolic computation: Substitution::subs_psi: 0.6309297536
 DEAL:Double:Logarithm with Base:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1914323370
 DEAL:Double:Logarithm with Base:Symbolic computation: Substitution::subs_dpsi_ds1: 0.4551196133
@@ -553,6 +627,8 @@ DEAL:Double:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_
 DEAL:Double:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 0.868588963806504*s1/s2
 DEAL:Double:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 0.868588963806504*s1/s2
 DEAL:Double:Logarithm base 10:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -0.434294481903252*s1**2/s2**2
+DEAL:Double:Logarithm base 10:Symbolic computation: Optimization::Logarithm base 10: About to optimise... done.
+DEAL:Double:Logarithm base 10:Symbolic computation: Substitution::Logarithm base 10: About to substitute... done.
 DEAL:Double:Logarithm base 10:Symbolic computation: Substitution::subs_psi: 2.709269961
 DEAL:Double:Logarithm base 10:Symbolic computation: Substitution::subs_dpsi_ds0: 1.806179974
 DEAL:Double:Logarithm base 10:Symbolic computation: Substitution::subs_dpsi_ds1: 1.954325169
@@ -569,6 +645,8 @@ DEAL:Double:Sin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s2**
 DEAL:Double:Sin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -cos(s2/s1)/s1**2 + s2*sin(s2/s1)/s1**3
 DEAL:Double:Sin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -cos(s2/s1)/s1**2 + s2*sin(s2/s1)/s1**3
 DEAL:Double:Sin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -sin(s2/s1)/s1**2
+DEAL:Double:Sin:Symbolic computation: Optimization::Sin: About to optimise... done.
+DEAL:Double:Sin:Symbolic computation: Substitution::Sin: About to substitute... done.
 DEAL:Double:Sin:Symbolic computation: Substitution::subs_psi: 0.6183698031
 DEAL:Double:Sin:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1746416135
 DEAL:Double:Sin:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2619624203
@@ -585,6 +663,8 @@ DEAL:Double:Cos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s2**
 DEAL:Double:Cos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: sin(s2/s1)/s1**2 + s2*cos(s2/s1)/s1**3
 DEAL:Double:Cos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: sin(s2/s1)/s1**2 + s2*cos(s2/s1)/s1**3
 DEAL:Double:Cos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -cos(s2/s1)/s1**2
+DEAL:Double:Cos:Symbolic computation: Optimization::Cos: About to optimise... done.
+DEAL:Double:Cos:Symbolic computation: Substitution::Cos: About to substitute... done.
 DEAL:Double:Cos:Symbolic computation: Substitution::subs_psi: 0.7858872608
 DEAL:Double:Cos:Symbolic computation: Substitution::subs_dpsi_ds0: 0.1374155118
 DEAL:Double:Cos:Symbolic computation: Substitution::subs_dpsi_ds1: -0.2061232677
@@ -601,6 +681,8 @@ DEAL:Double:Tan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*s2*
 DEAL:Double:Tan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -(1 + tan(s2/s1)**2)/s1**2 - 2*s2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**3
 DEAL:Double:Tan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -(1 + tan(s2/s1)**2)/s1**2 - 2*s2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**3
 DEAL:Double:Tan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2*tan(s2/s1)*(1 + tan(s2/s1)**2)/s1**2
+DEAL:Double:Tan:Symbolic computation: Optimization::Tan: About to optimise... done.
+DEAL:Double:Tan:Symbolic computation: Substitution::Tan: About to substitute... done.
 DEAL:Double:Tan:Symbolic computation: Substitution::subs_psi: 0.7868428895
 DEAL:Double:Tan:Symbolic computation: Substitution::subs_dpsi_ds0: -0.3598048295
 DEAL:Double:Tan:Symbolic computation: Substitution::subs_dpsi_ds1: 0.5397072442
@@ -617,6 +699,8 @@ DEAL:Double:ASin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 1.0*
 DEAL:Double:ASin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ASin:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ASin:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 1.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) + 2.0*1/(s1*s2**3*sqrt(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Double:ASin:Symbolic computation: Optimization::ASin: About to optimise... done.
+DEAL:Double:ASin:Symbolic computation: Substitution::ASin: About to substitute... done.
 DEAL:Double:ASin:Symbolic computation: Substitution::subs_psi: 0.1674480792
 DEAL:Double:ASin:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05634361698
 DEAL:Double:ASin:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08451542547
@@ -633,6 +717,8 @@ DEAL:Double:ACos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -1.0
 DEAL:Double:ACos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ACos:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 1.0*1/(s1**2*s2**2*sqrt(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ACos:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**(3/2)) - 2.0*1/(s1*s2**3*sqrt(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Double:ACos:Symbolic computation: Optimization::ACos: About to optimise... done.
+DEAL:Double:ACos:Symbolic computation: Substitution::ACos: About to substitute... done.
 DEAL:Double:ACos:Symbolic computation: Substitution::subs_psi: 1.403348248
 DEAL:Double:ACos:Symbolic computation: Substitution::subs_dpsi_ds0: 0.05634361698
 DEAL:Double:ACos:Symbolic computation: Substitution::subs_dpsi_ds1: 0.08451542547
@@ -649,6 +735,8 @@ DEAL:Double:ATan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -2.0
 DEAL:Double:ATan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -2.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ATan:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -2.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ATan:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -2.0*1/(s1**3*s2**5*(1 + 1.0*1/(s1**2*s2**2))**2) + 2.0*1/(s1*s2**3*(1 + 1.0*1/(s1**2*s2**2)))
+DEAL:Double:ATan:Symbolic computation: Optimization::ATan: About to optimise... done.
+DEAL:Double:ATan:Symbolic computation: Substitution::ATan: About to substitute... done.
 DEAL:Double:ATan:Symbolic computation: Substitution::subs_psi: 0.1651486774
 DEAL:Double:ATan:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05405405405
 DEAL:Double:ATan:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08108108108
@@ -665,6 +753,8 @@ DEAL:Double:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 6.0
 DEAL:Double:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 6.0*s2**2/(s1**2 + s2**2)**2 - 3.0*(s1**2 + s2**2)**(-1)
 DEAL:Double:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -6.0*s1**2/(s1**2 + s2**2)**2 + 3.0*(s1**2 + s2**2)**(-1)
 DEAL:Double:ATan2:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -6.0*s1*s2/(s1**2 + s2**2)**2
+DEAL:Double:ATan2:Symbolic computation: Optimization::ATan2: About to optimise... done.
+DEAL:Double:ATan2:Symbolic computation: Substitution::ATan2: About to substitute... done.
 DEAL:Double:ATan2:Symbolic computation: Substitution::subs_psi: 1.764007811
 DEAL:Double:ATan2:Symbolic computation: Substitution::subs_dpsi_ds0: -0.4615384615
 DEAL:Double:ATan2:Symbolic computation: Substitution::subs_dpsi_ds1: 0.6923076923
@@ -681,6 +771,8 @@ DEAL:Double:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: s2**
 DEAL:Double:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -cosh(s2/s1)/s1**2 - s2*sinh(s2/s1)/s1**3
 DEAL:Double:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -cosh(s2/s1)/s1**2 - s2*sinh(s2/s1)/s1**3
 DEAL:Double:Sinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: sinh(s2/s1)/s1**2
+DEAL:Double:Sinh:Symbolic computation: Optimization::Sinh: About to optimise... done.
+DEAL:Double:Sinh:Symbolic computation: Substitution::Sinh: About to substitute... done.
 DEAL:Double:Sinh:Symbolic computation: Substitution::subs_psi: 0.7171584610
 DEAL:Double:Sinh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.2734612400
 DEAL:Double:Sinh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.4101918600
@@ -697,6 +789,8 @@ DEAL:Double:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: s2**
 DEAL:Double:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -sinh(s2/s1)/s1**2 - s2*cosh(s2/s1)/s1**3
 DEAL:Double:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -sinh(s2/s1)/s1**2 - s2*cosh(s2/s1)/s1**3
 DEAL:Double:Cosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: cosh(s2/s1)/s1**2
+DEAL:Double:Cosh:Symbolic computation: Optimization::Cosh: About to optimise... done.
+DEAL:Double:Cosh:Symbolic computation: Substitution::Cosh: About to substitute... done.
 DEAL:Double:Cosh:Symbolic computation: Substitution::subs_psi: 1.230575580
 DEAL:Double:Cosh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1593685469
 DEAL:Double:Cosh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2390528203
@@ -713,6 +807,8 @@ DEAL:Double:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*s2
 DEAL:Double:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -(1 - tanh(s2/s1)**2)/s1**2 + 2*s2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**3
 DEAL:Double:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -(1 - tanh(s2/s1)**2)/s1**2 + 2*s2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**3
 DEAL:Double:Tanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -2*tanh(s2/s1)*(1 - tanh(s2/s1)**2)/s1**2
+DEAL:Double:Tanh:Symbolic computation: Optimization::Tanh: About to optimise... done.
+DEAL:Double:Tanh:Symbolic computation: Substitution::Tanh: About to substitute... done.
 DEAL:Double:Tanh:Symbolic computation: Substitution::subs_psi: 0.5827829453
 DEAL:Double:Tanh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.1467475641
 DEAL:Double:Tanh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.2201213462
@@ -729,6 +825,8 @@ DEAL:Double:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -1.
 DEAL:Double:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -1.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -1.0*1/(s1**4*s2**4*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 1.0*1/(s1**2*s2**2*sqrt(1 + 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ASinh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -1.0*1/(s1**3*s2**5*(1 + 1.0*1/(s1**2*s2**2))**(3/2)) + 2.0*1/(s1*s2**3*sqrt(1 + 1.0*1/(s1**2*s2**2)))
+DEAL:Double:ASinh:Symbolic computation: Optimization::ASinh: About to optimise... done.
+DEAL:Double:ASinh:Symbolic computation: Substitution::ASinh: About to substitute... done.
 DEAL:Double:ASinh:Symbolic computation: Substitution::subs_psi: 0.1659045503
 DEAL:Double:ASinh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05479966244
 DEAL:Double:ASinh:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08219949365
@@ -745,6 +843,8 @@ DEAL:Double:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: -s1
 DEAL:Double:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -s1**2*s2**2/(-1 + s1**2*s2**2)**(3/2) + (-1 + s1**2*s2**2)**(-1/2)
 DEAL:Double:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -s1**2*s2**2/(-1 + s1**2*s2**2)**(3/2) + (-1 + s1**2*s2**2)**(-1/2)
 DEAL:Double:ACosh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -s1**3*s2/(-1 + s1**2*s2**2)**(3/2)
+DEAL:Double:ACosh:Symbolic computation: Optimization::ACosh: About to optimise... done.
+DEAL:Double:ACosh:Symbolic computation: Substitution::ACosh: About to substitute... done.
 DEAL:Double:ACosh:Symbolic computation: Substitution::subs_psi: 2.477888730
 DEAL:Double:ACosh:Symbolic computation: Substitution::subs_dpsi_ds0: 0.3380617019
 DEAL:Double:ACosh:Symbolic computation: Substitution::subs_dpsi_ds1: 0.5070925528
@@ -761,6 +861,8 @@ DEAL:Double:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2.0
 DEAL:Double:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 2.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 2.0*1/(s1**4*s2**4*(1 - 1.0*1/(s1**2*s2**2))**2) + 1.0*1/(s1**2*s2**2*(1 - 1.0*1/(s1**2*s2**2)))
 DEAL:Double:ATanh:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 2.0*1/(s1**3*s2**5*(1 - 1.0*1/(s1**2*s2**2))**2) + 2.0*1/(s1*s2**3*(1 - 1.0*1/(s1**2*s2**2)))
+DEAL:Double:ATanh:Symbolic computation: Optimization::ATanh: About to optimise... done.
+DEAL:Double:ATanh:Symbolic computation: Substitution::ATanh: About to substitute... done.
 DEAL:Double:ATanh:Symbolic computation: Substitution::subs_psi: 0.1682361183
 DEAL:Double:ATanh:Symbolic computation: Substitution::subs_dpsi_ds0: -0.05714285714
 DEAL:Double:ATanh:Symbolic computation: Substitution::subs_dpsi_ds1: -0.08571428571
@@ -777,6 +879,8 @@ DEAL:Double:Erf:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*erf
 DEAL:Double:Erf:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: 4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Double:Erf:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: 4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Double:Erf:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: -4*exp(-s2**2)*s1**2*s2/sqrt(pi)
+DEAL:Double:Erf:Symbolic computation: Optimization::Erf: About to optimise... done.
+DEAL:Double:Erf:Symbolic computation: Substitution::Erf: About to substitute... done.
 DEAL:Double:Erf:Symbolic computation: Substitution::subs_psi: 8.957900385
 DEAL:Double:Erf:Symbolic computation: Substitution::subs_dpsi_ds0: 5.971933590
 DEAL:Double:Erf:Symbolic computation: Substitution::subs_dpsi_ds1: 0.1860028682
@@ -793,6 +897,8 @@ DEAL:Double:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds0_ds0: 2*er
 DEAL:Double:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds1_ds0: -4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Double:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds0_ds1: -4*exp(-s2**2)*s1/sqrt(pi)
 DEAL:Double:Erfc:Symbolic computation: Differentiation::symb_d2psi_ds1_ds1: 4*exp(-s2**2)*s1**2*s2/sqrt(pi)
+DEAL:Double:Erfc:Symbolic computation: Optimization::Erfc: About to optimise... done.
+DEAL:Double:Erfc:Symbolic computation: Substitution::Erfc: About to substitute... done.
 DEAL:Double:Erfc:Symbolic computation: Substitution::subs_psi: 0.04209961483
 DEAL:Double:Erfc:Symbolic computation: Substitution::subs_dpsi_ds0: 0.02806640989
 DEAL:Double:Erfc:Symbolic computation: Substitution::subs_dpsi_ds1: -0.1860028682
diff --git a/tests/symengine/batch_optimizer_02_1a.output.release b/tests/symengine/batch_optimizer_02_1a.output.release
new file mode 100644 (file)
index 0000000..2006ced
--- /dev/null
@@ -0,0 +1,36 @@
+
+DEAL:Scalar::Scalar values
+DEAL:Scalar:Float::--------------------------------------------------------------------------------
+DEAL:Scalar:Float::Scalar
+DEAL:Scalar:Float::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Float::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Float:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Float:Substitution::evaluation:   s: 11.28968525  ds_dx: 2.440623999
+DEAL:Scalar:Float:Optimisation + substitution::evaluation:   s: 11.28968525  ds_dx: 2.440623999
+DEAL:Scalar:Float::OK
+DEAL:Scalar:Double::--------------------------------------------------------------------------------
+DEAL:Scalar:Double::Scalar
+DEAL:Scalar:Double::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Double::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Double:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Double:Substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
+DEAL:Scalar:Double:Optimisation + substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
+DEAL:Scalar:Double::OK
+DEAL:Tensor::Tensor values
+DEAL:Tensor:Float::--------------------------------------------------------------------------------
+DEAL:Tensor:Float::Tensor: Dim: 2
+DEAL:Tensor:Float::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Float::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Float:Substitution::substitute:   s: 733.14522057327  ds_dx: 1803.37506684305  ds_dy: 4021.22089917923 -2721.50016323639  -2073.1341017929 2076.40974260677  ds_dz: 1.67424568153838 1.97850294809826  1.97850294809826 2.6968452216256
+DEAL:Tensor:Float:Substitution::evaluation:   s: 733.1452026  ds_dx: 1803.375122  ds_dy: 4021.220947 -2721.500244  -2073.134033 2076.409668  ds_dz: 1.674245715 1.978502989  1.978502989 2.696845293
+DEAL:Tensor:Float:Optimisation + substitution::evaluation:   s: 733.1452026  ds_dx: 1803.375122  ds_dy: 4021.220947 -2721.500244  -2073.134033 2076.409668  ds_dz: 1.674245715 1.978502989  1.978502989 2.696845293
+DEAL:Tensor:Float::OK
+DEAL:Tensor:Double::--------------------------------------------------------------------------------
+DEAL:Tensor:Double::Tensor: Dim: 2
+DEAL:Tensor:Double::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Double::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Double:Substitution::substitute:   s: 733.145160811002  ds_dx: 1803.37488238402  ds_dy: 4021.22068318209 -2721.5000393195  -2073.13387957114 2076.40954281835  ds_dz: 1.67425190103196 1.97851016739323  1.97851016739323 2.69685737072851
+DEAL:Tensor:Double:Substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
+DEAL:Tensor:Double:Optimisation + substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
+DEAL:Tensor:Double::OK
+DEAL::OK
diff --git a/tests/symengine/batch_optimizer_02_1b.output.debug b/tests/symengine/batch_optimizer_02_1b.output.debug
new file mode 100644 (file)
index 0000000..88daa96
--- /dev/null
@@ -0,0 +1,36 @@
+
+DEAL:Scalar::Scalar values
+DEAL:Scalar:Float::--------------------------------------------------------------------------------
+DEAL:Scalar:Float::Scalar
+DEAL:Scalar:Float::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Float::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Float:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Float:Substitution::evaluation:   s: 11.28968525  ds_dx: 2.440623999
+DEAL:Scalar:Float:Optimisation + substitution::evaluation:   s: 11.28968525  ds_dx: 2.440623999
+DEAL:Scalar:Float::OK
+DEAL:Scalar:Double::--------------------------------------------------------------------------------
+DEAL:Scalar:Double::Scalar
+DEAL:Scalar:Double::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Double::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Double:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Double:Substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
+DEAL:Scalar:Double:Optimisation + substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
+DEAL:Scalar:Double::OK
+DEAL:Tensor::Tensor values
+DEAL:Tensor:Float::--------------------------------------------------------------------------------
+DEAL:Tensor:Float::Tensor: Dim: 2
+DEAL:Tensor:Float::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Float::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Float:Substitution::substitute:   s: 733.14522057327  ds_dx: 1803.37506684305  ds_dy: 4021.22089917923 -2721.50016323639  -2073.1341017929 2076.40974260677  ds_dz: 1.67424568153838 1.97850294809826  1.97850294809826 2.6968452216256
+DEAL:Tensor:Float:Substitution::evaluation:   s: 733.1452026  ds_dx: 1803.375122  ds_dy: 4021.220947 -2721.500244  -2073.134033 2076.409668  ds_dz: 1.674245715 1.978502989  1.978502989 2.696845293
+DEAL:Tensor:Float:Optimisation + substitution::evaluation:   s: 733.1452026  ds_dx: 1803.375122  ds_dy: 4021.220703 -2721.500000  -2073.134033 2076.409668  ds_dz: 1.674248815 1.978506565  1.978506565 2.696851254
+DEAL:Tensor:Float::OK
+DEAL:Tensor:Double::--------------------------------------------------------------------------------
+DEAL:Tensor:Double::Tensor: Dim: 2
+DEAL:Tensor:Double::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Double::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Double:Substitution::substitute:   s: 733.145160811002  ds_dx: 1803.37488238402  ds_dy: 4021.22068318209 -2721.5000393195  -2073.13387957114 2076.40954281835  ds_dz: 1.67425190103196 1.97851016739323  1.97851016739323 2.69685737072851
+DEAL:Tensor:Double:Substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
+DEAL:Tensor:Double:Optimisation + substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
+DEAL:Tensor:Double::OK
+DEAL::OK
diff --git a/tests/symengine/batch_optimizer_02_2a.output.release b/tests/symengine/batch_optimizer_02_2a.output.release
new file mode 100644 (file)
index 0000000..2006ced
--- /dev/null
@@ -0,0 +1,36 @@
+
+DEAL:Scalar::Scalar values
+DEAL:Scalar:Float::--------------------------------------------------------------------------------
+DEAL:Scalar:Float::Scalar
+DEAL:Scalar:Float::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Float::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Float:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Float:Substitution::evaluation:   s: 11.28968525  ds_dx: 2.440623999
+DEAL:Scalar:Float:Optimisation + substitution::evaluation:   s: 11.28968525  ds_dx: 2.440623999
+DEAL:Scalar:Float::OK
+DEAL:Scalar:Double::--------------------------------------------------------------------------------
+DEAL:Scalar:Double::Scalar
+DEAL:Scalar:Double::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Double::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Double:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Double:Substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
+DEAL:Scalar:Double:Optimisation + substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
+DEAL:Scalar:Double::OK
+DEAL:Tensor::Tensor values
+DEAL:Tensor:Float::--------------------------------------------------------------------------------
+DEAL:Tensor:Float::Tensor: Dim: 2
+DEAL:Tensor:Float::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Float::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Float:Substitution::substitute:   s: 733.14522057327  ds_dx: 1803.37506684305  ds_dy: 4021.22089917923 -2721.50016323639  -2073.1341017929 2076.40974260677  ds_dz: 1.67424568153838 1.97850294809826  1.97850294809826 2.6968452216256
+DEAL:Tensor:Float:Substitution::evaluation:   s: 733.1452026  ds_dx: 1803.375122  ds_dy: 4021.220947 -2721.500244  -2073.134033 2076.409668  ds_dz: 1.674245715 1.978502989  1.978502989 2.696845293
+DEAL:Tensor:Float:Optimisation + substitution::evaluation:   s: 733.1452026  ds_dx: 1803.375122  ds_dy: 4021.220947 -2721.500244  -2073.134033 2076.409668  ds_dz: 1.674245715 1.978502989  1.978502989 2.696845293
+DEAL:Tensor:Float::OK
+DEAL:Tensor:Double::--------------------------------------------------------------------------------
+DEAL:Tensor:Double::Tensor: Dim: 2
+DEAL:Tensor:Double::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Double::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Double:Substitution::substitute:   s: 733.145160811002  ds_dx: 1803.37488238402  ds_dy: 4021.22068318209 -2721.5000393195  -2073.13387957114 2076.40954281835  ds_dz: 1.67425190103196 1.97851016739323  1.97851016739323 2.69685737072851
+DEAL:Tensor:Double:Substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
+DEAL:Tensor:Double:Optimisation + substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
+DEAL:Tensor:Double::OK
+DEAL::OK
diff --git a/tests/symengine/batch_optimizer_02_2b.output.debug b/tests/symengine/batch_optimizer_02_2b.output.debug
new file mode 100644 (file)
index 0000000..2006ced
--- /dev/null
@@ -0,0 +1,36 @@
+
+DEAL:Scalar::Scalar values
+DEAL:Scalar:Float::--------------------------------------------------------------------------------
+DEAL:Scalar:Float::Scalar
+DEAL:Scalar:Float::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Float::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Float:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Float:Substitution::evaluation:   s: 11.28968525  ds_dx: 2.440623999
+DEAL:Scalar:Float:Optimisation + substitution::evaluation:   s: 11.28968525  ds_dx: 2.440623999
+DEAL:Scalar:Float::OK
+DEAL:Scalar:Double::--------------------------------------------------------------------------------
+DEAL:Scalar:Double::Scalar
+DEAL:Scalar:Double::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Double::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Double:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Double:Substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
+DEAL:Scalar:Double:Optimisation + substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
+DEAL:Scalar:Double::OK
+DEAL:Tensor::Tensor values
+DEAL:Tensor:Float::--------------------------------------------------------------------------------
+DEAL:Tensor:Float::Tensor: Dim: 2
+DEAL:Tensor:Float::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Float::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Float:Substitution::substitute:   s: 733.14522057327  ds_dx: 1803.37506684305  ds_dy: 4021.22089917923 -2721.50016323639  -2073.1341017929 2076.40974260677  ds_dz: 1.67424568153838 1.97850294809826  1.97850294809826 2.6968452216256
+DEAL:Tensor:Float:Substitution::evaluation:   s: 733.1452026  ds_dx: 1803.375122  ds_dy: 4021.220947 -2721.500244  -2073.134033 2076.409668  ds_dz: 1.674245715 1.978502989  1.978502989 2.696845293
+DEAL:Tensor:Float:Optimisation + substitution::evaluation:   s: 733.1452026  ds_dx: 1803.375122  ds_dy: 4021.220947 -2721.500244  -2073.134033 2076.409668  ds_dz: 1.674245715 1.978502989  1.978502989 2.696845293
+DEAL:Tensor:Float::OK
+DEAL:Tensor:Double::--------------------------------------------------------------------------------
+DEAL:Tensor:Double::Tensor: Dim: 2
+DEAL:Tensor:Double::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Double::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Double:Substitution::substitute:   s: 733.145160811002  ds_dx: 1803.37488238402  ds_dy: 4021.22068318209 -2721.5000393195  -2073.13387957114 2076.40954281835  ds_dz: 1.67425190103196 1.97851016739323  1.97851016739323 2.69685737072851
+DEAL:Tensor:Double:Substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
+DEAL:Tensor:Double:Optimisation + substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
+DEAL:Tensor:Double::OK
+DEAL::OK
index fc25b8942e44593c18c236e372dd1fb41cec360e..767b29b641369e252c6b09794fdd788820e658d5 100644 (file)
@@ -1,6 +1,36 @@
 
+DEAL:Scalar::Scalar values
+DEAL:Scalar:Float::--------------------------------------------------------------------------------
+DEAL:Scalar:Float::Scalar
+DEAL:Scalar:Float::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Float::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Float:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Float:Substitution::evaluation:   s: 11.28968525  ds_dx: 2.440623999
+DEAL:Scalar:Float:Optimisation + substitution::evaluation:   s: 11.28968430  ds_dx: 2.440624714
 DEAL:Scalar:Float::OK
+DEAL:Scalar:Double::--------------------------------------------------------------------------------
+DEAL:Scalar:Double::Scalar
+DEAL:Scalar:Double::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Double::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Double:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Double:Substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
+DEAL:Scalar:Double:Optimisation + substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
 DEAL:Scalar:Double::OK
+DEAL:Tensor::Tensor values
+DEAL:Tensor:Float::--------------------------------------------------------------------------------
+DEAL:Tensor:Float::Tensor: Dim: 2
+DEAL:Tensor:Float::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Float::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Float:Substitution::substitute:   s: 733.14522057327  ds_dx: 1803.37506684305  ds_dy: 4021.22089917923 -2721.50016323639  -2073.1341017929 2076.40974260677  ds_dz: 1.67424568153838 1.97850294809826  1.97850294809826 2.6968452216256
+DEAL:Tensor:Float:Substitution::evaluation:   s: 733.1452026  ds_dx: 1803.375122  ds_dy: 4021.220947 -2721.500244  -2073.134033 2076.409668  ds_dz: 1.674245715 1.978502989  1.978502989 2.696845293
+DEAL:Tensor:Float:Optimisation + substitution::evaluation:   s: 733.1455078  ds_dx: 1803.375977  ds_dy: 4021.222656 -2721.501221  -2073.134766 2076.410645  ds_dz: 1.674248934 1.978506684  1.978506684 2.696851492
 DEAL:Tensor:Float::OK
+DEAL:Tensor:Double::--------------------------------------------------------------------------------
+DEAL:Tensor:Double::Tensor: Dim: 2
+DEAL:Tensor:Double::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Double::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Double:Substitution::substitute:   s: 733.145160811002  ds_dx: 1803.37488238402  ds_dy: 4021.22068318209 -2721.5000393195  -2073.13387957114 2076.40954281835  ds_dz: 1.67425190103196 1.97851016739323  1.97851016739323 2.69685737072851
+DEAL:Tensor:Double:Substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
+DEAL:Tensor:Double:Optimisation + substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
 DEAL:Tensor:Double::OK
 DEAL::OK
index fc25b8942e44593c18c236e372dd1fb41cec360e..7c8c6e2294ab71e5a92aef098a1ac3b19a7036ac 100644 (file)
@@ -1,6 +1,36 @@
 
+DEAL:Scalar::Scalar values
+DEAL:Scalar:Float::--------------------------------------------------------------------------------
+DEAL:Scalar:Float::Scalar
+DEAL:Scalar:Float::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Float::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Float:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Float:Substitution::evaluation:   s: 11.28968525  ds_dx: 2.440623999
+DEAL:Scalar:Float:Optimisation + substitution::evaluation:   s: 11.28968525  ds_dx: 2.440624237
 DEAL:Scalar:Float::OK
+DEAL:Scalar:Double::--------------------------------------------------------------------------------
+DEAL:Scalar:Double::Scalar
+DEAL:Scalar:Double::symb_s: 1.5 + 1.5*y/log(x) + 2.0*x**y + sin(0.666666666666667*x*y**2)
+DEAL:Scalar:Double::symb_ds_dx: 2.0*x**(-1 + y)*y + 0.666666666666667*y**2*cos(0.666666666666667*x*y**2) - 1.5*y/(x*log(x)**2)
+DEAL:Scalar:Double:Substitution::substitute:   s: 11.2896853345375  ds_dx: 2.44062401525716
+DEAL:Scalar:Double:Substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
+DEAL:Scalar:Double:Optimisation + substitution::evaluation:   s: 11.28968533  ds_dx: 2.440624015
 DEAL:Scalar:Double::OK
+DEAL:Tensor::Tensor values
+DEAL:Tensor:Float::--------------------------------------------------------------------------------
+DEAL:Tensor:Float::Tensor: Dim: 2
+DEAL:Tensor:Float::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Float::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Float:Substitution::substitute:   s: 733.14522057327  ds_dx: 1803.37506684305  ds_dy: 4021.22089917923 -2721.50016323639  -2073.1341017929 2076.40974260677  ds_dz: 1.67424568153838 1.97850294809826  1.97850294809826 2.6968452216256
+DEAL:Tensor:Float:Substitution::evaluation:   s: 733.1452026  ds_dx: 1803.375122  ds_dy: 4021.220947 -2721.500244  -2073.134033 2076.409668  ds_dz: 1.674245715 1.978502989  1.978502989 2.696845293
+DEAL:Tensor:Float:Optimisation + substitution::evaluation:   s: 733.1455078  ds_dx: 1803.375977  ds_dy: 4021.222656 -2721.501465  -2073.135010 2076.410645  ds_dz: 1.674248934 1.978506684  1.978506684 2.696851492
 DEAL:Tensor:Float::OK
+DEAL:Tensor:Double::--------------------------------------------------------------------------------
+DEAL:Tensor:Double::Tensor: Dim: 2
+DEAL:Tensor:Double::symb_s: 1.5 + 1.5*(y_00*y_11 - y_01*y_10)*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 2.0*x**(y_00*y_11 - y_01*y_10) + sin(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01))
+DEAL:Tensor:Double::symb_ds_dy: -1.5*y_11 + 1.5*y_11*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_00*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_11 1.5*y_10 - 1.5*y_10*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_10  1.5*y_01 - 1.5*y_01*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_01*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) - 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_01 -1.5*y_00 + 1.5*y_00*log((z_00**2 + 2*z_01**2 + z_11**2)/(y_00*y_11 - y_01*y_10)) + 0.666666666666667*z_11*cos(0.666666666666667*(z_00*y_00 + z_11*y_11 + 1.0*(y_01 + y_10)*z_01)) + 2.0*x**(y_00*y_11 - y_01*y_10)*log(x)*y_00
+DEAL:Tensor:Double:Substitution::substitute:   s: 733.145160811002  ds_dx: 1803.37488238402  ds_dy: 4021.22068318209 -2721.5000393195  -2073.13387957114 2076.40954281835  ds_dz: 1.67425190103196 1.97851016739323  1.97851016739323 2.69685737072851
+DEAL:Tensor:Double:Substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
+DEAL:Tensor:Double:Optimisation + substitution::evaluation:   s: 733.1451608  ds_dx: 1803.374882  ds_dy: 4021.220683 -2721.500039  -2073.133880 2076.409543  ds_dz: 1.674251901 1.978510167  1.978510167 2.696857371
 DEAL:Tensor:Double::OK
 DEAL::OK
index a8bb4b95c5da18208b555e0540a8a473716ec805..c82b9a35cff233e90d253e54d3488af8984c7b67 100644 (file)
@@ -13,6 +13,7 @@ DEAL:Standard::HH[1][1][0][0]: 0.000000000
 DEAL:Standard::HH[1][1][0][1]: 0.000000000
 DEAL:Standard::HH[1][1][1][1]: -10.00000000
 DEAL::--------------------------------------------------
+DEAL::--------------------------------------------------
 DEAL:Explicit::psi: 53.40812106
 DEAL:Explicit::S[0][0]: 34.00000000
 DEAL:Explicit::S[0][1]: 0.000000000
@@ -27,6 +28,7 @@ DEAL:Explicit::HH[1][1][0][0]: -1.000000000
 DEAL:Explicit::HH[1][1][0][1]: 0.000000000
 DEAL:Explicit::HH[1][1][1][1]: 10.00000000
 DEAL::--------------------------------------------------
+DEAL::--------------------------------------------------
 DEAL:Implicit::psi: 26.44646255
 DEAL:Implicit::S[0][0]: 21.50000000
 DEAL:Implicit::S[0][1]: 0.000000000
index a8bb4b95c5da18208b555e0540a8a473716ec805..c82b9a35cff233e90d253e54d3488af8984c7b67 100644 (file)
@@ -13,6 +13,7 @@ DEAL:Standard::HH[1][1][0][0]: 0.000000000
 DEAL:Standard::HH[1][1][0][1]: 0.000000000
 DEAL:Standard::HH[1][1][1][1]: -10.00000000
 DEAL::--------------------------------------------------
+DEAL::--------------------------------------------------
 DEAL:Explicit::psi: 53.40812106
 DEAL:Explicit::S[0][0]: 34.00000000
 DEAL:Explicit::S[0][1]: 0.000000000
@@ -27,6 +28,7 @@ DEAL:Explicit::HH[1][1][0][0]: -1.000000000
 DEAL:Explicit::HH[1][1][0][1]: 0.000000000
 DEAL:Explicit::HH[1][1][1][1]: 10.00000000
 DEAL::--------------------------------------------------
+DEAL::--------------------------------------------------
 DEAL:Implicit::psi: 26.44646255
 DEAL:Implicit::S[0][0]: 21.50000000
 DEAL:Implicit::S[0][1]: 0.000000000
index 9f2261a1cb31be816f7112a259634df943a5af83..72eb781ce9a60e57c35e191cad3c96aa4241dc6c 100644 (file)
 
-DEAL::Timestep 1: 0.00 -0.000129 0.00
-DEAL::Timestep 2: 0.00 -0.000256 0.00
-DEAL::Timestep 3: 0.00 -0.000376 0.00
-DEAL::Timestep 4: 0.00 -0.000485 0.00
-DEAL::Timestep 5: 0.00 -0.000580 0.00
-DEAL::Timestep 6: 0.00 -0.000660 0.00
-DEAL::Timestep 7: 0.00 -0.000727 0.00
-DEAL::Timestep 8: 0.00 -0.000783 0.00
-DEAL::Timestep 9: 0.00 -0.000828 0.00
-DEAL::Timestep 10: 0.00 -0.000866 0.00
+DEAL::Grid:
+        Reference volume: 1.00000e-09
+DEAL::Triangulation:
+        Number of active cells: 8
+        Number of degrees of freedom: 97
+DEAL::    Setting up quadrature point data...
+DEAL::
+DEAL::Timestep 1 @ 0.100000s
+DEAL::                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
+DEAL:: 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 40       1.522e-06  1.000e+00  1.000e+00  0.000e+00  0.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
+DEAL:: 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 28       8.121e-04  8.399e-02  8.399e-02  1.528e-12  7.365e-08  8.331e-02  2.912e-02  8.331e-02  8.333e-02  
+DEAL:: 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 34       2.259e-06  2.189e-04  2.189e-04  3.750e-15  3.924e-10  3.699e-04  1.858e-04  3.699e-04  3.699e-04  
+DEAL:: 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 46       5.765e-12  4.523e-07  4.523e-07  3.141e-20  1.149e-14  2.854e-07  2.360e-07  2.854e-07  2.854e-07  
+DEAL:: 4   ASM_R  CONVERGED! 
+DEAL::Relative errors:
+DEAL::Displacement:    2.360e-07
+DEAL::Force:           1.153e-12
+DEAL::Dilatation:      6.277e-07
+DEAL::v / V_0: 1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 1: 0.000e+00 -1.287e-04 0.000e+00
+DEAL::
+DEAL::Timestep 2 @ 2.000e-01s
+DEAL::                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
+DEAL:: 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 37       2.358e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
+DEAL:: 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 26       1.216e-03  8.160e-02  8.160e-02  3.267e+13  4.545e+06  8.746e-02  3.290e-02  8.746e-02  8.748e-02  
+DEAL:: 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 36       1.534e-06  4.179e-04  4.179e-04  8.316e+10  2.772e+04  4.907e-04  2.486e-04  4.907e-04  4.907e-04  
+DEAL:: 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 44       4.637e-11  3.087e-07  3.087e-07  1.609e+06  1.371e+00  2.458e-07  1.474e-07  2.458e-07  2.458e-07  
+DEAL:: 4   ASM_R  CONVERGED! 
+DEAL::Relative errors:
+DEAL::Displacement:    1.474e-07
+DEAL::Force:           9.274e-12
+DEAL::Dilatation:      1.242e-06
+DEAL::v / V_0: 1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 2: 0.000e+00 -2.559e-04 0.000e+00
+DEAL::
+DEAL::Timestep 3 @ 3.000e-01s
+DEAL::                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
+DEAL:: 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 37       4.596e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
+DEAL:: 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 26       1.352e-03  8.278e-02  8.278e-02  4.057e+13  3.280e+06  8.213e-02  4.290e-02  8.213e-02  8.215e-02  
+DEAL:: 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 36       1.736e-06  6.361e-04  6.361e-04  1.008e+11  1.817e+04  6.209e-04  3.569e-04  6.209e-04  6.209e-04  
+DEAL:: 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 44       6.946e-11  3.533e-07  3.533e-07  4.325e+06  1.370e+00  2.813e-07  1.869e-07  2.813e-07  2.814e-07  
+DEAL:: 4   ASM_R  CONVERGED! 
+DEAL::Relative errors:
+DEAL::Displacement:    1.869e-07
+DEAL::Force:           1.389e-11
+DEAL::Dilatation:      1.831e-06
+DEAL::v / V_0: 1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 3: 0.000e+00 -3.761e-04 0.000e+00
+DEAL::
+DEAL::Timestep 4 @ 4.000e-01s
+DEAL::                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
+DEAL:: 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 37       4.628e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
+DEAL:: 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 27       6.091e-04  8.865e-02  8.865e-02  3.667e+13  1.810e+06  6.839e-02  5.409e-02  6.839e-02  6.840e-02  
+DEAL:: 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 35       4.266e-06  1.001e-03  1.001e-03  1.159e+11  7.046e+03  8.439e-04  5.351e-04  8.439e-04  8.439e-04  
+DEAL:: 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 43       2.896e-10  8.669e-07  8.669e-07  1.057e+07  1.402e+00  9.019e-07  6.095e-07  9.019e-07  9.020e-07  
+DEAL:: 4   ASM_R  CONVERGED! 
+DEAL::Relative errors:
+DEAL::Displacement:    6.095e-07
+DEAL::Force:           5.792e-11
+DEAL::Dilatation:      2.385e-06
+DEAL::v / V_0: 1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 4: 0.000e+00 -4.848e-04 0.000e+00
+DEAL::
+DEAL::Timestep 5 @ 5.000e-01s
+DEAL::                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
+DEAL:: 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 37       3.395e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
+DEAL:: 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 27       5.273e-04  9.675e-02  9.675e-02  1.860e+12  1.313e+06  5.249e-02  6.308e-02  5.249e-02  5.250e-02  
+DEAL:: 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 35       3.278e-06  1.421e-03  1.421e-03  7.998e+09  2.845e+03  1.054e-03  7.475e-04  1.054e-03  1.054e-03  
+DEAL:: 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 43       2.584e-10  6.849e-07  6.849e-07  9.940e+05  1.359e+00  9.631e-07  4.771e-07  9.631e-07  9.632e-07  
+DEAL:: 4   ASM_R  CONVERGED! 
+DEAL::Relative errors:
+DEAL::Displacement:    4.771e-07
+DEAL::Force:           5.167e-11
+DEAL::Dilatation:      2.902e-06
+DEAL::v / V_0: 1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 5: 0.000e+00 -5.797e-04 0.000e+00
+DEAL::
+DEAL::Timestep 6 @ 6.000e-01s
+DEAL::                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
+DEAL:: 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 37       3.664e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
+DEAL:: 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 27       4.474e-04  1.037e-01  1.037e-01  3.016e+12  1.755e+06  4.051e-02  6.866e-02  4.051e-02  4.052e-02  
+DEAL:: 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 35       3.007e-06  1.719e-03  1.719e-03  1.549e+10  2.116e+03  1.180e-03  8.919e-04  1.180e-03  1.180e-03  
+DEAL:: 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 43       4.359e-10  7.312e-07  7.312e-07  2.205e+06  1.985e+00  1.036e-06  3.333e-07  1.036e-06  1.036e-06  
+DEAL:: 4   ASM_R  CONVERGED! 
+DEAL::Relative errors:
+DEAL::Displacement:    3.333e-07
+DEAL::Force:           8.719e-11
+DEAL::Dilatation:      3.380e-06
+DEAL::v / V_0: 1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 6: 0.000e+00 -6.602e-04 0.000e+00
+DEAL::
+DEAL::Timestep 7 @ 7.000e-01s
+DEAL::                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
+DEAL:: 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 38       4.601e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
+DEAL:: 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 27       5.286e-04  1.073e-01  1.073e-01  6.122e+12  1.165e+06  3.333e-02  7.079e-02  3.333e-02  3.334e-02  
+DEAL:: 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 33       1.854e-06  1.847e-03  1.847e-03  3.318e+10  1.056e+03  1.124e-03  9.342e-04  1.124e-03  1.124e-03  
+DEAL:: 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 44       4.945e-10  5.880e-07  5.880e-07  4.659e+06  1.215e+00  4.851e-07  3.729e-07  4.851e-07  4.851e-07  
+DEAL:: 4   ASM_R  CONVERGED! 
+DEAL::Relative errors:
+DEAL::Displacement:    3.729e-07
+DEAL::Force:           9.889e-11
+DEAL::Dilatation:      3.819e-06
+DEAL::v / V_0: 1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 7: 0.000e+00 -7.274e-04 0.000e+00
+DEAL::
+DEAL::Timestep 8 @ 8.000e-01s
+DEAL::                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
+DEAL:: 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 37       2.070e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
+DEAL:: 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 29       2.533e-04  1.070e-01  1.070e-01  1.020e+13  1.349e+06  2.834e-02  7.005e-02  2.834e-02  2.835e-02  
+DEAL:: 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 33       1.788e-06  1.769e-03  1.769e-03  5.342e+10  9.657e+02  9.171e-04  8.799e-04  9.171e-04  9.171e-04  
+DEAL:: 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 47       8.654e-11  5.626e-07  5.626e-07  6.449e+06  1.313e+00  4.618e-07  3.937e-07  4.618e-07  4.619e-07  
+DEAL:: 4   ASM_R  CONVERGED! 
+DEAL::Relative errors:
+DEAL::Displacement:    3.937e-07
+DEAL::Force:           1.730e-11
+DEAL::Dilatation:      4.221e-06
+DEAL::v / V_0: 1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 8: 0.000e+00 -7.829e-04 0.000e+00
+DEAL::
+DEAL::Timestep 9 @ 9.000e-01s
+DEAL::                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
+DEAL:: 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 36       4.773e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
+DEAL:: 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 29       2.663e-04  1.036e-01  1.036e-01  1.231e+13  1.190e+06  2.394e-02  6.727e-02  2.394e-02  2.395e-02  
+DEAL:: 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 35       1.501e-06  1.566e-03  1.566e-03  5.876e+10  6.105e+02  6.930e-04  7.695e-04  6.930e-04  6.930e-04  
+DEAL:: 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 47       6.536e-11  4.212e-07  4.212e-07  5.691e+06  1.357e+00  6.387e-07  2.896e-07  6.387e-07  6.388e-07  
+DEAL:: 4   ASM_R  CONVERGED! 
+DEAL::Relative errors:
+DEAL::Displacement:    2.896e-07
+DEAL::Force:           1.307e-11
+DEAL::Dilatation:      4.588e-06
+DEAL::v / V_0: 1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 9: 0.000e+00 -8.285e-04 0.000e+00
+DEAL::
+DEAL::Timestep 10 @ 1.000e+00s
+DEAL::                 SOLVER STEP                   |  LIN_IT   LIN_RES    RES_NORM     RES_U     RES_P      RES_J     NU_NORM      NU_U       NU_P       NU_J 
+DEAL:: 0   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 37       2.601e-06  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  1.000e+00  
+DEAL:: 1   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 29       2.712e-04  9.816e-02  9.816e-02  7.184e+12  9.463e+05  1.994e-02  6.328e-02  1.994e-02  1.994e-02  
+DEAL:: 2   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 36       8.124e-07  1.309e-03  1.309e-03  3.014e+10  3.207e+02  4.941e-04  6.399e-04  4.941e-04  4.941e-04  
+DEAL:: 3   ASM_R  ASM_K  CST  ASM_SC  SLV  PP  UQPH  | 47       6.126e-11  2.591e-07  2.591e-07  2.245e+06  9.672e-01  1.873e-07  1.957e-07  1.873e-07  1.873e-07  
+DEAL:: 4   ASM_R  CONVERGED! 
+DEAL::Relative errors:
+DEAL::Displacement:    1.957e-07
+DEAL::Force:           1.225e-11
+DEAL::Dilatation:      4.924e-06
+DEAL::v / V_0: 1.000e-09 / 1.000e-09 = 1.000e+00
+DEAL::Timestep 10: 0.000e+00 -8.660e-04 0.000e+00

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.