From 0b1f9545c77d3c21d43f1f0080fbe0bae95c3413 Mon Sep 17 00:00:00 2001
From: Jiri Pittner <jiri@pittnerovi.com>
Date: Fri, 23 Oct 2026 14:45:49 +0200
Subject: [PATCH] made compilable with gcc16

---
 t.cc      | 33 +++++++++++++++++++++++++++++----
 vecmat3.h |  3 +++
 2 files changed, 32 insertions(+), 4 deletions(-)

diff --git a/t.cc b/t.cc
index 169a746..de9d2fd 100644
--- a/t.cc
+++ b/t.cc
@@ -2461,7 +2461,7 @@ Quaternion<double> q({1,2,3,4});
 cout<<q<<endl;
 NRPerm<int> p({1,4,2,3});
 cout <<p<<endl;
-CyclePerm<int> c({NRVec_from1({1,2}),NRVec_from1({3,4,5})});
+CyclePerm<int> c({NRVec_from1<int>({1,2}),NRVec_from1<int>({3,4,5})});
 cout <<c;
 Polynomial<double> pp({1,2,3,4,5});
 cout<<pp;
@@ -2518,7 +2518,7 @@ NRMat<double> m(&mm[0][0],3,3);
 cout <<m<<"3 3\n"<<mm<<endl;
 
 double rr[2][3]={{1,2,3},{4,5,6}};
-NRMat r(rr);
+NRMat<double> r(rr);
 cout<<r;
 
 double d;
@@ -4645,7 +4645,7 @@ cout <<r;
 if(n<=20)
 	{
 	cout<<aa;
-	NRMat a(aa);
+	NRMat<double> a(aa);
 	NRVec<double> b(n);
 	cout <<a;
 	diagonalize(a,b);
@@ -4809,7 +4809,7 @@ if(p!=q) laerror("inverseme failed");
 else cout <<p<<"OK\n";
 }
 
-if(1)
+if(0)
 {
 int n;
 cin >>n;
@@ -4830,5 +4830,30 @@ cout <<"adjoint error = "<<(a*adj).norm(abs(det))<<endl;
 
 }
 
+if(1)
+{
+int seed;
+int f=open("/dev/random",O_RDONLY);
+if(sizeof(int)!=read(f,&seed,sizeof(int))) laerror("cannot read /dev/random");
+close(f);
+srand(seed);
+
+int n;
+cin >>n;
+NRVec<double> x(n);
+x.randomize(1.);
+testproblem<double> problem;
+BFGS_hyperparameters<double> param;
+param.my_maxstep=100;
+param.my_f_thr=1e-6;
+param.my_x_thr=1e-3;
+param.my_g_thr=1e-3;
+cout << "initial " <<x<<endl;
+BFGS_Optimizer<double> o;
+o.minimize(&problem,x,&param);
+cout << "optimized " <<x<<endl;
+double r=problem.funcgrad(x,NULL);
+cout <<r<<endl;
+}
 
 }//main
diff --git a/vecmat3.h b/vecmat3.h
index faa3abc..d1a3333 100644
--- a/vecmat3.h
+++ b/vecmat3.h
@@ -30,6 +30,9 @@
 #include <string.h>
 #include <math.h>
 #include <stdio.h>
+#ifndef NO_NUMERIC_LIMITS
+#include <limits>
+#endif
 
 namespace LA_Vecmat3 {