Print reference results

certik · Nov 4, 2022 · 0cd235c · 0cd235c
1 parent 123becc
commit 0cd235c
Showing 1 changed file with 8 additions and 1 deletion.
diff --git a/examples/example_hybrd.f90 b/examples/example_hybrd.f90
@@ -49,7 +49,14 @@ program example_hybrd
                mode, factor, nprint, info, nfev, fjac, ldfjac, &
                r, lr, qtf, wa1, wa2, wa3, wa4)
     fnorm = enorm(n, fvec)
-    write (nwrite, 1000) fnorm, nfev, info !, (x(j), j=1, n)
+    !write (nwrite, 1000) fnorm, nfev, info, (x(j), j=1, n)
+    print *, fnorm, nfev, info
+    print *, "Computed:"
+    print *, x
+    print *, "Reference:"
+    print *, -0.5706545D+00, -0.6816283D+00, -0.7017325D+00, &
+             -0.7042129D+00, -0.7013690D+00, -0.6918656D+00, &
+             -0.6657920D+00, -0.5960342D+00, -0.4164121D+00
 
 1000 format(5x, "FINAL L2 NORM OF THE RESIDUALS", d15.7// &
            5x, "NUMBER OF FUNCTION EVALUATIONS", i10// &