From 9bba9104e766913d9fae13c6d85bb8102ecdb532 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Fr=C3=A9d=C3=A9ric=20Devernay?= Date: Wed, 16 Dec 2020 14:22:48 -0800 Subject: [PATCH 1/2] also build for long double if available --- CMakeLists.txt | 36 +++++++++++++++++++++++++++++------- README.md | 2 ++ 2 files changed, 31 insertions(+), 7 deletions(-) diff --git a/CMakeLists.txt b/CMakeLists.txt index 6c76950..5764ee5 100644 --- a/CMakeLists.txt +++ b/CMakeLists.txt @@ -50,14 +50,32 @@ set (cminpack_srcs set (cminpack_hdrs cminpack.h minpack.h) +include (CheckTypeSize) +check_type_size ("double" SIZEOF_DOUBLE) +check_type_size ("long double" SIZEOF_LONG_DOUBLE) + if (${CMINPACK_PRECISION} STREQUAL s) + message (STATUS "Building for single precision (float).") set (cminpack_libs cminpacks) elseif (${CMINPACK_PRECISION} STREQUAL d) + message (STATUS "Building for double precision (double).") set (cminpack_libs cminpack) elseif (${CMINPACK_PRECISION} STREQUAL ld) - set (cminpack_libs cminpackld) + if (${SIZEOF_LONG_DOUBLE} GREATER ${SIZEOF_DOUBLE}) + message (STATUS "Building for extended precision (long double).") + set (cminpack_libs cminpackld) + else () + # Microsoft Visual C++ for x86 makes long double a synonym for double + message (FATAL_ERROR "Compiler does not support long double, or long double is a synonym of double") + endif () elseif (${CMINPACK_PRECISION} STREQUAL all) - set (cminpack_libs cminpacks;cminpack) + if (${SIZEOF_LONG_DOUBLE} GREATER ${SIZEOF_DOUBLE}) + message (STATUS "Building for single precision (float), double precision (double) and extended precision (long double).") + set (cminpack_libs cminpacks;cminpack;cminpackld) + else () + message (STATUS "Building for single precision (float) and double precision (double).") + set (cminpack_libs cminpacks;cminpack) + endif () endif () if (BUILD_SHARED_LIBS) @@ -114,11 +132,15 @@ foreach (cminpack_lib ${cminpack_libs}) # Link with CBLAS library if requested if (USE_BLAS) - find_package (CBLAS) - if (CBLAS_FOUND) - target_link_libraries (${cminpack_lib} PUBLIC ${CBLAS_LIBRARIES}) - set_target_properties (${cminpack_lib} PROPERTIES LINK_FLAGS "${CBLAS_LINKER_FLAGS}") - target_compile_definitions (${cminpack_lib} PUBLIC USE_CBLAS) + if (${cminpack_lib} STREQUAL cminpackld) + message (WARNING "BLAS cannot be used for the extended precision version") + else () + find_package (CBLAS) + if (CBLAS_FOUND) + target_link_libraries (${cminpack_lib} PUBLIC ${CBLAS_LIBRARIES}) + set_target_properties (${cminpack_lib} PROPERTIES LINK_FLAGS "${CBLAS_LINKER_FLAGS}") + target_compile_definitions (${cminpack_lib} PUBLIC USE_CBLAS) + endif () endif () endif () diff --git a/README.md b/README.md index 1ebeef7..efc8deb 100644 --- a/README.md +++ b/README.md @@ -29,6 +29,8 @@ History ------ * version 1.3.8 (//2021): + - cmake now builds by default the single-, double-, and extended-precision versions #45 + - Avoid promoting to doubles in all operations for the single-precision version #47 * version 1.3.7 (09/12/2020): - Makefile cleanups #11 From dc2e4de4bb3675b7cd329cf84323a8e69c0fa3a6 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Fr=C3=A9d=C3=A9ric=20Devernay?= Date: Wed, 16 Dec 2020 14:41:35 -0800 Subject: [PATCH 2/2] fix long double tests --- README.md | 2 +- examples/ref/ldchkdrvc.ref | 311 ++++++++ examples/ref/ldhybdrvc.ref | 1204 +++++++++++++++++++++++++++++++ examples/ref/ldhyjdrvc.ref | 1204 +++++++++++++++++++++++++++++++ examples/ref/ldibmdpdrc.ref | 42 ++ examples/ref/ldlmddrvc.ref | 1151 +++++++++++++++++++++++++++++ examples/ref/ldlmfdrvc.ref | 1151 +++++++++++++++++++++++++++++ examples/ref/ldlmsdrvc.ref | 1148 +++++++++++++++++++++++++++++ examples/ref/ldtchkderc.ref | 22 + examples/ref/ldtchkderc_box.ref | 22 + examples/ref/ldtenormc.ref | 6 + examples/ref/ldtfdjac2c.ref | 29 + examples/ref/ldthybrd1c.ref | 7 + examples/ref/ldthybrdc.ref | 11 + examples/ref/ldthybrj1c.ref | 9 + examples/ref/ldthybrjc.ref | 14 + examples/ref/ldthybrjc_box.ref | 14 + examples/ref/ldtlmder1c.ref | 7 + examples/ref/ldtlmderc.ref | 16 + examples/ref/ldtlmderc_box.ref | 16 + examples/ref/ldtlmdif1c.ref | 7 + examples/ref/ldtlmdifc.ref | 14 + examples/ref/ldtlmstr1c.ref | 7 + examples/ref/ldtlmstrc.ref | 16 + examples/tenorm_.c | 5 +- examples/tenormc.c | 4 +- 26 files changed, 6434 insertions(+), 5 deletions(-) create mode 100644 examples/ref/ldchkdrvc.ref create mode 100644 examples/ref/ldhybdrvc.ref create mode 100644 examples/ref/ldhyjdrvc.ref create mode 100644 examples/ref/ldibmdpdrc.ref create mode 100644 examples/ref/ldlmddrvc.ref create mode 100644 examples/ref/ldlmfdrvc.ref create mode 100644 examples/ref/ldlmsdrvc.ref create mode 100644 examples/ref/ldtchkderc.ref create mode 100644 examples/ref/ldtchkderc_box.ref create mode 100644 examples/ref/ldtenormc.ref create mode 100644 examples/ref/ldtfdjac2c.ref create mode 100644 examples/ref/ldthybrd1c.ref create mode 100644 examples/ref/ldthybrdc.ref create mode 100644 examples/ref/ldthybrj1c.ref create mode 100644 examples/ref/ldthybrjc.ref create mode 100644 examples/ref/ldthybrjc_box.ref create mode 100644 examples/ref/ldtlmder1c.ref create mode 100644 examples/ref/ldtlmderc.ref create mode 100644 examples/ref/ldtlmderc_box.ref create mode 100644 examples/ref/ldtlmdif1c.ref create mode 100644 examples/ref/ldtlmdifc.ref create mode 100644 examples/ref/ldtlmstr1c.ref create mode 100644 examples/ref/ldtlmstrc.ref diff --git a/README.md b/README.md index efc8deb..314ecca 100644 --- a/README.md +++ b/README.md @@ -29,7 +29,7 @@ History ------ * version 1.3.8 (//2021): - - cmake now builds by default the single-, double-, and extended-precision versions #45 + - cmake now builds by default the single-, double-, and extended-precision versions #45 #48 - Avoid promoting to doubles in all operations for the single-precision version #47 * version 1.3.7 (09/12/2020): diff --git a/examples/ref/ldchkdrvc.ref b/examples/ref/ldchkdrvc.ref new file mode 100644 index 0000000..8a24e18 --- /dev/null +++ b/examples/ref/ldchkdrvc.ref @@ -0,0 +1,311 @@ + + + + problem 1 with dimension 2 is F + + + first function vector + + 2.0770000e+00 -2.8292900e+00 + + + function difference vector + + -3.5462622e-10 1.0526366e-08 + + + error vector + + 3.0078824e-02 9.2433674e-01 + + + + problem 2 with dimension 4 is F + + + first function vector + + -8.1070000e+00 -1.6859953e+00 1.8741610e+00 1.5952160e+01 + + + function difference vector + + 4.7260447e-09 -5.5515146e-10 -7.9065710e-10 1.0505208e-08 + + + error vector + + 1.0000000e+00 1.0000000e+00 4.6869594e-02 9.9092222e-01 + + + + problem 3 with dimension 2 is F + + + first function vector + + 1.0777100e+03 3.0019279e-01 + + + function difference vector + + 7.1037855e-07 -1.5595059e-10 + + + error vector + + 0.0000000e+00 0.0000000e+00 + + + + problem 4 with dimension 4 is T + + + first function vector + + -5.4127112e+03 -1.9649458e+03 -4.8718278e+03 -1.7769432e+03 + + + function difference vector + + 5.1311172e-06 1.1789170e-06 4.6181002e-06 1.0625044e-06 + + + error vector + + 9.6545728e-01 1.0000000e+00 9.7164712e-01 1.0000000e+00 + + + + problem 5 with dimension 3 is F + + + first function vector + + -5.0987696e+01 -1.1441658e+00 1.2300000e-01 + + + function difference vector + + 4.0500487e-10 -2.9159797e-09 4.0500487e-11 + + + error vector + + 1.0000000e+00 0.0000000e+00 1.0000000e+00 + + + + problem 6 with dimension 9 is F + + + first function vector + + -5.7930631e+00 -3.3390751e+01 -3.4683179e+01 -3.6368725e+01 -3.7947668e+01 + -3.9431326e+01 -4.0834243e+01 -4.2170791e+01 -4.3453299e+01 + + + function difference vector + + 9.9768537e-09 1.8235860e-08 2.3137328e-08 2.6977852e-08 3.0131870e-08 + 3.2820793e-08 3.5178894e-08 3.7293188e-08 3.9222860e-08 + + + error vector + + 0.0000000e+00 3.6958717e-02 4.8655594e-02 5.8177865e-02 6.5994814e-02 + 7.2648566e-02 7.8448546e-02 8.3592559e-02 8.8215713e-02 + + + + problem 7 with dimension 7 is F + + + first function vector + + 3.5142857e-02 -4.5634667e-02 2.1936396e-01 2.1302816e-01 2.5865828e-01 + 2.3088814e-01 5.9195986e-02 + + + function difference vector + + 3.4084382e-10 4.5526348e-10 5.2523568e-10 1.1820072e-09 2.1443179e-09 + 3.4838991e-09 4.6673111e-09 + + + error vector + + 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 + 0.0000000e+00 0.0000000e+00 + + + + problem 8 with dimension 10 is T + + + first function vector + + -5.6230000e+00 -5.3770000e+00 -5.6230000e+00 -5.3770000e+00 -5.6230000e+00 + -5.3770000e+00 -5.6230000e+00 -5.3770000e+00 -5.6230000e+00 -9.9928526e-01 + + + function difference vector + + 1.7704969e-09 1.8514979e-09 1.7704969e-09 1.8514979e-09 1.7704969e-09 + 1.8514979e-09 1.7704969e-09 1.8514979e-09 1.7704969e-09 2.3534345e-12 + + + error vector + + 9.9846705e-01 9.9261857e-01 9.9846705e-01 9.9261857e-01 9.9846705e-01 + 9.9261857e-01 9.9846705e-01 9.9261857e-01 9.9846705e-01 1.0000000e+00 + + + + problem 9 with dimension 10 is F + + + first function vector + + -3.8266208e-01 4.8185552e-01 -5.0497059e-01 4.8364359e-01 -5.0327111e-01 + 4.8731586e-01 -4.9982317e-01 4.9409864e-01 -4.9324735e-01 3.8308496e-01 + + + function difference vector + + 1.2760182e-10 -1.5641890e-10 1.6863172e-10 -1.5586223e-10 1.6922236e-10 + -1.5553038e-10 1.6982190e-10 -1.5572015e-10 1.4352145e-10 -6.2281313e-11 + + + error vector + + 2.8219889e-02 1.3376113e-01 2.0444304e-02 6.8047375e-02 1.3641636e-02 + 6.1770151e-02 1.5282455e-02 8.5364343e-02 2.6626937e-02 1.0200768e-01 + + + + problem 10 with dimension 10 is F + + + first function vector + + -1.6796677e-01 4.6728538e-02 -2.2043167e-01 1.7378706e-02 -2.2845451e-01 + 2.8983382e-02 -2.0089459e-01 6.9050605e-02 -1.5510284e-01 1.1399106e-01 + + + function difference vector + + 7.2803384e-11 1.8004944e-11 1.1962540e-10 5.2614146e-11 1.4146512e-10 + 6.1093737e-11 1.3625274e-10 4.1589837e-11 1.0264708e-10 2.0182886e-11 + + + error vector + + 0.0000000e+00 6.8687797e-03 0.0000000e+00 0.0000000e+00 0.0000000e+00 + 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 3.6993619e-02 + + + + problem 11 with dimension 10 is F + + + first function vector + + 1.4839305e-01 -4.6502477e-02 1.4892203e-01 3.0207835e-03 1.4945101e-01 + 5.2544044e-02 1.4997998e-01 1.0206730e-01 1.5050896e-01 1.5159056e-01 + + + function difference vector + + 7.2579293e-11 4.1193132e-11 7.2231127e-11 7.3671153e-11 7.1882961e-11 + 1.0614917e-10 7.1534795e-11 1.3862719e-10 7.1186630e-11 1.7110521e-10 + + + error vector + + 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 + 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 + + + + problem 12 with dimension 10 is F + + + first function vector + + -1.0878876e+05 -2.1757715e+05 -3.2636603e+05 -4.3515443e+05 -5.4394331e+05 + -6.5273170e+05 -7.6152058e+05 -8.7030898e+05 -9.7909786e+05 -1.0878862e+06 + + + function difference vector + + 4.9710803e-05 9.9421397e-05 1.4913183e-04 1.9884242e-04 2.4855286e-04 + 2.9826345e-04 3.4797389e-04 3.9768448e-04 4.4739493e-04 4.9710551e-04 + + + error vector + + 1.6738365e-02 1.6738288e-02 1.6738331e-02 1.6738301e-02 1.6738324e-02 + 1.6738305e-02 1.6738321e-02 1.6738307e-02 1.6738319e-02 1.6738308e-02 + + + + problem 13 with dimension 10 is T + + + first function vector + + -3.1372580e+00 1.9974200e-01 -2.2602580e+00 1.9974200e-01 -2.2602580e+00 + 1.9974200e-01 -2.2602580e+00 1.9974200e-01 -2.2602580e+00 -2.0462580e+00 + + + function difference vector + + 2.1927938e-09 7.7000843e-10 1.9040221e-09 7.7000843e-10 1.9040221e-09 + 7.7000843e-10 1.9040221e-09 7.7000843e-10 1.9040221e-09 1.5095539e-09 + + + error vector + + 1.0000000e+00 9.6830903e-01 9.7919715e-01 9.6830903e-01 9.7919715e-01 + 9.6830903e-01 9.7919715e-01 9.6830903e-01 9.7919715e-01 1.0000000e+00 + + + + problem 14 with dimension 10 is F + + + first function vector + + -8.2193683e+00 -4.4028887e+00 -8.2496263e+00 -4.4331467e+00 -8.2798843e+00 + -4.4634047e+00 -8.1720133e+00 -4.4634047e+00 -8.1720133e+00 -4.3252757e+00 + + + function difference vector + + 7.9522463e-09 4.8305583e-09 8.6307170e-09 5.5090291e-09 9.3091878e-09 + 6.1874998e-09 9.5269217e-09 6.1874998e-09 9.5269217e-09 5.7267630e-09 + + + error vector + + 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 + 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 0.0000000e+00 + summary of 14 tests of chkder + + nprob n status errmin errmax + + 1 2 F 3.0078824e-02 9.2433674e-01 + 2 4 F 4.6869594e-02 1.0000000e+00 + 3 2 F 0.0000000e+00 0.0000000e+00 + 4 4 T 9.6545728e-01 1.0000000e+00 + 5 3 F 0.0000000e+00 1.0000000e+00 + 6 9 F 0.0000000e+00 8.8215713e-02 + 7 7 F 0.0000000e+00 0.0000000e+00 + 8 10 T 9.9261857e-01 1.0000000e+00 + 9 10 F 1.3641636e-02 1.3376113e-01 + 10 10 F 0.0000000e+00 3.6993619e-02 + 11 10 F 0.0000000e+00 0.0000000e+00 + 12 10 F 1.6738288e-02 1.6738365e-02 + 13 10 T 9.6830903e-01 1.0000000e+00 + 14 10 F 0.0000000e+00 0.0000000e+00 diff --git a/examples/ref/ldhybdrvc.ref b/examples/ref/ldhybdrvc.ref new file mode 100644 index 0000000..040aa80 --- /dev/null +++ b/examples/ref/ldhybdrvc.ref @@ -0,0 +1,1204 @@ + + + + + problem 1 dimension 2 + + + initial l2 norm of the residuals 4.9193496e+00 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 16 + + number of Jacobian evaluations 3 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 1 dimension 2 + + + initial l2 norm of the residuals 1.3400631e+03 + + final l2 norm of the residuals 1.4094628e-17 + + number of function evaluations 7 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 1 dimension 2 + + + initial l2 norm of the residuals 1.4300005e+05 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 8 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 2 dimension 4 + + + initial l2 norm of the residuals 1.4662878e+01 + + final l2 norm of the residuals 5.6208452e-41 + + number of function evaluations 121 + + number of Jacobian evaluations 5 + + exit parameter 4 + + final approximate solution + + 3.1663427e-21 -3.1663427e-22 3.5902019e-21 3.5902019e-21 + + + + + problem 2 dimension 4 + + + initial l2 norm of the residuals 1.2709839e+03 + + final l2 norm of the residuals 3.4399657e-40 + + number of function evaluations 114 + + number of Jacobian evaluations 3 + + exit parameter 4 + + final approximate solution + + 7.5145392e-21 -7.5145392e-22 8.8971294e-21 8.8971294e-21 + + + + + problem 2 dimension 4 + + + initial l2 norm of the residuals 1.2688790e+05 + + final l2 norm of the residuals 2.3005507e-41 + + number of function evaluations 123 + + number of Jacobian evaluations 3 + + exit parameter 4 + + final approximate solution + + 4.1789519e-22 -4.1789519e-23 5.8356510e-22 5.8356510e-22 + + + + + problem 3 dimension 2 + + + initial l2 norm of the residuals 1.0654866e+00 + + final l2 norm of the residuals 1.1626313e-10 + + number of function evaluations 170 + + number of Jacobian evaluations 6 + + exit parameter 1 + + final approximate solution + + 1.0981593e-05 9.1061467e+00 + + + + + problem 3 dimension 2 + + + initial l2 norm of the residuals 1.0000000e+00 + + final l2 norm of the residuals 2.2204460e-16 + + number of function evaluations 11 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 1.0981593e-05 9.1061467e+00 + + + + + problem 4 dimension 4 + + + initial l2 norm of the residuals 8.5505574e+03 + + final l2 norm of the residuals 1.0617427e-11 + + number of function evaluations 87 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -9.6797402e-01 9.4713914e-01 -9.6951631e-01 9.5124767e-01 + + + + + problem 4 dimension 4 + + + initial l2 norm of the residuals 7.3498230e+06 + + final l2 norm of the residuals 4.0879794e-14 + + number of function evaluations 204 + + number of Jacobian evaluations 8 + + exit parameter 1 + + final approximate solution + + -9.6797402e-01 9.4713914e-01 -9.6951631e-01 9.5124767e-01 + + + + + problem 4 dimension 4 + + + initial l2 norm of the residuals 7.2730700e+09 + + final l2 norm of the residuals 1.4814296e-11 + + number of function evaluations 369 + + number of Jacobian evaluations 35 + + exit parameter 1 + + final approximate solution + + -9.6797402e-01 9.4713914e-01 -9.6951631e-01 9.5124767e-01 + + + + + problem 5 dimension 3 + + + initial l2 norm of the residuals 5.0000000e+01 + + final l2 norm of the residuals 2.7534517e-13 + + number of function evaluations 18 + + number of Jacobian evaluations 3 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 -1.6121109e-14 -9.2946876e-36 + + + + + problem 5 dimension 3 + + + initial l2 norm of the residuals 1.0295630e+02 + + final l2 norm of the residuals 3.7079701e-19 + + number of function evaluations 20 + + number of Jacobian evaluations 4 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 2.3297863e-20 1.5616214e-49 + + + + + problem 5 dimension 3 + + + initial l2 norm of the residuals 9.9126182e+02 + + final l2 norm of the residuals 1.6220693e-13 + + number of function evaluations 28 + + number of Jacobian evaluations 4 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.2499837e-15 5.0763888e-48 + + + + + problem 6 dimension 6 + + + initial l2 norm of the residuals 6.8485872e+01 + + final l2 norm of the residuals 5.6066258e-16 + + number of function evaluations 61 + + number of Jacobian evaluations 6 + + exit parameter 1 + + final approximate solution + + -1.5725086e-02 1.0124349e+00 -2.3299163e-01 1.2604301e+00 -1.5137289e+00 + 9.9299643e-01 + + + + + problem 6 dimension 6 + + + initial l2 norm of the residuals 3.5312586e+06 + + final l2 norm of the residuals 4.9884892e-15 + + number of function evaluations 140 + + number of Jacobian evaluations 7 + + exit parameter 1 + + final approximate solution + + -1.5725086e-02 1.0124349e+00 -2.3299163e-01 1.2604301e+00 -1.5137289e+00 + 9.9299643e-01 + + + + + problem 6 dimension 9 + + + initial l2 norm of the residuals 8.8789552e+01 + + final l2 norm of the residuals 1.4027649e-12 + + number of function evaluations 92 + + number of Jacobian evaluations 4 + + exit parameter 1 + + final approximate solution + + -1.5307037e-05 9.9978970e-01 1.4763964e-02 1.4634233e-01 1.0008211e+00 + -2.6177311e+00 4.1044032e+00 -3.1436123e+00 1.0526264e+00 + + + + + problem 6 dimension 9 + + + initial l2 norm of the residuals 1.0151080e+07 + + final l2 norm of the residuals 3.5259553e-02 + + number of function evaluations 313 + + number of Jacobian evaluations 13 + + exit parameter 4 + + final approximate solution + + 3.0955000e-01 1.0887169e+00 1.3364306e+00 -1.1286099e+01 6.0292577e+01 + -1.5749196e+02 2.2409538e+02 -1.6366396e+02 4.8983186e+01 + + + + + problem 7 dimension 5 + + + initial l2 norm of the residuals 2.2570657e-01 + + final l2 norm of the residuals 5.1051901e-14 + + number of function evaluations 12 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 8.3751256e-02 3.1272930e-01 5.0000000e-01 6.8727070e-01 9.1624874e-01 + + + + + problem 7 dimension 5 + + + initial l2 norm of the residuals 4.1172432e+06 + + final l2 norm of the residuals 1.0075086e-11 + + number of function evaluations 130 + + number of Jacobian evaluations 25 + + exit parameter 1 + + final approximate solution + + 8.3751256e-02 5.0000000e-01 3.1272930e-01 9.1624874e-01 6.8727070e-01 + + + + + problem 7 dimension 5 + + + initial l2 norm of the residuals 5.6361303e+11 + + final l2 norm of the residuals 6.1979157e-14 + + number of function evaluations 244 + + number of Jacobian evaluations 49 + + exit parameter 1 + + final approximate solution + + 5.0000000e-01 8.3751256e-02 3.1272930e-01 9.1624874e-01 6.8727070e-01 + + + + + problem 7 dimension 6 + + + initial l2 norm of the residuals 2.1547198e-01 + + final l2 norm of the residuals 3.7765156e-13 + + number of function evaluations 15 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 6.6876591e-02 3.6668230e-01 2.8874067e-01 7.1125933e-01 6.3331770e-01 + 9.3312341e-01 + + + + + problem 7 dimension 6 + + + initial l2 norm of the residuals 1.3079247e+08 + + final l2 norm of the residuals 7.9417332e-13 + + number of function evaluations 96 + + number of Jacobian evaluations 13 + + exit parameter 1 + + final approximate solution + + 9.3312341e-01 3.6668230e-01 6.6876591e-02 7.1125933e-01 2.8874067e-01 + 6.3331770e-01 + + + + + problem 7 dimension 6 + + + initial l2 norm of the residuals 1.8755789e+14 + + final l2 norm of the residuals 7.7595534e-14 + + number of function evaluations 136 + + number of Jacobian evaluations 21 + + exit parameter 1 + + final approximate solution + + 3.6668230e-01 6.3331770e-01 9.3312341e-01 6.6876591e-02 7.1125933e-01 + 2.8874067e-01 + + + + + problem 7 dimension 7 + + + initial l2 norm of the residuals 1.8376789e-01 + + final l2 norm of the residuals 1.3208654e-11 + + number of function evaluations 14 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 5.8069150e-02 2.3517161e-01 3.3804409e-01 5.0000000e-01 6.6195591e-01 + 7.6482839e-01 9.4193085e-01 + + + + + problem 7 dimension 7 + + + initial l2 norm of the residuals 4.2693282e+09 + + final l2 norm of the residuals 1.1808623e-11 + + number of function evaluations 328 + + number of Jacobian evaluations 43 + + exit parameter 1 + + final approximate solution + + 5.0000000e-01 3.3804409e-01 2.3517161e-01 5.8069150e-02 6.6195591e-01 + 9.4193085e-01 7.6482839e-01 + + + + + problem 7 dimension 7 + + + initial l2 norm of the residuals 6.4143166e+16 + + final l2 norm of the residuals 6.5823425e+11 + + number of function evaluations 57 + + number of Jacobian evaluations 7 + + exit parameter 4 + + final approximate solution + + -2.0419317e+00 -3.1262717e+01 2.8910888e+01 2.4452585e+01 2.7755867e+01 + 2.1278397e+01 1.6794390e+01 + + + + + problem 7 dimension 8 + + + initial l2 norm of the residuals 1.9651386e-01 + + final l2 norm of the residuals 6.4405121e-02 + + number of function evaluations 56 + + number of Jacobian evaluations 8 + + exit parameter 4 + + final approximate solution + + 4.9856400e-02 1.9863513e-01 2.6982882e-01 4.9927230e-01 5.0072770e-01 + 7.3017118e-01 8.0136487e-01 9.5014360e-01 + + + + + problem 7 dimension 9 + + + initial l2 norm of the residuals 1.6994993e-01 + + final l2 norm of the residuals 2.7932597e-12 + + number of function evaluations 26 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 4.4205346e-02 1.9949067e-01 2.3561911e-01 4.1604691e-01 5.0000000e-01 + 5.8395309e-01 7.6438089e-01 8.0050933e-01 9.5579465e-01 + + + + + problem 8 dimension 10 + + + initial l2 norm of the residuals 1.6530216e+01 + + final l2 norm of the residuals 5.4210109e-20 + + number of function evaluations 12 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 8 dimension 10 + + + initial l2 norm of the residuals 9.7656240e+06 + + final l2 norm of the residuals 3.4285483e-19 + + number of function evaluations 12 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 8 dimension 10 + + + initial l2 norm of the residuals 9.7656250e+16 + + final l2 norm of the residuals 8.3307375e-18 + + number of function evaluations 15 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 8 dimension 30 + + + initial l2 norm of the residuals 8.3476044e+01 + + final l2 norm of the residuals 6.4848524e-18 + + number of function evaluations 24 + + number of Jacobian evaluations 3 + + exit parameter 1 + + final approximate solution + + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 1.0673735e+00 + + + + + problem 8 dimension 40 + + + initial l2 norm of the residuals 1.2802636e+02 + + final l2 norm of the residuals 9.5808925e-17 + + number of function evaluations 13 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 9 dimension 10 + + + initial l2 norm of the residuals 2.8080582e-02 + + final l2 norm of the residuals 2.5152500e-15 + + number of function evaluations 6 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 9 dimension 10 + + + initial l2 norm of the residuals 5.2555258e-01 + + final l2 norm of the residuals 1.7259752e-13 + + number of function evaluations 9 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 9 dimension 10 + + + initial l2 norm of the residuals 1.0657390e+02 + + final l2 norm of the residuals 7.1342011e-12 + + number of function evaluations 44 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 10 dimension 1 + + + initial l2 norm of the residuals 1.2792969e-01 + + final l2 norm of the residuals 4.0657581e-20 + + number of function evaluations 6 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -1.5281388e-01 + + + + + problem 10 dimension 1 + + + initial l2 norm of the residuals 2.5625000e+00 + + final l2 norm of the residuals 4.0657581e-20 + + number of function evaluations 8 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -1.5281388e-01 + + + + + problem 10 dimension 1 + + + initial l2 norm of the residuals 8.3611719e+02 + + final l2 norm of the residuals 2.7105054e-20 + + number of function evaluations 16 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -1.5281388e-01 + + + + + problem 10 dimension 10 + + + initial l2 norm of the residuals 2.5182701e-01 + + final l2 norm of the residuals 5.0324391e-15 + + number of function evaluations 6 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 10 dimension 10 + + + initial l2 norm of the residuals 6.1168330e+00 + + final l2 norm of the residuals 2.1884220e-13 + + number of function evaluations 9 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 10 dimension 10 + + + initial l2 norm of the residuals 1.2693089e+03 + + final l2 norm of the residuals 3.0402481e-15 + + number of function evaluations 19 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 11 dimension 10 + + + initial l2 norm of the residuals 8.4117534e-02 + + final l2 norm of the residuals 5.3041556e-03 + + number of function evaluations 48 + + number of Jacobian evaluations 9 + + exit parameter 4 + + final approximate solution + + 5.5311057e-02 5.7008719e-02 5.8941652e-02 6.1171890e-02 6.3845908e-02 + 6.7098188e-02 2.0789713e-01 1.6420765e-01 8.6862485e-02 9.1085023e-02 + + + + + problem 11 dimension 10 + + + initial l2 norm of the residuals 2.0305195e+01 + + final l2 norm of the residuals 3.3326566e-13 + + number of function evaluations 36 + + number of Jacobian evaluations 5 + + exit parameter 1 + + final approximate solution + + 3.4396289e-02 3.5032316e-02 3.5719196e-02 3.6465224e-02 3.7280912e-02 + 3.8179863e-02 3.9180141e-02 4.0306503e-02 1.7972019e-01 1.5624088e-01 + + + + + problem 11 dimension 10 + + + initial l2 norm of the residuals 9.3369375e+01 + + final l2 norm of the residuals 1.3332552e-09 + + number of function evaluations 44 + + number of Jacobian evaluations 4 + + exit parameter 1 + + final approximate solution + + 1.8883952e+01 2.5167774e+01 1.8885275e+01 1.8886021e+01 1.8886837e+01 + 1.8887736e+01 1.8888736e+01 1.8889862e+01 1.9029276e+01 1.9005797e+01 + + + + + problem 12 dimension 10 + + + initial l2 norm of the residuals 2.2402135e+06 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 22 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 12 dimension 10 + + + initial l2 norm of the residuals 5.2234376e+07 + + final l2 norm of the residuals 1.8007939e-14 + + number of function evaluations 26 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 12 dimension 10 + + + initial l2 norm of the residuals 1.5923646e+11 + + final l2 norm of the residuals 6.8293667e-11 + + number of function evaluations 46 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 13 dimension 10 + + + initial l2 norm of the residuals 4.5825757e+00 + + final l2 norm of the residuals 7.6228700e-11 + + number of function evaluations 13 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -5.7072213e-01 -6.8180695e-01 -7.0221008e-01 -7.0551063e-01 -7.0490616e-01 + -7.0149661e-01 -6.9188932e-01 -6.6579651e-01 -5.9603511e-01 -4.1641226e-01 + + + + + problem 13 dimension 10 + + + initial l2 norm of the residuals 6.3910093e+02 + + final l2 norm of the residuals 2.3214215e-12 + + number of function evaluations 53 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -5.7072213e-01 -6.8180695e-01 -7.0221008e-01 -7.0551063e-01 -7.0490616e-01 + -7.0149661e-01 -6.9188932e-01 -6.6579651e-01 -5.9603511e-01 -4.1641226e-01 + + + + + problem 13 dimension 10 + + + initial l2 norm of the residuals 6.3337583e+04 + + final l2 norm of the residuals 3.2465706e-12 + + number of function evaluations 23 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -5.7072213e-01 -6.8180695e-01 -7.0221008e-01 -7.0551063e-01 -7.0490616e-01 + -7.0149661e-01 -6.9188932e-01 -6.6579651e-01 -5.9603511e-01 -4.1641226e-01 + + + + + problem 14 dimension 10 + + + initial l2 norm of the residuals 1.8973666e+01 + + final l2 norm of the residuals 8.2517941e-11 + + number of function evaluations 23 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.2830286e-01 -4.7659642e-01 -5.1965246e-01 -5.5809932e-01 -5.9250616e-01 + -6.2450368e-01 -6.2323947e-01 -6.2139384e-01 -6.2045360e-01 -5.8646927e-01 + + + + + problem 14 dimension 10 + + + initial l2 norm of the residuals 1.7130922e+04 + + final l2 norm of the residuals 1.3555351e-10 + + number of function evaluations 27 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -4.2830286e-01 -4.7659642e-01 -5.1965246e-01 -5.5809932e-01 -5.9250616e-01 + -6.2450368e-01 -6.2323947e-01 -6.2139384e-01 -6.2045360e-01 -5.8646927e-01 + + + + + problem 14 dimension 10 + + + initial l2 norm of the residuals 1.5949860e+07 + + final l2 norm of the residuals 8.7598526e-11 + + number of function evaluations 39 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -4.2830286e-01 -4.7659642e-01 -5.1965246e-01 -5.5809932e-01 -5.9250616e-01 + -6.2450368e-01 -6.2323947e-01 -6.2139384e-01 -6.2045360e-01 -5.8646927e-01 + summary of 55 calls to hybrd1 + + nprob n nfev njev info final l2 norm + + 1 2 16 3 1 0.0000000e+00 + 1 2 7 1 1 1.4094628e-17 + 1 2 8 1 1 0.0000000e+00 + 2 4 121 5 4 5.6208452e-41 + 2 4 114 3 4 3.4399657e-40 + 2 4 123 3 4 2.3005507e-41 + 3 2 170 6 1 1.1626313e-10 + 3 2 11 1 1 2.2204460e-16 + 4 4 87 2 1 1.0617427e-11 + 4 4 204 8 1 4.0879794e-14 + 4 4 369 35 1 1.4814296e-11 + 5 3 18 3 1 2.7534517e-13 + 5 3 20 4 1 3.7079701e-19 + 5 3 28 4 1 1.6220693e-13 + 6 6 61 6 1 5.6066258e-16 + 6 6 140 7 1 4.9884892e-15 + 6 9 92 4 1 1.4027649e-12 + 6 9 313 13 4 3.5259553e-02 + 7 5 12 1 1 5.1051901e-14 + 7 5 130 25 1 1.0075086e-11 + 7 5 244 49 1 6.1979157e-14 + 7 6 15 2 1 3.7765156e-13 + 7 6 96 13 1 7.9417332e-13 + 7 6 136 21 1 7.7595534e-14 + 7 7 14 1 1 1.3208654e-11 + 7 7 328 43 1 1.1808623e-11 + 7 7 57 7 4 6.5823425e+11 + 7 8 56 8 4 6.4405121e-02 + 7 9 26 2 1 2.7932597e-12 + 8 10 12 2 1 5.4210109e-20 + 8 10 12 2 1 3.4285483e-19 + 8 10 15 2 1 8.3307375e-18 + 8 30 24 3 1 6.4848524e-18 + 8 40 13 2 1 9.5808925e-17 + 9 10 6 1 1 2.5152500e-15 + 9 10 9 1 1 1.7259752e-13 + 9 10 44 1 1 7.1342011e-12 + 10 1 6 1 1 4.0657581e-20 + 10 1 8 1 1 4.0657581e-20 + 10 1 16 1 1 2.7105054e-20 + 10 10 6 1 1 5.0324391e-15 + 10 10 9 1 1 2.1884220e-13 + 10 10 19 2 1 3.0402481e-15 + 11 10 48 9 4 5.3041556e-03 + 11 10 36 5 1 3.3326566e-13 + 11 10 44 4 1 1.3332552e-09 + 12 10 22 1 1 0.0000000e+00 + 12 10 26 1 1 1.8007939e-14 + 12 10 46 2 1 6.8293667e-11 + 13 10 13 1 1 7.6228700e-11 + 13 10 53 1 1 2.3214215e-12 + 13 10 23 2 1 3.2465706e-12 + 14 10 23 1 1 8.2517941e-11 + 14 10 27 2 1 1.3555351e-10 + 14 10 39 2 1 8.7598526e-11 diff --git a/examples/ref/ldhyjdrvc.ref b/examples/ref/ldhyjdrvc.ref new file mode 100644 index 0000000..7dd5801 --- /dev/null +++ b/examples/ref/ldhyjdrvc.ref @@ -0,0 +1,1204 @@ + + + + + problem 1 dimension 2 + + + initial l2 norm of the residuals 4.9193496e+00 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 16 + + number of jacobian evaluations 3 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 1 dimension 2 + + + initial l2 norm of the residuals 1.3400631e+03 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 7 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 1 dimension 2 + + + initial l2 norm of the residuals 1.4300005e+05 + + final l2 norm of the residuals 2.6237693e-16 + + number of function evaluations 7 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 2 dimension 4 + + + initial l2 norm of the residuals 1.4662878e+01 + + final l2 norm of the residuals 1.0965660e-43 + + number of function evaluations 126 + + number of jacobian evaluations 5 + + exit parameter 4 + + final approximate solution + + 9.2134084e-23 -9.2134084e-24 1.6022069e-22 1.6022069e-22 + + + + + problem 2 dimension 4 + + + initial l2 norm of the residuals 1.2709839e+03 + + final l2 norm of the residuals 6.0574582e-42 + + number of function evaluations 125 + + number of jacobian evaluations 4 + + exit parameter 4 + + final approximate solution + + 2.9423593e-22 -2.9423593e-23 -9.5371000e-22 -9.5371000e-22 + + + + + problem 2 dimension 4 + + + initial l2 norm of the residuals 1.2688790e+05 + + final l2 norm of the residuals 3.6871121e-41 + + number of function evaluations 169 + + number of jacobian evaluations 9 + + exit parameter 4 + + final approximate solution + + 4.4226549e-21 -4.4226549e-22 1.0331732e-21 1.0331732e-21 + + + + + problem 3 dimension 2 + + + initial l2 norm of the residuals 1.0654866e+00 + + final l2 norm of the residuals 1.1585631e-10 + + number of function evaluations 170 + + number of jacobian evaluations 6 + + exit parameter 1 + + final approximate solution + + 1.0981593e-05 9.1061467e+00 + + + + + problem 3 dimension 2 + + + initial l2 norm of the residuals 1.0000000e+00 + + final l2 norm of the residuals 1.0842022e-19 + + number of function evaluations 11 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 1.0981593e-05 9.1061467e+00 + + + + + problem 4 dimension 4 + + + initial l2 norm of the residuals 8.5505574e+03 + + final l2 norm of the residuals 1.0569600e-11 + + number of function evaluations 87 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -9.6797402e-01 9.4713914e-01 -9.6951631e-01 9.5124767e-01 + + + + + problem 4 dimension 4 + + + initial l2 norm of the residuals 7.3498230e+06 + + final l2 norm of the residuals 4.3101277e-14 + + number of function evaluations 204 + + number of jacobian evaluations 8 + + exit parameter 1 + + final approximate solution + + -9.6797402e-01 9.4713914e-01 -9.6951631e-01 9.5124767e-01 + + + + + problem 4 dimension 4 + + + initial l2 norm of the residuals 7.2730700e+09 + + final l2 norm of the residuals 3.4779835e-11 + + number of function evaluations 386 + + number of jacobian evaluations 35 + + exit parameter 1 + + final approximate solution + + -9.6797402e-01 9.4713914e-01 -9.6951631e-01 9.5124767e-01 + + + + + problem 5 dimension 3 + + + initial l2 norm of the residuals 5.0000000e+01 + + final l2 norm of the residuals 2.7534492e-13 + + number of function evaluations 18 + + number of jacobian evaluations 3 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 -1.6121092e-14 0.0000000e+00 + + + + + problem 5 dimension 3 + + + initial l2 norm of the residuals 1.0295630e+02 + + final l2 norm of the residuals 1.4599744e-12 + + number of function evaluations 21 + + number of jacobian evaluations 4 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 -9.0059944e-14 0.0000000e+00 + + + + + problem 5 dimension 3 + + + initial l2 norm of the residuals 9.9126182e+02 + + final l2 norm of the residuals 4.7530476e-13 + + number of function evaluations 38 + + number of jacobian evaluations 8 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 2.1305562e-14 0.0000000e+00 + + + + + problem 6 dimension 6 + + + initial l2 norm of the residuals 6.8485872e+01 + + final l2 norm of the residuals 5.7711673e-16 + + number of function evaluations 61 + + number of jacobian evaluations 6 + + exit parameter 1 + + final approximate solution + + -1.5725086e-02 1.0124349e+00 -2.3299163e-01 1.2604301e+00 -1.5137289e+00 + 9.9299643e-01 + + + + + problem 6 dimension 6 + + + initial l2 norm of the residuals 3.5312586e+06 + + final l2 norm of the residuals 8.2109861e-16 + + number of function evaluations 140 + + number of jacobian evaluations 6 + + exit parameter 1 + + final approximate solution + + -1.5725086e-02 1.0124349e+00 -2.3299163e-01 1.2604301e+00 -1.5137289e+00 + 9.9299643e-01 + + + + + problem 6 dimension 9 + + + initial l2 norm of the residuals 8.8789552e+01 + + final l2 norm of the residuals 4.7333350e-13 + + number of function evaluations 93 + + number of jacobian evaluations 4 + + exit parameter 1 + + final approximate solution + + -1.5307037e-05 9.9978970e-01 1.4763964e-02 1.4634233e-01 1.0008211e+00 + -2.6177311e+00 4.1044032e+00 -3.1436123e+00 1.0526264e+00 + + + + + problem 6 dimension 9 + + + initial l2 norm of the residuals 1.0151080e+07 + + final l2 norm of the residuals 3.3676416e-02 + + number of function evaluations 301 + + number of jacobian evaluations 13 + + exit parameter 4 + + final approximate solution + + 3.1018658e-01 1.0889663e+00 1.3539024e+00 -1.1469368e+01 6.1188578e+01 + -1.5974966e+02 2.2719004e+02 -1.6585335e+02 4.9614166e+01 + + + + + problem 7 dimension 5 + + + initial l2 norm of the residuals 2.2570657e-01 + + final l2 norm of the residuals 6.0446772e-13 + + number of function evaluations 10 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 8.3751256e-02 3.1272930e-01 5.0000000e-01 6.8727070e-01 9.1624874e-01 + + + + + problem 7 dimension 5 + + + initial l2 norm of the residuals 4.1172432e+06 + + final l2 norm of the residuals 8.7599009e-13 + + number of function evaluations 142 + + number of jacobian evaluations 25 + + exit parameter 1 + + final approximate solution + + 8.3751256e-02 5.0000000e-01 6.8727070e-01 9.1624874e-01 3.1272930e-01 + + + + + problem 7 dimension 5 + + + initial l2 norm of the residuals 5.6361303e+11 + + final l2 norm of the residuals 2.7446469e-14 + + number of function evaluations 252 + + number of jacobian evaluations 50 + + exit parameter 1 + + final approximate solution + + 5.0000000e-01 3.1272930e-01 8.3751256e-02 6.8727070e-01 9.1624874e-01 + + + + + problem 7 dimension 6 + + + initial l2 norm of the residuals 2.1547198e-01 + + final l2 norm of the residuals 6.5096426e-13 + + number of function evaluations 14 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 6.6876591e-02 3.6668230e-01 2.8874067e-01 7.1125933e-01 6.3331770e-01 + 9.3312341e-01 + + + + + problem 7 dimension 6 + + + initial l2 norm of the residuals 1.3079247e+08 + + final l2 norm of the residuals 1.1293830e-12 + + number of function evaluations 92 + + number of jacobian evaluations 13 + + exit parameter 1 + + final approximate solution + + 9.3312341e-01 3.6668230e-01 6.6876591e-02 7.1125933e-01 2.8874067e-01 + 6.3331770e-01 + + + + + problem 7 dimension 6 + + + initial l2 norm of the residuals 1.8755789e+14 + + final l2 norm of the residuals 6.1076684e-13 + + number of function evaluations 137 + + number of jacobian evaluations 20 + + exit parameter 1 + + final approximate solution + + 9.3312341e-01 6.3331770e-01 2.8874067e-01 3.6668230e-01 6.6876591e-02 + 7.1125933e-01 + + + + + problem 7 dimension 7 + + + initial l2 norm of the residuals 1.8376789e-01 + + final l2 norm of the residuals 2.2709397e-12 + + number of function evaluations 12 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 5.8069150e-02 2.3517161e-01 3.3804409e-01 5.0000000e-01 6.6195591e-01 + 7.6482839e-01 9.4193085e-01 + + + + + problem 7 dimension 7 + + + initial l2 norm of the residuals 4.2693282e+09 + + final l2 norm of the residuals 6.3507454e-13 + + number of function evaluations 279 + + number of jacobian evaluations 47 + + exit parameter 1 + + final approximate solution + + 2.3517161e-01 7.6482839e-01 5.8069150e-02 6.6195591e-01 9.4193085e-01 + 3.3804409e-01 5.0000000e-01 + + + + + problem 7 dimension 7 + + + initial l2 norm of the residuals 6.4143166e+16 + + final l2 norm of the residuals 6.7782527e+11 + + number of function evaluations 67 + + number of jacobian evaluations 13 + + exit parameter 4 + + final approximate solution + + -2.1047582e+00 -3.1408739e+01 2.9011612e+01 2.3034471e+01 2.8047581e+01 + 2.2965665e+01 1.7978208e+01 + + + + + problem 7 dimension 8 + + + initial l2 norm of the residuals 1.9651386e-01 + + final l2 norm of the residuals 6.4405122e-02 + + number of function evaluations 56 + + number of jacobian evaluations 8 + + exit parameter 4 + + final approximate solution + + 4.9856400e-02 1.9863513e-01 2.6982882e-01 4.9927230e-01 5.0072770e-01 + 7.3017118e-01 8.0136487e-01 9.5014360e-01 + + + + + problem 7 dimension 9 + + + initial l2 norm of the residuals 1.6994993e-01 + + final l2 norm of the residuals 8.3793390e-14 + + number of function evaluations 25 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 4.4205346e-02 1.9949067e-01 2.3561911e-01 4.1604691e-01 5.0000000e-01 + 5.8395309e-01 7.6438089e-01 8.0050933e-01 9.5579465e-01 + + + + + problem 8 dimension 10 + + + initial l2 norm of the residuals 1.6530216e+01 + + final l2 norm of the residuals 1.4636729e-18 + + number of function evaluations 12 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 8 dimension 10 + + + initial l2 norm of the residuals 9.7656240e+06 + + final l2 norm of the residuals 4.0661195e-18 + + number of function evaluations 12 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 8 dimension 10 + + + initial l2 norm of the residuals 9.7656250e+16 + + final l2 norm of the residuals 9.6343123e-18 + + number of function evaluations 15 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 8 dimension 30 + + + initial l2 norm of the residuals 8.3476044e+01 + + final l2 norm of the residuals 1.3604685e-16 + + number of function evaluations 12 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 8 dimension 40 + + + initial l2 norm of the residuals 1.2802636e+02 + + final l2 norm of the residuals 8.9376145e-17 + + number of function evaluations 13 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 9 dimension 10 + + + initial l2 norm of the residuals 2.8080582e-02 + + final l2 norm of the residuals 2.5152585e-15 + + number of function evaluations 6 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 9 dimension 10 + + + initial l2 norm of the residuals 5.2555258e-01 + + final l2 norm of the residuals 1.7259752e-13 + + number of function evaluations 9 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 9 dimension 10 + + + initial l2 norm of the residuals 1.0657390e+02 + + final l2 norm of the residuals 7.1343184e-12 + + number of function evaluations 44 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 10 dimension 1 + + + initial l2 norm of the residuals 1.2792969e-01 + + final l2 norm of the residuals 4.0657581e-20 + + number of function evaluations 7 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -1.5281388e-01 + + + + + problem 10 dimension 1 + + + initial l2 norm of the residuals 2.5625000e+00 + + final l2 norm of the residuals 2.7105054e-20 + + number of function evaluations 8 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -1.5281388e-01 + + + + + problem 10 dimension 1 + + + initial l2 norm of the residuals 8.3611719e+02 + + final l2 norm of the residuals 4.0657581e-20 + + number of function evaluations 25 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -1.5281388e-01 + + + + + problem 10 dimension 10 + + + initial l2 norm of the residuals 2.5182701e-01 + + final l2 norm of the residuals 8.6329069e-15 + + number of function evaluations 9 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 10 dimension 10 + + + initial l2 norm of the residuals 6.1168330e+00 + + final l2 norm of the residuals 4.0860096e-15 + + number of function evaluations 10 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 10 dimension 10 + + + initial l2 norm of the residuals 1.2693089e+03 + + final l2 norm of the residuals 2.7254806e-17 + + number of function evaluations 51 + + number of jacobian evaluations 3 + + exit parameter 1 + + final approximate solution + + -4.3164983e-02 -8.1577157e-02 -1.1448571e-01 -1.4097358e-01 -1.5990870e-01 + -1.6987720e-01 -1.6908998e-01 -1.5524954e-01 -1.2535589e-01 -7.5416534e-02 + + + + + problem 11 dimension 10 + + + initial l2 norm of the residuals 8.4117534e-02 + + final l2 norm of the residuals 5.2963921e-03 + + number of function evaluations 60 + + number of jacobian evaluations 7 + + exit parameter 4 + + final approximate solution + + 5.5261544e-02 5.6957555e-02 5.8890667e-02 6.1136063e-02 6.3778558e-02 + 6.7004825e-02 2.0794103e-01 1.6426711e-01 8.6439479e-02 9.1335063e-02 + + + + + problem 11 dimension 10 + + + initial l2 norm of the residuals 2.0305195e+01 + + final l2 norm of the residuals 3.3454654e-13 + + number of function evaluations 36 + + number of jacobian evaluations 5 + + exit parameter 1 + + final approximate solution + + 3.4396289e-02 3.5032316e-02 3.5719196e-02 3.6465224e-02 3.7280912e-02 + 3.8179863e-02 3.9180141e-02 4.0306503e-02 1.7972019e-01 1.5624088e-01 + + + + + problem 11 dimension 10 + + + initial l2 norm of the residuals 9.3369375e+01 + + final l2 norm of the residuals 1.5143551e-09 + + number of function evaluations 44 + + number of jacobian evaluations 4 + + exit parameter 1 + + final approximate solution + + 1.8883952e+01 2.5167774e+01 1.8885275e+01 1.8886021e+01 1.8886837e+01 + 1.8887736e+01 1.8888736e+01 1.8889862e+01 1.9029276e+01 1.9005797e+01 + + + + + problem 12 dimension 10 + + + initial l2 norm of the residuals 2.2402135e+06 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 22 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 12 dimension 10 + + + initial l2 norm of the residuals 5.2234376e+07 + + final l2 norm of the residuals 1.2183033e-14 + + number of function evaluations 26 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 12 dimension 10 + + + initial l2 norm of the residuals 1.5923646e+11 + + final l2 norm of the residuals 9.9855502e-12 + + number of function evaluations 35 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 13 dimension 10 + + + initial l2 norm of the residuals 4.5825757e+00 + + final l2 norm of the residuals 7.6228700e-11 + + number of function evaluations 13 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -5.7072213e-01 -6.8180695e-01 -7.0221008e-01 -7.0551063e-01 -7.0490616e-01 + -7.0149661e-01 -6.9188932e-01 -6.6579651e-01 -5.9603511e-01 -4.1641226e-01 + + + + + problem 13 dimension 10 + + + initial l2 norm of the residuals 6.3910093e+02 + + final l2 norm of the residuals 2.3156807e-12 + + number of function evaluations 53 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -5.7072213e-01 -6.8180695e-01 -7.0221008e-01 -7.0551063e-01 -7.0490616e-01 + -7.0149661e-01 -6.9188932e-01 -6.6579651e-01 -5.9603511e-01 -4.1641226e-01 + + + + + problem 13 dimension 10 + + + initial l2 norm of the residuals 6.3337583e+04 + + final l2 norm of the residuals 3.2465706e-12 + + number of function evaluations 23 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -5.7072213e-01 -6.8180695e-01 -7.0221008e-01 -7.0551063e-01 -7.0490616e-01 + -7.0149661e-01 -6.9188932e-01 -6.6579651e-01 -5.9603511e-01 -4.1641226e-01 + + + + + problem 14 dimension 10 + + + initial l2 norm of the residuals 1.8973666e+01 + + final l2 norm of the residuals 8.2517941e-11 + + number of function evaluations 23 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + -4.2830286e-01 -4.7659642e-01 -5.1965246e-01 -5.5809932e-01 -5.9250616e-01 + -6.2450368e-01 -6.2323947e-01 -6.2139384e-01 -6.2045360e-01 -5.8646927e-01 + + + + + problem 14 dimension 10 + + + initial l2 norm of the residuals 1.7130922e+04 + + final l2 norm of the residuals 1.3555351e-10 + + number of function evaluations 27 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -4.2830286e-01 -4.7659642e-01 -5.1965246e-01 -5.5809932e-01 -5.9250616e-01 + -6.2450368e-01 -6.2323947e-01 -6.2139384e-01 -6.2045360e-01 -5.8646927e-01 + + + + + problem 14 dimension 10 + + + initial l2 norm of the residuals 1.5949860e+07 + + final l2 norm of the residuals 8.7598527e-11 + + number of function evaluations 39 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -4.2830286e-01 -4.7659642e-01 -5.1965246e-01 -5.5809932e-01 -5.9250616e-01 + -6.2450368e-01 -6.2323947e-01 -6.2139384e-01 -6.2045360e-01 -5.8646927e-01 + summary of 55 calls to hybrj1 + + nprob n nfev njev info final l2 norm + + 1 2 16 3 1 0.0000000e+00 + 1 2 7 1 1 0.0000000e+00 + 1 2 7 1 1 2.6237693e-16 + 2 4 126 5 4 1.0965660e-43 + 2 4 125 4 4 6.0574582e-42 + 2 4 169 9 4 3.6871121e-41 + 3 2 170 6 1 1.1585631e-10 + 3 2 11 1 1 1.0842022e-19 + 4 4 87 2 1 1.0569600e-11 + 4 4 204 8 1 4.3101277e-14 + 4 4 386 35 1 3.4779835e-11 + 5 3 18 3 1 2.7534492e-13 + 5 3 21 4 1 1.4599744e-12 + 5 3 38 8 1 4.7530476e-13 + 6 6 61 6 1 5.7711673e-16 + 6 6 140 6 1 8.2109861e-16 + 6 9 93 4 1 4.7333350e-13 + 6 9 301 13 4 3.3676416e-02 + 7 5 10 1 1 6.0446772e-13 + 7 5 142 25 1 8.7599009e-13 + 7 5 252 50 1 2.7446469e-14 + 7 6 14 2 1 6.5096426e-13 + 7 6 92 13 1 1.1293830e-12 + 7 6 137 20 1 6.1076684e-13 + 7 7 12 1 1 2.2709397e-12 + 7 7 279 47 1 6.3507454e-13 + 7 7 67 13 4 6.7782527e+11 + 7 8 56 8 4 6.4405122e-02 + 7 9 25 2 1 8.3793390e-14 + 8 10 12 2 1 1.4636729e-18 + 8 10 12 2 1 4.0661195e-18 + 8 10 15 2 1 9.6343123e-18 + 8 30 12 2 1 1.3604685e-16 + 8 40 13 2 1 8.9376145e-17 + 9 10 6 1 1 2.5152585e-15 + 9 10 9 1 1 1.7259752e-13 + 9 10 44 1 1 7.1343184e-12 + 10 1 7 1 1 4.0657581e-20 + 10 1 8 1 1 2.7105054e-20 + 10 1 25 2 1 4.0657581e-20 + 10 10 9 1 1 8.6329069e-15 + 10 10 10 1 1 4.0860096e-15 + 10 10 51 3 1 2.7254806e-17 + 11 10 60 7 4 5.2963921e-03 + 11 10 36 5 1 3.3454654e-13 + 11 10 44 4 1 1.5143551e-09 + 12 10 22 1 1 0.0000000e+00 + 12 10 26 1 1 1.2183033e-14 + 12 10 35 1 1 9.9855502e-12 + 13 10 13 1 1 7.6228700e-11 + 13 10 53 1 1 2.3156807e-12 + 13 10 23 2 1 3.2465706e-12 + 14 10 23 1 1 8.2517941e-11 + 14 10 27 2 1 1.3555351e-10 + 14 10 39 2 1 8.7598527e-11 diff --git a/examples/ref/ldibmdpdrc.ref b/examples/ref/ldibmdpdrc.ref new file mode 100644 index 0000000..c2747f4 --- /dev/null +++ b/examples/ref/ldibmdpdrc.ref @@ -0,0 +1,42 @@ + MACHAR constants + + + ibeta = 2 + + it = 64 + + irnd = 5 + + ngrd = 0 + + machep = -63 + + negep = -64 + + iexp = 15 + + minexp =-16382 + + maxexp = 16384 + + eps = 1.0842022e-19 + + epsneg = 5.4210109e-20 + + xmin = 0.0000000e+00 + + xmax = inf + + + + DPMPAR constants and relative differences + + + epsmch = 1.0842022e-19 + rerr(1) = 0.0000000e+00 + + dwarf = 0.0000000e+00 + rerr(2) = 0.0000000e+00 + + giant = inf + rerr(3) = 0.0000000e+00 diff --git a/examples/ref/ldlmddrvc.ref b/examples/ref/ldlmddrvc.ref new file mode 100644 index 0000000..4841eea --- /dev/null +++ b/examples/ref/ldlmddrvc.ref @@ -0,0 +1,1151 @@ + + + + + problem 1 dimensions 5 10 + + + initial l2 norm of the residuals 5.0000000e+00 + + final l2 norm of the residuals 2.2360680e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 3 + + final approximate solution + + -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 + + + + + problem 1 dimensions 5 50 + + + initial l2 norm of the residuals 8.0622577e+00 + + final l2 norm of the residuals 6.7082039e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 3 + + final approximate solution + + -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 + + + + + problem 2 dimensions 5 10 + + + initial l2 norm of the residuals 2.9152187e+02 + + final l2 norm of the residuals 1.4638501e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 2.6727336e+02 1.3413668e+02 -2.0345282e+02 6.7568339e+01 -3.9063749e+01 + + + + + problem 2 dimensions 5 50 + + + initial l2 norm of the residuals 3.1016004e+03 + + final l2 norm of the residuals 3.4826302e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -2.0674605e+01 -9.8373024e+00 1.8372733e+02 -4.4186512e+00 -9.8625693e+01 + + + + + problem 3 dimensions 5 10 + + + initial l2 norm of the residuals 1.2603968e+02 + + final l2 norm of the residuals 1.9097274e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0963164e+02 -1.4678337e+02 5.5315821e+01 1.0000000e+00 + + + + + problem 3 dimensions 5 50 + + + initial l2 norm of the residuals 1.7489500e+03 + + final l2 norm of the residuals 3.6917294e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 -1.8284307e+02 1.3169090e+02 -7.3389069e+00 1.0000000e+00 + + + + + problem 4 dimensions 2 2 + + + initial l2 norm of the residuals 4.9193496e+00 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 21 + + number of jacobian evaluations 16 + + exit parameter 4 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 4 dimensions 2 2 + + + initial l2 norm of the residuals 1.3400631e+03 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 8 + + number of jacobian evaluations 5 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 4 dimensions 2 2 + + + initial l2 norm of the residuals 1.4300005e+05 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 6 + + number of jacobian evaluations 4 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 5 dimensions 3 3 + + + initial l2 norm of the residuals 5.0000000e+01 + + final l2 norm of the residuals 1.4499893e-32 + + number of function evaluations 12 + + number of jacobian evaluations 9 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 -9.1105514e-34 0.0000000e+00 + + + + + problem 5 dimensions 3 3 + + + initial l2 norm of the residuals 1.0295630e+02 + + final l2 norm of the residuals 1.0446754e-19 + + number of function evaluations 20 + + number of jacobian evaluations 15 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 6.5638892e-21 0.0000000e+00 + + + + + problem 5 dimensions 3 3 + + + initial l2 norm of the residuals 9.9126182e+02 + + final l2 norm of the residuals 2.6406807e-29 + + number of function evaluations 19 + + number of jacobian evaluations 16 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 -1.6591886e-30 0.0000000e+00 + + + + + problem 6 dimensions 4 4 + + + initial l2 norm of the residuals 1.4662878e+01 + + final l2 norm of the residuals 1.4565773e-40 + + number of function evaluations 69 + + number of jacobian evaluations 68 + + exit parameter 4 + + final approximate solution + + 8.0669805e-21 -8.0669805e-22 1.2907169e-21 1.2907169e-21 + + + + + problem 6 dimensions 4 4 + + + initial l2 norm of the residuals 1.2709839e+03 + + final l2 norm of the residuals 2.9624937e-169 + + number of function evaluations 500 + + number of jacobian evaluations 499 + + exit parameter 5 + + final approximate solution + + 1.4547316e-149-1.4547316e-150 2.3275706e-150 2.3275706e-150 + + + + + problem 6 dimensions 4 4 + + + initial l2 norm of the residuals 1.2688790e+05 + + final l2 norm of the residuals 8.8902421e-41 + + number of function evaluations 76 + + number of jacobian evaluations 75 + + exit parameter 4 + + final approximate solution + + 6.3023285e-21 -6.3023285e-22 1.0083726e-21 1.0083726e-21 + + + + + problem 7 dimensions 2 2 + + + initial l2 norm of the residuals 2.0012496e+01 + + final l2 norm of the residuals 6.9988752e+00 + + number of function evaluations 18 + + number of jacobian evaluations 11 + + exit parameter 1 + + final approximate solution + + 1.1412779e+01 -8.9680502e-01 + + + + + problem 7 dimensions 2 2 + + + initial l2 norm of the residuals 1.2432834e+04 + + final l2 norm of the residuals 6.9988752e+00 + + number of function evaluations 21 + + number of jacobian evaluations 13 + + exit parameter 1 + + final approximate solution + + 1.1412738e+01 -8.9680828e-01 + + + + + problem 7 dimensions 2 2 + + + initial l2 norm of the residuals 1.1426455e+07 + + final l2 norm of the residuals 6.9988752e+00 + + number of function evaluations 25 + + number of jacobian evaluations 17 + + exit parameter 1 + + final approximate solution + + 1.1412782e+01 -8.9680511e-01 + + + + + problem 8 dimensions 3 15 + + + initial l2 norm of the residuals 6.4561363e+00 + + final l2 norm of the residuals 9.0635960e-02 + + number of function evaluations 7 + + number of jacobian evaluations 6 + + exit parameter 1 + + final approximate solution + + 8.2410560e-02 1.1330361e+00 2.3436952e+00 + + + + + problem 8 dimensions 3 15 + + + initial l2 norm of the residuals 3.6141853e+01 + + final l2 norm of the residuals 4.1747687e+00 + + number of function evaluations 47 + + number of jacobian evaluations 46 + + exit parameter 1 + + final approximate solution + + 8.4066667e-01 -5.0830897e+09 -5.2601486e+09 + + + + + problem 8 dimensions 3 15 + + + initial l2 norm of the residuals 3.8411468e+02 + + final l2 norm of the residuals 4.1747687e+00 + + number of function evaluations 24 + + number of jacobian evaluations 23 + + exit parameter 1 + + final approximate solution + + 8.4066667e-01 -4.9730046e+09 -5.1462126e+09 + + + + + problem 9 dimensions 4 11 + + + initial l2 norm of the residuals 7.2891510e-02 + + final l2 norm of the residuals 1.7535838e-02 + + number of function evaluations 22 + + number of jacobian evaluations 20 + + exit parameter 1 + + final approximate solution + + 1.9280707e-01 1.9127927e-01 1.2305593e-01 1.3606091e-01 + + + + + problem 9 dimensions 4 11 + + + initial l2 norm of the residuals 2.9793701e+00 + + final l2 norm of the residuals 3.2052193e-02 + + number of function evaluations 88 + + number of jacobian evaluations 80 + + exit parameter 1 + + final approximate solution + + 2.3317387e+07 -1.4075883e+01 -1.0552795e+09 -6.5828474e+08 + + + + + problem 9 dimensions 4 11 + + + initial l2 norm of the residuals 2.9959062e+01 + + final l2 norm of the residuals 1.7535838e-02 + + number of function evaluations 499 + + number of jacobian evaluations 385 + + exit parameter 1 + + final approximate solution + + 1.9280705e-01 1.9127981e-01 1.2305603e-01 1.3606116e-01 + + + + + problem 10 dimensions 3 16 + + + initial l2 norm of the residuals 4.1153467e+04 + + final l2 norm of the residuals 9.3779451e+00 + + number of function evaluations 126 + + number of jacobian evaluations 116 + + exit parameter 2 + + final approximate solution + + 5.6096365e-03 6.1813463e+03 3.4522363e+02 + + + + + problem 10 dimensions 3 16 + + + initial l2 norm of the residuals 4.1682169e+06 + + final l2 norm of the residuals 7.9694345e+02 + + number of function evaluations 400 + + number of jacobian evaluations 353 + + exit parameter 5 + + final approximate solution + + 1.1178221e-11 3.4169657e+04 9.0810387e+02 + + + + + problem 11 dimensions 6 31 + + + initial l2 norm of the residuals 5.4772256e+00 + + final l2 norm of the residuals 4.7829594e-02 + + number of function evaluations 9 + + number of jacobian evaluations 8 + + exit parameter 1 + + final approximate solution + + -1.5725102e-02 1.0124349e+00 -2.3299161e-01 1.2604300e+00 -1.5137287e+00 + 9.9299633e-01 + + + + + problem 11 dimensions 6 31 + + + initial l2 norm of the residuals 6.4331258e+03 + + final l2 norm of the residuals 4.7829594e-02 + + number of function evaluations 15 + + number of jacobian evaluations 14 + + exit parameter 1 + + final approximate solution + + -1.5725073e-02 1.0124349e+00 -2.3299164e-01 1.2604302e+00 -1.5137291e+00 + 9.9299652e-01 + + + + + problem 11 dimensions 6 31 + + + initial l2 norm of the residuals 6.7425604e+05 + + final l2 norm of the residuals 4.7829594e-02 + + number of function evaluations 16 + + number of jacobian evaluations 15 + + exit parameter 1 + + final approximate solution + + -1.5725135e-02 1.0124349e+00 -2.3299159e-01 1.2604297e+00 -1.5137284e+00 + 9.9299610e-01 + + + + + problem 11 dimensions 9 31 + + + initial l2 norm of the residuals 5.4772256e+00 + + final l2 norm of the residuals 1.1831146e-03 + + number of function evaluations 8 + + number of jacobian evaluations 7 + + exit parameter 1 + + final approximate solution + + -1.5307037e-05 9.9978970e-01 1.4763964e-02 1.4634233e-01 1.0008211e+00 + -2.6177311e+00 4.1044032e+00 -3.1436123e+00 1.0526264e+00 + + + + + problem 11 dimensions 9 31 + + + initial l2 norm of the residuals 1.2088127e+04 + + final l2 norm of the residuals 1.1831146e-03 + + number of function evaluations 19 + + number of jacobian evaluations 15 + + exit parameter 1 + + final approximate solution + + -1.5307036e-05 9.9978970e-01 1.4763964e-02 1.4634233e-01 1.0008211e+00 + -2.6177311e+00 4.1044032e+00 -3.1436123e+00 1.0526264e+00 + + + + + problem 11 dimensions 9 31 + + + initial l2 norm of the residuals 1.2691093e+06 + + final l2 norm of the residuals 1.1831146e-03 + + number of function evaluations 18 + + number of jacobian evaluations 15 + + exit parameter 1 + + final approximate solution + + -1.5307037e-05 9.9978970e-01 1.4763964e-02 1.4634233e-01 1.0008211e+00 + -2.6177311e+00 4.1044032e+00 -3.1436123e+00 1.0526264e+00 + + + + + problem 11 dimensions 12 31 + + + initial l2 norm of the residuals 5.4772256e+00 + + final l2 norm of the residuals 2.1731040e-05 + + number of function evaluations 9 + + number of jacobian evaluations 8 + + exit parameter 1 + + final approximate solution + + -6.6380605e-09 1.0000016e+00 -5.6393221e-04 3.4782054e-01 -1.5673150e-01 + 1.0528152e+00 -3.2472712e+00 7.2884349e+00 -1.0271848e+01 9.0741136e+00 + -4.5413755e+00 1.0120119e+00 + + + + + problem 11 dimensions 12 31 + + + initial l2 norm of the residuals 1.9220759e+04 + + final l2 norm of the residuals 2.1731040e-05 + + number of function evaluations 14 + + number of jacobian evaluations 13 + + exit parameter 3 + + final approximate solution + + -6.6380605e-09 1.0000016e+00 -5.6393221e-04 3.4782054e-01 -1.5673150e-01 + 1.0528152e+00 -3.2472712e+00 7.2884349e+00 -1.0271848e+01 9.0741136e+00 + -4.5413755e+00 1.0120119e+00 + + + + + problem 11 dimensions 12 31 + + + initial l2 norm of the residuals 2.0189180e+06 + + final l2 norm of the residuals 2.1731040e-05 + + number of function evaluations 34 + + number of jacobian evaluations 28 + + exit parameter 2 + + final approximate solution + + -6.6380605e-09 1.0000016e+00 -5.6393221e-04 3.4782054e-01 -1.5673150e-01 + 1.0528152e+00 -3.2472712e+00 7.2884349e+00 -1.0271848e+01 9.0741136e+00 + -4.5413755e+00 1.0120119e+00 + + + + + problem 12 dimensions 3 10 + + + initial l2 norm of the residuals 3.2111584e+01 + + final l2 norm of the residuals 1.7167030e-16 + + number of function evaluations 8 + + number of jacobian evaluations 7 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+01 1.0000000e+00 + + + + + problem 13 dimensions 2 10 + + + initial l2 norm of the residuals 6.4585650e+01 + + final l2 norm of the residuals 1.1151779e+01 + + number of function evaluations 24 + + number of jacobian evaluations 13 + + exit parameter 1 + + final approximate solution + + 2.5782487e-01 2.5782554e-01 + + + + + problem 14 dimensions 4 20 + + + initial l2 norm of the residuals 2.8154384e+03 + + final l2 norm of the residuals 2.9295427e+02 + + number of function evaluations 300 + + number of jacobian evaluations 278 + + exit parameter 1 + + final approximate solution + + -1.1594164e+01 1.3203531e+01 -4.0345234e-01 2.3677411e-01 + + + + + problem 14 dimensions 4 20 + + + initial l2 norm of the residuals 5.5507335e+05 + + final l2 norm of the residuals 2.9295427e+02 + + number of function evaluations 68 + + number of jacobian evaluations 54 + + exit parameter 1 + + final approximate solution + + -1.1594593e+01 1.3203687e+01 -4.0341581e-01 2.3679328e-01 + + + + + problem 14 dimensions 4 20 + + + initial l2 norm of the residuals 6.1211252e+07 + + final l2 norm of the residuals 2.9295427e+02 + + number of function evaluations 310 + + number of jacobian evaluations 289 + + exit parameter 1 + + final approximate solution + + -1.1593993e+01 1.3203462e+01 -4.0345948e-01 2.3677112e-01 + + + + + problem 15 dimensions 1 8 + + + initial l2 norm of the residuals 1.8862380e+00 + + final l2 norm of the residuals 1.8862380e+00 + + number of function evaluations 1 + + number of jacobian evaluations 1 + + exit parameter 4 + + final approximate solution + + 5.0000000e-01 + + + + + problem 15 dimensions 1 8 + + + initial l2 norm of the residuals 5.3833444e+09 + + final l2 norm of the residuals 1.8842482e+00 + + number of function evaluations 30 + + number of jacobian evaluations 29 + + exit parameter 1 + + final approximate solution + + 9.8173132e-01 + + + + + problem 15 dimensions 1 8 + + + initial l2 norm of the residuals 1.1808873e+18 + + final l2 norm of the residuals 1.8842482e+00 + + number of function evaluations 48 + + number of jacobian evaluations 47 + + exit parameter 1 + + final approximate solution + + 9.8173132e-01 + + + + + problem 15 dimensions 8 8 + + + initial l2 norm of the residuals 1.9651386e-01 + + final l2 norm of the residuals 5.9303235e-02 + + number of function evaluations 45 + + number of jacobian evaluations 23 + + exit parameter 1 + + final approximate solution + + 4.3152867e-02 1.9309093e-01 2.6632869e-01 5.0000008e-01 4.9999992e-01 + 7.3367131e-01 8.0690907e-01 9.5684713e-01 + + + + + problem 15 dimensions 9 9 + + + initial l2 norm of the residuals 1.6994993e-01 + + final l2 norm of the residuals 3.3033507e-19 + + number of function evaluations 12 + + number of jacobian evaluations 9 + + exit parameter 2 + + final approximate solution + + 4.4205346e-02 1.9949067e-01 2.3561911e-01 4.1604691e-01 5.0000000e-01 + 5.8395309e-01 7.6438089e-01 8.0050933e-01 9.5579465e-01 + + + + + problem 15 dimensions 10 10 + + + initial l2 norm of the residuals 1.8374783e-01 + + final l2 norm of the residuals 8.0647100e-02 + + number of function evaluations 29 + + number of jacobian evaluations 14 + + exit parameter 1 + + final approximate solution + + 5.9619952e-02 1.6670835e-01 2.3917071e-01 3.9888443e-01 3.9888421e-01 + 6.0111579e-01 6.0111557e-01 7.6082929e-01 8.3329165e-01 9.4038005e-01 + + + + + problem 16 dimensions 10 10 + + + initial l2 norm of the residuals 1.6530216e+01 + + final l2 norm of the residuals 6.8570966e-19 + + number of function evaluations 15 + + number of jacobian evaluations 13 + + exit parameter 2 + + final approximate solution + + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 1.2056970e+00 + + + + + problem 16 dimensions 10 10 + + + initial l2 norm of the residuals 9.7656240e+06 + + final l2 norm of the residuals 7.2730493e-19 + + number of function evaluations 13 + + number of jacobian evaluations 8 + + exit parameter 2 + + final approximate solution + + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 1.2056970e+00 + + + + + problem 16 dimensions 10 10 + + + initial l2 norm of the residuals 9.7656250e+16 + + final l2 norm of the residuals 2.9293530e-18 + + number of function evaluations 30 + + number of jacobian evaluations 27 + + exit parameter 2 + + final approximate solution + + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 1.2056970e+00 + + + + + problem 16 dimensions 30 30 + + + initial l2 norm of the residuals 8.3476044e+01 + + final l2 norm of the residuals 3.8012846e-18 + + number of function evaluations 19 + + number of jacobian evaluations 14 + + exit parameter 2 + + final approximate solution + + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 1.0673735e+00 + + + + + problem 16 dimensions 40 40 + + + initial l2 norm of the residuals 1.2802636e+02 + + final l2 norm of the residuals 1.1855046e-17 + + number of function evaluations 19 + + number of jacobian evaluations 14 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 17 dimensions 5 33 + + + initial l2 norm of the residuals 9.3756402e-01 + + final l2 norm of the residuals 7.3924926e-03 + + number of function evaluations 19 + + number of jacobian evaluations 16 + + exit parameter 1 + + final approximate solution + + 3.7541005e-01 1.9358469e+00 -1.4646871e+00 1.2867535e-02 2.2122700e-02 + + + + + problem 18 dimensions 11 65 + + + initial l2 norm of the residuals 1.4468654e+00 + + final l2 norm of the residuals 2.0034404e-01 + + number of function evaluations 17 + + number of jacobian evaluations 13 + + exit parameter 1 + + final approximate solution + + 1.3099771e+00 4.3155362e-01 6.3366163e-01 5.9943049e-01 7.5418274e-01 + 9.0428996e-01 1.3658118e+00 4.8236970e+00 2.3986850e+00 4.5688745e+00 + 5.6753415e+00 + summary of 53 calls to lmder1: + + + + nprob n m nfev njev info final L2 norm + + 1 5 10 3 2 3 2.2360680e+00 + 1 5 50 3 2 3 6.7082039e+00 + 2 5 10 3 2 1 1.4638501e+00 + 2 5 50 3 2 1 3.4826302e+00 + 3 5 10 3 2 1 1.9097274e+00 + 3 5 50 3 2 1 3.6917294e+00 + 4 2 2 21 16 4 0.0000000e+00 + 4 2 2 8 5 2 0.0000000e+00 + 4 2 2 6 4 2 0.0000000e+00 + 5 3 3 12 9 2 1.4499893e-32 + 5 3 3 20 15 2 1.0446754e-19 + 5 3 3 19 16 2 2.6406807e-29 + 6 4 4 69 68 4 1.4565773e-40 + 6 4 4 500 499 5 2.9624937e-169 + 6 4 4 76 75 4 8.8902421e-41 + 7 2 2 18 11 1 6.9988752e+00 + 7 2 2 21 13 1 6.9988752e+00 + 7 2 2 25 17 1 6.9988752e+00 + 8 3 15 7 6 1 9.0635960e-02 + 8 3 15 47 46 1 4.1747687e+00 + 8 3 15 24 23 1 4.1747687e+00 + 9 4 11 22 20 1 1.7535838e-02 + 9 4 11 88 80 1 3.2052193e-02 + 9 4 11 499 385 1 1.7535838e-02 + 10 3 16 126 116 2 9.3779451e+00 + 10 3 16 400 353 5 7.9694345e+02 + 11 6 31 9 8 1 4.7829594e-02 + 11 6 31 15 14 1 4.7829594e-02 + 11 6 31 16 15 1 4.7829594e-02 + 11 9 31 8 7 1 1.1831146e-03 + 11 9 31 19 15 1 1.1831146e-03 + 11 9 31 18 15 1 1.1831146e-03 + 11 12 31 9 8 1 2.1731040e-05 + 11 12 31 14 13 3 2.1731040e-05 + 11 12 31 34 28 2 2.1731040e-05 + 12 3 10 8 7 2 1.7167030e-16 + 13 2 10 24 13 1 1.1151779e+01 + 14 4 20 300 278 1 2.9295427e+02 + 14 4 20 68 54 1 2.9295427e+02 + 14 4 20 310 289 1 2.9295427e+02 + 15 1 8 1 1 4 1.8862380e+00 + 15 1 8 30 29 1 1.8842482e+00 + 15 1 8 48 47 1 1.8842482e+00 + 15 8 8 45 23 1 5.9303235e-02 + 15 9 9 12 9 2 3.3033507e-19 + 15 10 10 29 14 1 8.0647100e-02 + 16 10 10 15 13 2 6.8570966e-19 + 16 10 10 13 8 2 7.2730493e-19 + 16 10 10 30 27 2 2.9293530e-18 + 16 30 30 19 14 2 3.8012846e-18 + 16 40 40 19 14 2 1.1855046e-17 + 17 5 33 19 16 1 7.3924926e-03 + 18 11 65 17 13 1 2.0034404e-01 diff --git a/examples/ref/ldlmfdrvc.ref b/examples/ref/ldlmfdrvc.ref new file mode 100644 index 0000000..1a4e504 --- /dev/null +++ b/examples/ref/ldlmfdrvc.ref @@ -0,0 +1,1151 @@ + + + + + problem 1 dimensions 5 10 + + + initial l2 norm of the residuals 5.0000000e+00 + + final l2 norm of the residuals 2.2360680e+00 + + number of function evaluations 3 + + number of Jacobian evaluations 2 + + exit parameter 3 + + final approximate solution + + -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 + + + + + problem 1 dimensions 5 50 + + + initial l2 norm of the residuals 8.0622577e+00 + + final l2 norm of the residuals 6.7082039e+00 + + number of function evaluations 3 + + number of Jacobian evaluations 2 + + exit parameter 3 + + final approximate solution + + -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 + + + + + problem 2 dimensions 5 10 + + + initial l2 norm of the residuals 2.9152187e+02 + + final l2 norm of the residuals 1.4638501e+00 + + number of function evaluations 3 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -3.7789553e+02 1.4736529e+02 -3.1027917e+02 1.3964860e+02 9.1110183e+01 + + + + + problem 2 dimensions 5 50 + + + initial l2 norm of the residuals 3.1016004e+03 + + final l2 norm of the residuals 3.4826302e+00 + + number of function evaluations 5 + + number of Jacobian evaluations 3 + + exit parameter 1 + + final approximate solution + + 1.1735027e+02 2.7070435e+02 -8.4272019e+01 -1.0001270e+02 -1.1724834e+00 + + + + + problem 3 dimensions 5 10 + + + initial l2 norm of the residuals 1.2603968e+02 + + final l2 norm of the residuals 1.9097274e+00 + + number of function evaluations 3 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0946948e+02 7.3313314e+01 -1.0967561e+02 1.0000000e+00 + + + + + problem 3 dimensions 5 50 + + + initial l2 norm of the residuals 1.7489500e+03 + + final l2 norm of the residuals 3.6917294e+00 + + number of function evaluations 3 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 3.0278467e+02 9.7301759e+01 -2.2436092e+02 1.0000000e+00 + + + + + problem 4 dimensions 2 2 + + + initial l2 norm of the residuals 4.9193496e+00 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 22 + + number of Jacobian evaluations 16 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 4 dimensions 2 2 + + + initial l2 norm of the residuals 1.3400631e+03 + + final l2 norm of the residuals 5.4210109e-19 + + number of function evaluations 8 + + number of Jacobian evaluations 5 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 4 dimensions 2 2 + + + initial l2 norm of the residuals 1.4300005e+05 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 7 + + number of Jacobian evaluations 5 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 5 dimensions 3 3 + + + initial l2 norm of the residuals 5.0000000e+01 + + final l2 norm of the residuals 1.9983203e-21 + + number of function evaluations 12 + + number of Jacobian evaluations 9 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 -1.2555817e-22 -1.0779180e-61 + + + + + problem 5 dimensions 3 3 + + + initial l2 norm of the residuals 1.0295630e+02 + + final l2 norm of the residuals 3.6115952e-16 + + number of function evaluations 20 + + number of Jacobian evaluations 15 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 -2.2692322e-17 6.2685579e-43 + + + + + problem 5 dimensions 3 3 + + + initial l2 norm of the residuals 9.9126182e+02 + + final l2 norm of the residuals 2.6296676e-20 + + number of function evaluations 19 + + number of Jacobian evaluations 16 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.6522689e-21 7.1313297e-94 + + + + + problem 6 dimensions 4 4 + + + initial l2 norm of the residuals 1.4662878e+01 + + final l2 norm of the residuals 2.2046973e-24 + + number of function evaluations 238 + + number of Jacobian evaluations 191 + + exit parameter 5 + + final approximate solution + + 1.2564350e-12 -1.2564350e-13 5.8689177e-13 5.8689177e-13 + + + + + problem 6 dimensions 4 4 + + + initial l2 norm of the residuals 1.2709839e+03 + + final l2 norm of the residuals 3.6276298e-26 + + number of function evaluations 233 + + number of Jacobian evaluations 192 + + exit parameter 5 + + final approximate solution + + 1.6092406e-13 -1.6092406e-14 7.5810191e-14 7.5810191e-14 + + + + + problem 6 dimensions 4 4 + + + initial l2 norm of the residuals 1.2688790e+05 + + final l2 norm of the residuals 8.1822173e-30 + + number of function evaluations 98 + + number of Jacobian evaluations 79 + + exit parameter 2 + + final approximate solution + + 2.4171418e-15 -2.4171418e-16 1.1378998e-15 1.1378998e-15 + + + + + problem 7 dimensions 2 2 + + + initial l2 norm of the residuals 2.0012496e+01 + + final l2 norm of the residuals 6.9988752e+00 + + number of function evaluations 18 + + number of Jacobian evaluations 11 + + exit parameter 1 + + final approximate solution + + 1.1412779e+01 -8.9680502e-01 + + + + + problem 7 dimensions 2 2 + + + initial l2 norm of the residuals 1.2432834e+04 + + final l2 norm of the residuals 6.9988752e+00 + + number of function evaluations 21 + + number of Jacobian evaluations 13 + + exit parameter 1 + + final approximate solution + + 1.1412738e+01 -8.9680828e-01 + + + + + problem 7 dimensions 2 2 + + + initial l2 norm of the residuals 1.1426455e+07 + + final l2 norm of the residuals 6.9988752e+00 + + number of function evaluations 25 + + number of Jacobian evaluations 17 + + exit parameter 1 + + final approximate solution + + 1.1412782e+01 -8.9680511e-01 + + + + + problem 8 dimensions 3 15 + + + initial l2 norm of the residuals 6.4561363e+00 + + final l2 norm of the residuals 9.0635960e-02 + + number of function evaluations 7 + + number of Jacobian evaluations 6 + + exit parameter 1 + + final approximate solution + + 8.2410560e-02 1.1330361e+00 2.3436952e+00 + + + + + problem 8 dimensions 3 15 + + + initial l2 norm of the residuals 3.6141853e+01 + + final l2 norm of the residuals 4.1747687e+00 + + number of function evaluations 50 + + number of Jacobian evaluations 49 + + exit parameter 1 + + final approximate solution + + 8.4066667e-01 -6.1162285e+09 -6.2475746e+09 + + + + + problem 8 dimensions 3 15 + + + initial l2 norm of the residuals 3.8411468e+02 + + final l2 norm of the residuals 4.1747687e+00 + + number of function evaluations 25 + + number of Jacobian evaluations 24 + + exit parameter 1 + + final approximate solution + + 8.4066667e-01 -3.7863397e+09 -3.8816500e+09 + + + + + problem 9 dimensions 4 11 + + + initial l2 norm of the residuals 7.2891510e-02 + + final l2 norm of the residuals 1.7535838e-02 + + number of function evaluations 22 + + number of Jacobian evaluations 20 + + exit parameter 1 + + final approximate solution + + 1.9280707e-01 1.9127927e-01 1.2305593e-01 1.3606091e-01 + + + + + problem 9 dimensions 4 11 + + + initial l2 norm of the residuals 2.9793701e+00 + + final l2 norm of the residuals 3.2052193e-02 + + number of function evaluations 93 + + number of Jacobian evaluations 85 + + exit parameter 1 + + final approximate solution + + 1.5871311e+07 -1.4075883e+01 -7.1829097e+08 -4.4807087e+08 + + + + + problem 9 dimensions 4 11 + + + initial l2 norm of the residuals 2.9959062e+01 + + final l2 norm of the residuals 3.0718476e-02 + + number of function evaluations 252 + + number of Jacobian evaluations 187 + + exit parameter 5 + + final approximate solution + + 2.3491837e-03 1.0196061e+03 7.3473117e+00 4.9602506e+00 + + + + + problem 10 dimensions 3 16 + + + initial l2 norm of the residuals 4.1153467e+04 + + final l2 norm of the residuals 9.3779452e+00 + + number of function evaluations 129 + + number of Jacobian evaluations 116 + + exit parameter 3 + + final approximate solution + + 5.6097045e-03 6.1813362e+03 3.4522329e+02 + + + + + problem 10 dimensions 3 16 + + + initial l2 norm of the residuals 4.1682169e+06 + + final l2 norm of the residuals 6.2375992e+04 + + number of function evaluations 4 + + number of Jacobian evaluations 3 + + exit parameter 4 + + final approximate solution + + 1.3765544e+00 -3.7421368e+03 -3.2123769e+01 + + + + + problem 11 dimensions 6 31 + + + initial l2 norm of the residuals 5.4772256e+00 + + final l2 norm of the residuals 5.1046620e-02 + + number of function evaluations 8 + + number of Jacobian evaluations 7 + + exit parameter 1 + + final approximate solution + + -1.5275022e-29 1.0136384e+00 -2.4430325e-01 1.3737704e+00 -1.6856389e+00 + 1.0980968e+00 + + + + + problem 11 dimensions 6 31 + + + initial l2 norm of the residuals 6.4331258e+03 + + final l2 norm of the residuals 4.7829594e-02 + + number of function evaluations 15 + + number of Jacobian evaluations 14 + + exit parameter 1 + + final approximate solution + + -1.5725072e-02 1.0124349e+00 -2.3299164e-01 1.2604302e+00 -1.5137291e+00 + 9.9299653e-01 + + + + + problem 11 dimensions 6 31 + + + initial l2 norm of the residuals 6.7425604e+05 + + final l2 norm of the residuals 4.7829594e-02 + + number of function evaluations 16 + + number of Jacobian evaluations 15 + + exit parameter 1 + + final approximate solution + + -1.5725136e-02 1.0124349e+00 -2.3299159e-01 1.2604297e+00 -1.5137284e+00 + 9.9299610e-01 + + + + + problem 11 dimensions 9 31 + + + initial l2 norm of the residuals 5.4772256e+00 + + final l2 norm of the residuals 1.1832136e-03 + + number of function evaluations 7 + + number of Jacobian evaluations 6 + + exit parameter 1 + + final approximate solution + + -2.7023764e-28 9.9978967e-01 1.4781655e-02 1.4631229e-01 1.0009910e+00 + -2.6181665e+00 4.1050509e+00 -3.1441034e+00 1.0527847e+00 + + + + + problem 11 dimensions 9 31 + + + initial l2 norm of the residuals 1.2088127e+04 + + final l2 norm of the residuals 1.1831146e-03 + + number of function evaluations 19 + + number of Jacobian evaluations 15 + + exit parameter 1 + + final approximate solution + + -1.5314467e-05 9.9978970e-01 1.4763968e-02 1.4634221e-01 1.0008216e+00 + -2.6177322e+00 4.1044044e+00 -3.1436130e+00 1.0526266e+00 + + + + + problem 11 dimensions 9 31 + + + initial l2 norm of the residuals 1.2691093e+06 + + final l2 norm of the residuals 1.1831146e-03 + + number of function evaluations 19 + + number of Jacobian evaluations 16 + + exit parameter 1 + + final approximate solution + + -1.5282761e-05 9.9978970e-01 1.4764045e-02 1.4634184e-01 1.0008231e+00 + -2.6177354e+00 4.1044082e+00 -3.1436154e+00 1.0526272e+00 + + + + + problem 11 dimensions 12 31 + + + initial l2 norm of the residuals 5.4772256e+00 + + final l2 norm of the residuals 2.1731041e-05 + + number of function evaluations 19 + + number of Jacobian evaluations 9 + + exit parameter 2 + + final approximate solution + + -7.3656773e-27 1.0000016e+00 -5.6389147e-04 3.4782034e-01 -1.5673230e-01 + 1.0528279e+00 -3.2473286e+00 7.2885719e+00 -1.0272040e+01 9.0742728e+00 + -4.5414477e+00 1.0120258e+00 + + + + + problem 11 dimensions 12 31 + + + initial l2 norm of the residuals 1.9220759e+04 + + final l2 norm of the residuals 2.1731044e-05 + + number of function evaluations 30 + + number of Jacobian evaluations 13 + + exit parameter 2 + + final approximate solution + + -4.1716399e-09 1.0000016e+00 -5.6442758e-04 3.4783002e-01 -1.5681262e-01 + 1.0532021e+00 -3.2483976e+00 7.2905183e+00 -1.0274308e+01 9.0759099e+00 + -4.5421149e+00 1.0121431e+00 + + + + + problem 11 dimensions 12 31 + + + initial l2 norm of the residuals 2.0189180e+06 + + final l2 norm of the residuals 2.1731042e-05 + + number of function evaluations 50 + + number of Jacobian evaluations 29 + + exit parameter 2 + + final approximate solution + + -9.1636503e-09 1.0000016e+00 -5.6360827e-04 3.4781495e-01 -1.5668678e-01 + 1.0526119e+00 -3.2466999e+00 7.2874069e+00 -1.0270660e+01 9.0732598e+00 + -4.5410283e+00 1.0119508e+00 + + + + + problem 12 dimensions 3 10 + + + initial l2 norm of the residuals 3.2111584e+01 + + final l2 norm of the residuals 1.4671279e-16 + + number of function evaluations 8 + + number of Jacobian evaluations 7 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+01 1.0000000e+00 + + + + + problem 13 dimensions 2 10 + + + initial l2 norm of the residuals 6.4585650e+01 + + final l2 norm of the residuals 1.1151779e+01 + + number of function evaluations 24 + + number of Jacobian evaluations 13 + + exit parameter 1 + + final approximate solution + + 2.5782487e-01 2.5782554e-01 + + + + + problem 14 dimensions 4 20 + + + initial l2 norm of the residuals 2.8154384e+03 + + final l2 norm of the residuals 2.9295453e+02 + + number of function evaluations 216 + + number of Jacobian evaluations 197 + + exit parameter 5 + + final approximate solution + + -1.1583810e+01 1.3199835e+01 -4.0532902e-01 2.3530088e-01 + + + + + problem 14 dimensions 4 20 + + + initial l2 norm of the residuals 5.5507335e+05 + + final l2 norm of the residuals 2.9295427e+02 + + number of function evaluations 67 + + number of Jacobian evaluations 53 + + exit parameter 1 + + final approximate solution + + -1.1594602e+01 1.3203691e+01 -4.0346539e-01 2.3675676e-01 + + + + + problem 14 dimensions 4 20 + + + initial l2 norm of the residuals 6.1211252e+07 + + final l2 norm of the residuals 2.9295444e+02 + + number of function evaluations 212 + + number of Jacobian evaluations 197 + + exit parameter 5 + + final approximate solution + + -1.1585880e+01 1.3200427e+01 -4.0504695e-01 2.3579221e-01 + + + + + problem 15 dimensions 1 8 + + + initial l2 norm of the residuals 1.8862380e+00 + + final l2 norm of the residuals 1.8862380e+00 + + number of function evaluations 2 + + number of Jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + 5.0000000e-01 + + + + + problem 15 dimensions 1 8 + + + initial l2 norm of the residuals 5.3833444e+09 + + final l2 norm of the residuals 1.8842482e+00 + + number of function evaluations 30 + + number of Jacobian evaluations 29 + + exit parameter 1 + + final approximate solution + + 9.8173132e-01 + + + + + problem 15 dimensions 1 8 + + + initial l2 norm of the residuals 1.1808873e+18 + + final l2 norm of the residuals 1.8842482e+00 + + number of function evaluations 48 + + number of Jacobian evaluations 47 + + exit parameter 1 + + final approximate solution + + 9.8173132e-01 + + + + + problem 15 dimensions 8 8 + + + initial l2 norm of the residuals 1.9651386e-01 + + final l2 norm of the residuals 5.9303235e-02 + + number of function evaluations 45 + + number of Jacobian evaluations 23 + + exit parameter 1 + + final approximate solution + + 4.3152867e-02 1.9309093e-01 2.6632869e-01 5.0000008e-01 4.9999992e-01 + 7.3367131e-01 8.0690907e-01 9.5684713e-01 + + + + + problem 15 dimensions 9 9 + + + initial l2 norm of the residuals 1.6994993e-01 + + final l2 norm of the residuals 3.9694447e-19 + + number of function evaluations 12 + + number of Jacobian evaluations 9 + + exit parameter 2 + + final approximate solution + + 4.4205346e-02 1.9949067e-01 2.3561911e-01 4.1604691e-01 5.0000000e-01 + 5.8395309e-01 7.6438089e-01 8.0050933e-01 9.5579465e-01 + + + + + problem 15 dimensions 10 10 + + + initial l2 norm of the residuals 1.8374783e-01 + + final l2 norm of the residuals 8.0647100e-02 + + number of function evaluations 29 + + number of Jacobian evaluations 14 + + exit parameter 1 + + final approximate solution + + 5.9619952e-02 1.6670835e-01 2.3917071e-01 3.9888443e-01 3.9888421e-01 + 6.0111579e-01 6.0111557e-01 7.6082929e-01 8.3329165e-01 9.4038005e-01 + + + + + problem 16 dimensions 10 10 + + + initial l2 norm of the residuals 1.6530216e+01 + + final l2 norm of the residuals 6.1801901e-18 + + number of function evaluations 15 + + number of Jacobian evaluations 13 + + exit parameter 2 + + final approximate solution + + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 1.2056970e+00 + + + + + problem 16 dimensions 10 10 + + + initial l2 norm of the residuals 9.7656240e+06 + + final l2 norm of the residuals 4.3923571e-18 + + number of function evaluations 13 + + number of Jacobian evaluations 8 + + exit parameter 2 + + final approximate solution + + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 1.2056970e+00 + + + + + problem 16 dimensions 10 10 + + + initial l2 norm of the residuals 9.7656250e+16 + + final l2 norm of the residuals 1.0000000e+00 + + number of function evaluations 22 + + number of Jacobian evaluations 20 + + exit parameter 3 + + final approximate solution + + 2.3769773e-02 2.3769773e-02 2.3769773e-02 2.3769773e-02 2.3769773e-02 + 2.3769773e-02 2.3769773e-02 2.3769773e-02 2.3769773e-02 1.0762302e+01 + + + + + problem 16 dimensions 30 30 + + + initial l2 norm of the residuals 8.3476044e+01 + + final l2 norm of the residuals 1.0000000e+00 + + number of function evaluations 4 + + number of Jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 5.8080005e-02 5.8080005e-02 5.8080005e-02 5.8080005e-02 5.8080005e-02 + 5.8080005e-02 5.8080005e-02 5.8080005e-02 5.8080005e-02 5.8080005e-02 + 5.8080005e-02 5.8080005e-02 5.8080005e-02 5.8080005e-02 5.8080005e-02 + 5.8080005e-02 5.8080005e-02 5.8080005e-02 5.8080005e-02 5.8080005e-02 + 5.8080005e-02 5.8080005e-02 5.8080005e-02 5.8080005e-02 5.8080005e-02 + 5.8080005e-02 5.8080005e-02 5.8080005e-02 5.8080005e-02 2.9257600e+01 + + + + + problem 16 dimensions 40 40 + + + initial l2 norm of the residuals 1.2802636e+02 + + final l2 norm of the residuals 1.0258455e-16 + + number of function evaluations 11 + + number of Jacobian evaluations 10 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 17 dimensions 5 33 + + + initial l2 norm of the residuals 9.3756402e-01 + + final l2 norm of the residuals 7.3924926e-03 + + number of function evaluations 19 + + number of Jacobian evaluations 16 + + exit parameter 1 + + final approximate solution + + 3.7541004e-01 1.9358454e+00 -1.4646856e+00 1.2867532e-02 2.2122706e-02 + + + + + problem 18 dimensions 11 65 + + + initial l2 norm of the residuals 1.4468654e+00 + + final l2 norm of the residuals 2.0034404e-01 + + number of function evaluations 17 + + number of Jacobian evaluations 13 + + exit parameter 1 + + final approximate solution + + 1.3099771e+00 4.3155361e-01 6.3366163e-01 5.9943051e-01 7.5418270e-01 + 9.0429004e-01 1.3658119e+00 4.8236967e+00 2.3986850e+00 4.5688744e+00 + 5.6753414e+00 + summary of 53 calls to lmdif1: + + + + nprob n m nfev njev info final L2 norm + + 1 5 10 3 2 3 2.2360680e+00 + 1 5 50 3 2 3 6.7082039e+00 + 2 5 10 3 2 1 1.4638501e+00 + 2 5 50 5 3 1 3.4826302e+00 + 3 5 10 3 2 1 1.9097274e+00 + 3 5 50 3 2 1 3.6917294e+00 + 4 2 2 22 16 2 0.0000000e+00 + 4 2 2 8 5 2 5.4210109e-19 + 4 2 2 7 5 2 0.0000000e+00 + 5 3 3 12 9 2 1.9983203e-21 + 5 3 3 20 15 2 3.6115952e-16 + 5 3 3 19 16 2 2.6296676e-20 + 6 4 4 238 191 5 2.2046973e-24 + 6 4 4 233 192 5 3.6276298e-26 + 6 4 4 98 79 2 8.1822173e-30 + 7 2 2 18 11 1 6.9988752e+00 + 7 2 2 21 13 1 6.9988752e+00 + 7 2 2 25 17 1 6.9988752e+00 + 8 3 15 7 6 1 9.0635960e-02 + 8 3 15 50 49 1 4.1747687e+00 + 8 3 15 25 24 1 4.1747687e+00 + 9 4 11 22 20 1 1.7535838e-02 + 9 4 11 93 85 1 3.2052193e-02 + 9 4 11 252 187 5 3.0718476e-02 + 10 3 16 129 116 3 9.3779452e+00 + 10 3 16 4 3 4 6.2375992e+04 + 11 6 31 8 7 1 5.1046620e-02 + 11 6 31 15 14 1 4.7829594e-02 + 11 6 31 16 15 1 4.7829594e-02 + 11 9 31 7 6 1 1.1832136e-03 + 11 9 31 19 15 1 1.1831146e-03 + 11 9 31 19 16 1 1.1831146e-03 + 11 12 31 19 9 2 2.1731041e-05 + 11 12 31 30 13 2 2.1731044e-05 + 11 12 31 50 29 2 2.1731042e-05 + 12 3 10 8 7 2 1.4671279e-16 + 13 2 10 24 13 1 1.1151779e+01 + 14 4 20 216 197 5 2.9295453e+02 + 14 4 20 67 53 1 2.9295427e+02 + 14 4 20 212 197 5 2.9295444e+02 + 15 1 8 2 1 1 1.8862380e+00 + 15 1 8 30 29 1 1.8842482e+00 + 15 1 8 48 47 1 1.8842482e+00 + 15 8 8 45 23 1 5.9303235e-02 + 15 9 9 12 9 2 3.9694447e-19 + 15 10 10 29 14 1 8.0647100e-02 + 16 10 10 15 13 2 6.1801901e-18 + 16 10 10 13 8 2 4.3923571e-18 + 16 10 10 22 20 3 1.0000000e+00 + 16 30 30 4 2 1 1.0000000e+00 + 16 40 40 11 10 2 1.0258455e-16 + 17 5 33 19 16 1 7.3924926e-03 + 18 11 65 17 13 1 2.0034404e-01 diff --git a/examples/ref/ldlmsdrvc.ref b/examples/ref/ldlmsdrvc.ref new file mode 100644 index 0000000..e7e640c --- /dev/null +++ b/examples/ref/ldlmsdrvc.ref @@ -0,0 +1,1148 @@ + + + + + problem 1 dimensions 5 10 + + + initial l2 norm of the residuals 5.0000000e+00 + + final l2 norm of the residuals 2.2360680e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 2 + + final approximate solution + + -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 + + + + + problem 1 dimensions 5 50 + + + initial l2 norm of the residuals 8.0622577e+00 + + final l2 norm of the residuals 6.7082039e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 3 + + final approximate solution + + -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 -1.0000000e+00 + + + + + problem 2 dimensions 5 10 + + + initial l2 norm of the residuals 2.9152187e+02 + + final l2 norm of the residuals 1.4638501e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 8.0365911e+01 1.4011782e+02 -2.1423539e+02 7.0558910e+01 2.3706088e-03 + + + + + problem 2 dimensions 5 50 + + + initial l2 norm of the residuals 3.1016004e+03 + + final l2 norm of the residuals 3.4826302e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + -2.1546072e+02 -7.6406364e+01 -4.8102820e+01 -3.7703182e+01 1.3268487e+02 + + + + + problem 3 dimensions 5 10 + + + initial l2 norm of the residuals 1.2603968e+02 + + final l2 norm of the residuals 1.9097274e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 1.0954106e+02 -1.4666259e+02 5.5270531e+01 1.0000000e+00 + + + + + problem 3 dimensions 5 50 + + + initial l2 norm of the residuals 1.7489500e+03 + + final l2 norm of the residuals 3.6917294e+00 + + number of function evaluations 3 + + number of jacobian evaluations 2 + + exit parameter 1 + + final approximate solution + + 1.0000000e+00 -9.5807303e+01 -8.2600869e+01 1.0986204e+02 1.0000000e+00 + + + + + problem 4 dimensions 2 2 + + + initial l2 norm of the residuals 4.9193496e+00 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 21 + + number of jacobian evaluations 16 + + exit parameter 4 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 4 dimensions 2 2 + + + initial l2 norm of the residuals 1.3400631e+03 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 8 + + number of jacobian evaluations 5 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 4 dimensions 2 2 + + + initial l2 norm of the residuals 1.4300005e+05 + + final l2 norm of the residuals 0.0000000e+00 + + number of function evaluations 6 + + number of jacobian evaluations 4 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 + + + + + problem 5 dimensions 3 3 + + + initial l2 norm of the residuals 5.0000000e+01 + + final l2 norm of the residuals 1.4386145e-32 + + number of function evaluations 12 + + number of jacobian evaluations 9 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 -9.0390814e-34 0.0000000e+00 + + + + + problem 5 dimensions 3 3 + + + initial l2 norm of the residuals 1.0295630e+02 + + final l2 norm of the residuals 1.0446754e-19 + + number of function evaluations 20 + + number of jacobian evaluations 15 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 6.5638891e-21 0.0000000e+00 + + + + + problem 5 dimensions 3 3 + + + initial l2 norm of the residuals 9.9126182e+02 + + final l2 norm of the residuals 2.6163123e-29 + + number of function evaluations 19 + + number of jacobian evaluations 16 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 -1.6438775e-30 0.0000000e+00 + + + + + problem 6 dimensions 4 4 + + + initial l2 norm of the residuals 1.4662878e+01 + + final l2 norm of the residuals 1.4565773e-40 + + number of function evaluations 69 + + number of jacobian evaluations 68 + + exit parameter 4 + + final approximate solution + + 8.0669805e-21 -8.0669805e-22 1.2907169e-21 1.2907169e-21 + + + + + problem 6 dimensions 4 4 + + + initial l2 norm of the residuals 1.2709839e+03 + + final l2 norm of the residuals 2.2759020e-40 + + number of function evaluations 72 + + number of jacobian evaluations 71 + + exit parameter 4 + + final approximate solution + + 1.0083726e-20 -1.0083726e-21 1.6133961e-21 1.6133961e-21 + + + + + problem 6 dimensions 4 4 + + + initial l2 norm of the residuals 1.2688790e+05 + + final l2 norm of the residuals 2.2225605e-41 + + number of function evaluations 78 + + number of jacobian evaluations 77 + + exit parameter 4 + + final approximate solution + + 3.1511642e-21 -3.1511642e-22 5.0418628e-22 5.0418628e-22 + + + + + problem 7 dimensions 2 2 + + + initial l2 norm of the residuals 2.0012496e+01 + + final l2 norm of the residuals 6.9988752e+00 + + number of function evaluations 18 + + number of jacobian evaluations 11 + + exit parameter 1 + + final approximate solution + + 1.1412779e+01 -8.9680502e-01 + + + + + problem 7 dimensions 2 2 + + + initial l2 norm of the residuals 1.2432834e+04 + + final l2 norm of the residuals 6.9988752e+00 + + number of function evaluations 21 + + number of jacobian evaluations 13 + + exit parameter 1 + + final approximate solution + + 1.1412738e+01 -8.9680828e-01 + + + + + problem 7 dimensions 2 2 + + + initial l2 norm of the residuals 1.1426455e+07 + + final l2 norm of the residuals 6.9988752e+00 + + number of function evaluations 25 + + number of jacobian evaluations 17 + + exit parameter 1 + + final approximate solution + + 1.1412782e+01 -8.9680511e-01 + + + + + problem 8 dimensions 3 15 + + + initial l2 norm of the residuals 6.4561363e+00 + + final l2 norm of the residuals 9.0635960e-02 + + number of function evaluations 7 + + number of jacobian evaluations 6 + + exit parameter 1 + + final approximate solution + + 8.2410560e-02 1.1330361e+00 2.3436952e+00 + + + + + problem 8 dimensions 3 15 + + + initial l2 norm of the residuals 3.6141853e+01 + + final l2 norm of the residuals 4.1747687e+00 + + number of function evaluations 47 + + number of jacobian evaluations 46 + + exit parameter 1 + + final approximate solution + + 8.4066667e-01 -5.0830897e+09 -5.2601486e+09 + + + + + problem 8 dimensions 3 15 + + + initial l2 norm of the residuals 3.8411468e+02 + + final l2 norm of the residuals 4.1747687e+00 + + number of function evaluations 24 + + number of jacobian evaluations 23 + + exit parameter 1 + + final approximate solution + + 8.4066667e-01 -4.9730046e+09 -5.1462126e+09 + + + + + problem 9 dimensions 4 11 + + + initial l2 norm of the residuals 7.2891510e-02 + + final l2 norm of the residuals 1.7535838e-02 + + number of function evaluations 22 + + number of jacobian evaluations 20 + + exit parameter 1 + + final approximate solution + + 1.9280707e-01 1.9127927e-01 1.2305593e-01 1.3606091e-01 + + + + + problem 9 dimensions 4 11 + + + initial l2 norm of the residuals 2.9793701e+00 + + final l2 norm of the residuals 3.2052193e-02 + + number of function evaluations 88 + + number of jacobian evaluations 80 + + exit parameter 1 + + final approximate solution + + 2.3317387e+07 -1.4075883e+01 -1.0552795e+09 -6.5828474e+08 + + + + + problem 9 dimensions 4 11 + + + initial l2 norm of the residuals 2.9959062e+01 + + final l2 norm of the residuals 1.7535838e-02 + + number of function evaluations 496 + + number of jacobian evaluations 381 + + exit parameter 1 + + final approximate solution + + 1.9280678e-01 1.9128581e-01 1.2305716e-01 1.3606394e-01 + + + + + problem 10 dimensions 3 16 + + + initial l2 norm of the residuals 4.1153467e+04 + + final l2 norm of the residuals 9.3779451e+00 + + number of function evaluations 126 + + number of jacobian evaluations 116 + + exit parameter 3 + + final approximate solution + + 5.6096365e-03 6.1813463e+03 3.4522363e+02 + + + + + problem 10 dimensions 3 16 + + + initial l2 norm of the residuals 4.1682169e+06 + + final l2 norm of the residuals 7.9968310e+02 + + number of function evaluations 400 + + number of jacobian evaluations 347 + + exit parameter 5 + + final approximate solution + + 9.5646267e-12 3.4476558e+04 9.1250720e+02 + + + + + problem 11 dimensions 6 31 + + + initial l2 norm of the residuals 5.4772256e+00 + + final l2 norm of the residuals 4.7829594e-02 + + number of function evaluations 9 + + number of jacobian evaluations 8 + + exit parameter 1 + + final approximate solution + + -1.5725102e-02 1.0124349e+00 -2.3299161e-01 1.2604300e+00 -1.5137287e+00 + 9.9299633e-01 + + + + + problem 11 dimensions 6 31 + + + initial l2 norm of the residuals 6.4331258e+03 + + final l2 norm of the residuals 4.7829594e-02 + + number of function evaluations 15 + + number of jacobian evaluations 14 + + exit parameter 1 + + final approximate solution + + -1.5725073e-02 1.0124349e+00 -2.3299164e-01 1.2604302e+00 -1.5137291e+00 + 9.9299652e-01 + + + + + problem 11 dimensions 6 31 + + + initial l2 norm of the residuals 6.7425604e+05 + + final l2 norm of the residuals 4.7829594e-02 + + number of function evaluations 16 + + number of jacobian evaluations 15 + + exit parameter 1 + + final approximate solution + + -1.5725135e-02 1.0124349e+00 -2.3299159e-01 1.2604297e+00 -1.5137284e+00 + 9.9299610e-01 + + + + + problem 11 dimensions 9 31 + + + initial l2 norm of the residuals 5.4772256e+00 + + final l2 norm of the residuals 1.1831146e-03 + + number of function evaluations 8 + + number of jacobian evaluations 7 + + exit parameter 1 + + final approximate solution + + -1.5307064e-05 9.9978970e-01 1.4763963e-02 1.4634233e-01 1.0008211e+00 + -2.6177311e+00 4.1044031e+00 -3.1436123e+00 1.0526264e+00 + + + + + problem 11 dimensions 9 31 + + + initial l2 norm of the residuals 1.2088127e+04 + + final l2 norm of the residuals 1.1831146e-03 + + number of function evaluations 19 + + number of jacobian evaluations 15 + + exit parameter 1 + + final approximate solution + + -1.5307036e-05 9.9978970e-01 1.4763964e-02 1.4634233e-01 1.0008211e+00 + -2.6177311e+00 4.1044032e+00 -3.1436123e+00 1.0526264e+00 + + + + + problem 11 dimensions 9 31 + + + initial l2 norm of the residuals 1.2691093e+06 + + final l2 norm of the residuals 1.1831146e-03 + + number of function evaluations 18 + + number of jacobian evaluations 15 + + exit parameter 1 + + final approximate solution + + -1.5307037e-05 9.9978970e-01 1.4763964e-02 1.4634233e-01 1.0008211e+00 + -2.6177311e+00 4.1044032e+00 -3.1436123e+00 1.0526264e+00 + + + + + problem 11 dimensions 12 31 + + + initial l2 norm of the residuals 5.4772256e+00 + + final l2 norm of the residuals 2.1731040e-05 + + number of function evaluations 9 + + number of jacobian evaluations 8 + + exit parameter 1 + + final approximate solution + + -6.6380605e-09 1.0000016e+00 -5.6393221e-04 3.4782054e-01 -1.5673150e-01 + 1.0528152e+00 -3.2472712e+00 7.2884349e+00 -1.0271848e+01 9.0741136e+00 + -4.5413755e+00 1.0120119e+00 + + + + + problem 11 dimensions 12 31 + + + initial l2 norm of the residuals 1.9220759e+04 + + final l2 norm of the residuals 2.1731040e-05 + + number of function evaluations 14 + + number of jacobian evaluations 13 + + exit parameter 2 + + final approximate solution + + -6.6380605e-09 1.0000016e+00 -5.6393221e-04 3.4782054e-01 -1.5673150e-01 + 1.0528152e+00 -3.2472712e+00 7.2884349e+00 -1.0271848e+01 9.0741136e+00 + -4.5413755e+00 1.0120119e+00 + + + + + problem 11 dimensions 12 31 + + + initial l2 norm of the residuals 2.0189180e+06 + + final l2 norm of the residuals 2.1731040e-05 + + number of function evaluations 34 + + number of jacobian evaluations 28 + + exit parameter 3 + + final approximate solution + + -6.6379511e-09 1.0000016e+00 -5.6393221e-04 3.4782054e-01 -1.5673150e-01 + 1.0528152e+00 -3.2472712e+00 7.2884349e+00 -1.0271848e+01 9.0741136e+00 + -4.5413755e+00 1.0120119e+00 + + + + + problem 12 dimensions 3 10 + + + initial l2 norm of the residuals 3.2111584e+01 + + final l2 norm of the residuals 1.1345334e-16 + + number of function evaluations 8 + + number of jacobian evaluations 7 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+01 1.0000000e+00 + + + + + problem 13 dimensions 2 10 + + + initial l2 norm of the residuals 6.4585650e+01 + + final l2 norm of the residuals 1.1151779e+01 + + number of function evaluations 24 + + number of jacobian evaluations 13 + + exit parameter 1 + + final approximate solution + + 2.5782487e-01 2.5782554e-01 + + + + + problem 14 dimensions 4 20 + + + initial l2 norm of the residuals 2.8154384e+03 + + final l2 norm of the residuals 2.9295427e+02 + + number of function evaluations 301 + + number of jacobian evaluations 279 + + exit parameter 1 + + final approximate solution + + -1.1594175e+01 1.3203535e+01 -4.0345582e-01 2.3677057e-01 + + + + + problem 14 dimensions 4 20 + + + initial l2 norm of the residuals 5.5507335e+05 + + final l2 norm of the residuals 2.9295427e+02 + + number of function evaluations 68 + + number of jacobian evaluations 54 + + exit parameter 1 + + final approximate solution + + -1.1594593e+01 1.3203687e+01 -4.0341581e-01 2.3679328e-01 + + + + + problem 14 dimensions 4 20 + + + initial l2 norm of the residuals 6.1211252e+07 + + final l2 norm of the residuals 2.9295427e+02 + + number of function evaluations 310 + + number of jacobian evaluations 289 + + exit parameter 1 + + final approximate solution + + -1.1593993e+01 1.3203462e+01 -4.0345948e-01 2.3677112e-01 + + + + + problem 15 dimensions 1 8 + + + initial l2 norm of the residuals 1.8862380e+00 + + final l2 norm of the residuals 1.8862380e+00 + + number of function evaluations 1 + + number of jacobian evaluations 1 + + exit parameter 4 + + final approximate solution + + 5.0000000e-01 + + + + + problem 15 dimensions 1 8 + + + initial l2 norm of the residuals 5.3833444e+09 + + final l2 norm of the residuals 1.8842482e+00 + + number of function evaluations 30 + + number of jacobian evaluations 29 + + exit parameter 1 + + final approximate solution + + 9.8173132e-01 + + + + + problem 15 dimensions 1 8 + + + initial l2 norm of the residuals 1.1808873e+18 + + final l2 norm of the residuals 1.8842482e+00 + + number of function evaluations 48 + + number of jacobian evaluations 47 + + exit parameter 1 + + final approximate solution + + 9.8173132e-01 + + + + + problem 15 dimensions 8 8 + + + initial l2 norm of the residuals 1.9651386e-01 + + final l2 norm of the residuals 5.9303235e-02 + + number of function evaluations 45 + + number of jacobian evaluations 23 + + exit parameter 1 + + final approximate solution + + 4.3152867e-02 1.9309093e-01 2.6632869e-01 5.0000008e-01 4.9999992e-01 + 7.3367131e-01 8.0690907e-01 9.5684713e-01 + + + + + problem 15 dimensions 9 9 + + + initial l2 norm of the residuals 1.6994993e-01 + + final l2 norm of the residuals 1.9152345e-19 + + number of function evaluations 12 + + number of jacobian evaluations 9 + + exit parameter 2 + + final approximate solution + + 4.4205346e-02 1.9949067e-01 2.3561911e-01 4.1604691e-01 5.0000000e-01 + 5.8395309e-01 7.6438089e-01 8.0050933e-01 9.5579465e-01 + + + + + problem 15 dimensions 10 10 + + + initial l2 norm of the residuals 1.8374783e-01 + + final l2 norm of the residuals 8.0647100e-02 + + number of function evaluations 29 + + number of jacobian evaluations 14 + + exit parameter 1 + + final approximate solution + + 5.9619952e-02 1.6670835e-01 2.3917071e-01 3.9888443e-01 3.9888421e-01 + 6.0111579e-01 6.0111557e-01 7.6082929e-01 8.3329165e-01 9.4038005e-01 + + + + + problem 16 dimensions 10 10 + + + initial l2 norm of the residuals 1.6530216e+01 + + final l2 norm of the residuals 6.1801901e-18 + + number of function evaluations 15 + + number of jacobian evaluations 13 + + exit parameter 2 + + final approximate solution + + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 1.2056970e+00 + + + + + problem 16 dimensions 10 10 + + + initial l2 norm of the residuals 9.7656240e+06 + + final l2 norm of the residuals 2.6043430e-18 + + number of function evaluations 13 + + number of jacobian evaluations 8 + + exit parameter 2 + + final approximate solution + + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 1.2056970e+00 + + + + + problem 16 dimensions 10 10 + + + initial l2 norm of the residuals 9.7656250e+16 + + final l2 norm of the residuals 1.8407505e-18 + + number of function evaluations 30 + + number of jacobian evaluations 27 + + exit parameter 2 + + final approximate solution + + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 + 9.7943030e-01 9.7943030e-01 9.7943030e-01 9.7943030e-01 1.2056970e+00 + + + + + problem 16 dimensions 30 30 + + + initial l2 norm of the residuals 8.3476044e+01 + + final l2 norm of the residuals 5.3715626e-17 + + number of function evaluations 19 + + number of jacobian evaluations 14 + + exit parameter 2 + + final approximate solution + + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 + 9.9775422e-01 9.9775422e-01 9.9775422e-01 9.9775422e-01 1.0673735e+00 + + + + + problem 16 dimensions 40 40 + + + initial l2 norm of the residuals 1.2802636e+02 + + final l2 norm of the residuals 1.4828176e-16 + + number of function evaluations 19 + + number of jacobian evaluations 14 + + exit parameter 2 + + final approximate solution + + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 1.0000000e+00 + + + + + problem 17 dimensions 5 33 + + + initial l2 norm of the residuals 9.3756402e-01 + + final l2 norm of the residuals 7.3924926e-03 + + number of function evaluations 19 + + number of jacobian evaluations 16 + + exit parameter 1 + + final approximate solution + + 3.7541005e-01 1.9358469e+00 -1.4646871e+00 1.2867535e-02 2.2122700e-02 + + + + + problem 18 dimensions 11 65 + + + initial l2 norm of the residuals 1.4468654e+00 + + final l2 norm of the residuals 2.0034404e-01 + + number of function evaluations 17 + + number of jacobian evaluations 13 + + exit parameter 1 + + final approximate solution + + 1.3099771e+00 4.3155362e-01 6.3366163e-01 5.9943049e-01 7.5418274e-01 + 9.0428996e-01 1.3658118e+00 4.8236970e+00 2.3986850e+00 4.5688745e+00 + 5.6753415e+00 + summary of 53 calls to lmstr1 + nprob n m nfev njev info final L2 norm + + 1 5 10 3 2 2 2.2360680e+00 + 1 5 50 3 2 3 6.7082039e+00 + 2 5 10 3 2 1 1.4638501e+00 + 2 5 50 3 2 1 3.4826302e+00 + 3 5 10 3 2 1 1.9097274e+00 + 3 5 50 3 2 1 3.6917294e+00 + 4 2 2 21 16 4 0.0000000e+00 + 4 2 2 8 5 2 0.0000000e+00 + 4 2 2 6 4 2 0.0000000e+00 + 5 3 3 12 9 2 1.4386145e-32 + 5 3 3 20 15 2 1.0446754e-19 + 5 3 3 19 16 2 2.6163123e-29 + 6 4 4 69 68 4 1.4565773e-40 + 6 4 4 72 71 4 2.2759020e-40 + 6 4 4 78 77 4 2.2225605e-41 + 7 2 2 18 11 1 6.9988752e+00 + 7 2 2 21 13 1 6.9988752e+00 + 7 2 2 25 17 1 6.9988752e+00 + 8 3 15 7 6 1 9.0635960e-02 + 8 3 15 47 46 1 4.1747687e+00 + 8 3 15 24 23 1 4.1747687e+00 + 9 4 11 22 20 1 1.7535838e-02 + 9 4 11 88 80 1 3.2052193e-02 + 9 4 11 496 381 1 1.7535838e-02 + 10 3 16 126 116 3 9.3779451e+00 + 10 3 16 400 347 5 7.9968310e+02 + 11 6 31 9 8 1 4.7829594e-02 + 11 6 31 15 14 1 4.7829594e-02 + 11 6 31 16 15 1 4.7829594e-02 + 11 9 31 8 7 1 1.1831146e-03 + 11 9 31 19 15 1 1.1831146e-03 + 11 9 31 18 15 1 1.1831146e-03 + 11 12 31 9 8 1 2.1731040e-05 + 11 12 31 14 13 2 2.1731040e-05 + 11 12 31 34 28 3 2.1731040e-05 + 12 3 10 8 7 2 1.1345334e-16 + 13 2 10 24 13 1 1.1151779e+01 + 14 4 20 301 279 1 2.9295427e+02 + 14 4 20 68 54 1 2.9295427e+02 + 14 4 20 310 289 1 2.9295427e+02 + 15 1 8 1 1 4 1.8862380e+00 + 15 1 8 30 29 1 1.8842482e+00 + 15 1 8 48 47 1 1.8842482e+00 + 15 8 8 45 23 1 5.9303235e-02 + 15 9 9 12 9 2 1.9152345e-19 + 15 10 10 29 14 1 8.0647100e-02 + 16 10 10 15 13 2 6.1801901e-18 + 16 10 10 13 8 2 2.6043430e-18 + 16 10 10 30 27 2 1.8407505e-18 + 16 30 30 19 14 2 5.3715626e-17 + 16 40 40 19 14 2 1.4828176e-16 + 17 5 33 19 16 1 7.3924926e-03 + 18 11 65 17 13 1 2.0034404e-01 diff --git a/examples/ref/ldtchkderc.ref b/examples/ref/ldtchkderc.ref new file mode 100644 index 0000000..32faf40 --- /dev/null +++ b/examples/ref/ldtchkderc.ref @@ -0,0 +1,22 @@ + + fvec + + -1.181606 -1.429655 -1.606344 + -1.745269 -1.840654 -1.921586 + -1.984141 -2.022537 -2.468977 + -2.827562 -3.473582 -4.437612 + -6.047662 -9.267761 -18.91806 + fvecp - fvec + + -1.706926e-10 -7.584616e-11 -4.496407e-12 + 5.112571e-11 9.570422e-11 1.322311e-10 + 1.627071e-10 1.885207e-10 3.289353e-10 + 5.161547e-10 7.78262e-10 1.171423e-09 + 1.826691e-09 3.137227e-09 7.068836e-09 + err + + 0.09421056 0.08207347 0.0798528 + 0.08237829 0.08857469 0.1007351 + 0.1260227 1 1 + 1 1 1 + 1 1 1 diff --git a/examples/ref/ldtchkderc_box.ref b/examples/ref/ldtchkderc_box.ref new file mode 100644 index 0000000..5a80825 --- /dev/null +++ b/examples/ref/ldtchkderc_box.ref @@ -0,0 +1,22 @@ + + fvec + + -0.9829917 -1.111746 -1.214351 + -1.306448 -1.372163 -1.434587 + -1.486078 -1.518468 -1.820838 + -1.987333 -2.364425 -2.925063 + -3.862792 -5.738252 -11.35463 + fvecp - fvec + + -2.85884e-10 -2.71445e-10 -2.579249e-10 + -2.453894e-10 -2.338206e-10 -2.231626e-10 + -2.133445e-10 -2.042922e-10 -1.766614e-10 + -1.398203e-10 -8.824282e-11 -1.087659e-11 + 1.180671e-10 3.759546e-10 1.149617e-09 + err + + 0.1975415 0.1758616 0.1672297 + 0.1651403 0.1676895 0.1769975 + 0.1999344 0.789874 0.7864705 + 0.7788527 0.7739175 0.7697708 + 0.765575 0.7618356 0.7582316 diff --git a/examples/ref/ldtenormc.ref b/examples/ref/ldtenormc.ref new file mode 100644 index 0000000..2073dd6 --- /dev/null +++ b/examples/ref/ldtenormc.ref @@ -0,0 +1,6 @@ +dpmpar(2) = 3.362103e-4932, dpmpar(3) = 1.189731e+4932 +rdwarf = 2.245696932951581572e-2465, rgiant = 1.090748135619415929e+2465 +norm/rdwarf (naive) = 0.999999999999997851 +norm/rdwarf (enorm) = 1 +norm/agiant (naive) = 1.00000000000000677 +norm/agiant (enorm) = 0.999999999999998332 diff --git a/examples/ref/ldtfdjac2c.ref b/examples/ref/ldtfdjac2c.ref new file mode 100644 index 0000000..cc712e2 --- /dev/null +++ b/examples/ref/ldtfdjac2c.ref @@ -0,0 +1,29 @@ + + fvec + + -1.181606 -1.429655 -1.606344 + -1.745269 -1.840654 -1.921586 + -1.984141 -2.022537 -2.468977 + -2.827562 -3.473582 -4.437612 + -6.047662 -9.267761 -18.91806 + fvecp - fvec + + -1.706926e-10 -7.584616e-11 -4.496407e-12 + 5.112571e-11 9.570422e-11 1.322311e-10 + 1.627071e-10 1.885207e-10 3.289353e-10 + 5.161547e-10 7.78262e-10 1.171423e-09 + 1.826691e-09 3.137227e-09 7.068836e-09 + errd + + 1 1 1 + 1 0.9961974 0.9981682 + 1 1 1 + 1 1 1 + 1 1 1 + err + + 1 1 1 + 1 1 1 + 1 1 1 + 1 1 1 + 1 1 1 diff --git a/examples/ref/ldthybrd1c.ref b/examples/ref/ldthybrd1c.ref new file mode 100644 index 0000000..12023ec --- /dev/null +++ b/examples/ref/ldthybrd1c.ref @@ -0,0 +1,7 @@ + final L2 norm of the residuals 6.178398e-11 + exit parameter 1 + final approximates solution + + -0.5706545 -0.6816283 -0.7017325 + -0.7042129 -0.701369 -0.6918656 + -0.665792 -0.5960342 -0.4164121 diff --git a/examples/ref/ldthybrdc.ref b/examples/ref/ldthybrdc.ref new file mode 100644 index 0000000..377f3f1 --- /dev/null +++ b/examples/ref/ldthybrdc.ref @@ -0,0 +1,11 @@ + final l2 norm of the residuals 6.178398e-11 + + number of function evaluations 16 + + exit parameter 1 + + final approximate solution + + -0.5706545 -0.6816283 -0.7017325 + -0.7042129 -0.701369 -0.6918656 + -0.665792 -0.5960342 -0.4164121 diff --git a/examples/ref/ldthybrj1c.ref b/examples/ref/ldthybrj1c.ref new file mode 100644 index 0000000..0a6e591 --- /dev/null +++ b/examples/ref/ldthybrj1c.ref @@ -0,0 +1,9 @@ + final l2 norm of the residuals 6.178398e-11 + + exit parameter 1 + + final approximate solution + + -0.5706545 -0.6816283 -0.7017325 + -0.7042129 -0.701369 -0.6918656 + -0.665792 -0.5960342 -0.4164121 diff --git a/examples/ref/ldthybrjc.ref b/examples/ref/ldthybrjc.ref new file mode 100644 index 0000000..98f6871 --- /dev/null +++ b/examples/ref/ldthybrjc.ref @@ -0,0 +1,14 @@ + final l2 norm of the residuals 6.178398e-11 + + number of function evaluations 13 + + number of jacobian evaluations 1 + + exit parameter 1 + + final approximate solution + + + -0.5706545 -0.6816283 -0.7017325 + -0.7042129 -0.701369 -0.6918656 + -0.665792 -0.5960342 -0.4164121 diff --git a/examples/ref/ldthybrjc_box.ref b/examples/ref/ldthybrjc_box.ref new file mode 100644 index 0000000..fda26ac --- /dev/null +++ b/examples/ref/ldthybrjc_box.ref @@ -0,0 +1,14 @@ + final l2 norm of the residuals 0.7585268 + + number of function evaluations 21 + + number of jacobian evaluations 2 + + exit parameter 5 + + final approximate solution + + + -0.4715775 -0.5 -0.6075997 + -0.6912395 -0.715058 -0.7148161 + -0.689358 -0.613326 -0.4219724 diff --git a/examples/ref/ldtlmder1c.ref b/examples/ref/ldtlmder1c.ref new file mode 100644 index 0000000..68f478c --- /dev/null +++ b/examples/ref/ldtlmder1c.ref @@ -0,0 +1,7 @@ + final l2 norm of the residuals 0.09063596 + + exit parameter 1 + + final approximate solution + + 0.08241056 1.133036 2.343695 diff --git a/examples/ref/ldtlmderc.ref b/examples/ref/ldtlmderc.ref new file mode 100644 index 0000000..2e8f321 --- /dev/null +++ b/examples/ref/ldtlmderc.ref @@ -0,0 +1,16 @@ + final l2 norm of the residuals 0.09063596 + + number of function evaluations 7 + + number of Jacobian evaluations 6 + + exit parameter 1 + + final approximate solution + + 0.08241056 1.133036 2.343695 + covariance + + 0.0001531199 0.002869831 -0.002656551 + 0.002869831 0.09480247 -0.09098321 + -0.002656551 -0.09098321 0.08778068 diff --git a/examples/ref/ldtlmderc_box.ref b/examples/ref/ldtlmderc_box.ref new file mode 100644 index 0000000..d1a91e5 --- /dev/null +++ b/examples/ref/ldtlmderc_box.ref @@ -0,0 +1,16 @@ + final l2 norm of the residuals 0.09071386 + + number of function evaluations 10 + + number of Jacobian evaluations 9 + + exit parameter 1 + + final approximate solution + + 0.08371095 1.178309 2.3 + covariance + + 0.0001536684 0.002994926 -0 + 0.002994926 0.1026389 -0 + -0 -0 0 diff --git a/examples/ref/ldtlmdif1c.ref b/examples/ref/ldtlmdif1c.ref new file mode 100644 index 0000000..faa893f --- /dev/null +++ b/examples/ref/ldtlmdif1c.ref @@ -0,0 +1,7 @@ + final l2 norm of the residuals 0.09063596 + + exit parameter 1 + + final approximate solution + + 0.08241056 1.133036 2.343695 diff --git a/examples/ref/ldtlmdifc.ref b/examples/ref/ldtlmdifc.ref new file mode 100644 index 0000000..9e20c02 --- /dev/null +++ b/examples/ref/ldtlmdifc.ref @@ -0,0 +1,14 @@ + final l2 norm of the residuals 0.09063596 + + number of function evaluations 25 + + exit parameter 1 + + final approximate solution + + 0.08241056 1.133036 2.343695 + covariance + + 0.0001531199 0.002869831 -0.002656551 + 0.002869831 0.09480247 -0.09098321 + -0.002656551 -0.09098321 0.08778068 diff --git a/examples/ref/ldtlmstr1c.ref b/examples/ref/ldtlmstr1c.ref new file mode 100644 index 0000000..68f478c --- /dev/null +++ b/examples/ref/ldtlmstr1c.ref @@ -0,0 +1,7 @@ + final l2 norm of the residuals 0.09063596 + + exit parameter 1 + + final approximate solution + + 0.08241056 1.133036 2.343695 diff --git a/examples/ref/ldtlmstrc.ref b/examples/ref/ldtlmstrc.ref new file mode 100644 index 0000000..2e8f321 --- /dev/null +++ b/examples/ref/ldtlmstrc.ref @@ -0,0 +1,16 @@ + final l2 norm of the residuals 0.09063596 + + number of function evaluations 7 + + number of Jacobian evaluations 6 + + exit parameter 1 + + final approximate solution + + 0.08241056 1.133036 2.343695 + covariance + + 0.0001531199 0.002869831 -0.002656551 + 0.002869831 0.09480247 -0.09098321 + -0.002656551 -0.09098321 0.08778068 diff --git a/examples/tenorm_.c b/examples/tenorm_.c index ecbbfb4..79fe675 100644 --- a/examples/tenorm_.c +++ b/examples/tenorm_.c @@ -17,8 +17,9 @@ #ifdef __cminpack_long_double__ #define N1 3 #define N2 367 -#define SQRTFAC 1.5 -#define FAC 10. +#define SQRTFAC 1.5L +#define FAC 10.L +#define sqrt(x) sqrtl(x) #endif #ifdef __cminpack_float__ #define N1 3 diff --git a/examples/tenormc.c b/examples/tenormc.c index f3e36b0..bd3941b 100644 --- a/examples/tenormc.c +++ b/examples/tenormc.c @@ -17,8 +17,8 @@ #ifdef __cminpack_long_double__ #define N1 3 #define N2 367 -#define SQRTFAC 1.5 -#define FAC 10. +#define SQRTFAC 1.5L +#define FAC 10.L #define sqrt(x) sqrtl(x) #endif #ifdef __cminpack_float__