From 2621f444e1791199ec19ec35fb786b5ea303ada1 Mon Sep 17 00:00:00 2001
From: Jiri Pittner <jiri@pittnerovi.com>
Date: Tue, 21 Jun 2022 17:33:16 +0200
Subject: [PATCH] fixed bug in nrvec::concatme

---
 t.cc  | 13 +++++++++++--
 vec.h |  8 ++++++--
 2 files changed, 17 insertions(+), 4 deletions(-)

diff --git a/t.cc b/t.cc
index 68f67c9..18738e3 100644
--- a/t.cc
+++ b/t.cc
@@ -2539,7 +2539,7 @@ cout<<"big "<<xx1<<endl;
 cout << "TEST "<<(xx<xx1) <<" "<<(xx>xx2) <<endl;
 }
 
-if(1)
+if(0)
 {
 int seed;
 int f=open("/dev/random",O_RDONLY);
@@ -2574,6 +2574,15 @@ cout <<test;
 cout <<"Error = "<<(expitszsz-test).norm()<<endl;
 }
 
-
+if(1)
+{
+NRVec<double> x({1,2,3});
+NRVec<double> y({4,5,6});
+//cout <<x.concat(y);
+x.append(10.);
+cout <<x;
+x.concatme(y);
+cout <<x;
+}
 
 }
diff --git a/vec.h b/vec.h
index 583bf0c..f9789b1 100644
--- a/vec.h
+++ b/vec.h
@@ -244,14 +244,18 @@ public:
 	//!concatenate vector into *this 
 	void concatme(const NRVec &rhs) 
 		{
+		NOT_GPU(*this);
 		if(rhs.nn==0) return;
+		int nnold=nn;
 		resize(nn+rhs.nn,true);
-		for(int i=0; i<rhs.nn; ++i) *(this)[nn+i] = rhs[i];
+		for(int i=0; i<rhs.nn; ++i) v[nnold+i] = rhs[i];
 		}
 	void append(const T &a) //not efficient if done repeatedly
 		{
+		NOT_GPU(*this);
+		int nnold=nn;
                 resize(nn+1,true);
-                (*this)[nn] = a;
+                v[nnold] = a;
                 }
 
 	//! determine the actual value of the reference counter