From 3bfa3007cee237035dbb8bd324d03d4120c4163f Mon Sep 17 00:00:00 2001
From: Jiri Pittner <jiri@pittnerovi.com>
Date: Mon, 1 Jun 2026 15:45:36 +0200
Subject: [PATCH] bugfix for compatibility with gcc16

---
 mat.h  | 2 +-
 smat.h | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/mat.h b/mat.h
index 8355df5..06d8862 100644
--- a/mat.h
+++ b/mat.h
@@ -1634,7 +1634,7 @@ for(int i=0; i<b.nrows(); ++i)
 template<typename T>
 void NRMat<T>::printsorted(std::ostream &s, int direction, bool absvalue, bool stable, typename LA_traits<T>::normtype thr, int topmax) const                                                                                                                                {
 NRMat<int> indexmap(nn*mm,2);
-NRVec tmp(*this);
+NRVec<T> tmp(*this);
 int ii=0;
 for(int i=0; i<nn; ++i) for(int j=0; j<mm; ++j) {indexmap[ii][0]=i; indexmap[ii][1]=j; ii++;}
 tmp.printsorted(s,direction,absvalue,stable,thr,topmax,&indexmap);
diff --git a/smat.h b/smat.h
index 94152d9..9347040 100644
--- a/smat.h
+++ b/smat.h
@@ -1293,7 +1293,7 @@ template<typename T>
 void NRSMat<T>::printsorted(std::ostream &s, int direction, bool absvalue, bool stable, typename LA_traits<T>::normtype thr, int topmax) const
 {
 NRMat<int> indexmap(NN2,2);
-NRVec tmp(*this);
+NRVec<T> tmp(*this);
 int ii=0;
 for(int i=0; i<nn; ++i) for(int j=0; j<=i; ++j) {indexmap[ii][1]=i; indexmap[ii][0]=j; ii++;}
 tmp.printsorted(s,direction,absvalue,stable,thr,topmax,&indexmap);