diff -r -b -B -c1 cpp/disk.id ../cpp/disk.id
*** cpp/disk.id	Sat Dec  8 16:49:48 2001
--- ../cpp/disk.id	Sun Oct 20 12:42:08 2002
***************
*** 1,3 ****
  Numerical Recipes in ANSI C++
! 2.10
  Disk 1 of 1
--- 1,3 ----
  Numerical Recipes in ANSI C++
! 2.11
  Disk 1 of 1
diff -r -b -B -c1 cpp/examples/xfactrl.cpp ../cpp/examples/xfactrl.cpp
*** cpp/examples/xfactrl.cpp	Sun Nov 25 15:27:40 2001
--- ../cpp/examples/xfactrl.cpp	Sun Oct  6 17:43:10 2002
***************
*** 32,40 ****
            fp >> n >> actual;
-           if (actual < 1.0e10) {
              cout << setw(6) << n << setw(16) << actual;
              cout << setw(16) << NR::factrl(n) << endl;
-           } else {
-             cout << setw(6) << n << setw(16) << actual;
-             cout << setw(16) << NR::factrl(n) << endl;
-           }
          }
--- 32,35 ----
diff -r -b -B -c1 cpp/examples/xfrenel.cpp ../cpp/examples/xfrenel.cpp
*** cpp/examples/xfrenel.cpp	Sun Nov 25 15:27:40 2001
--- ../cpp/examples/xfrenel.cpp	Sun Oct  6 17:43:18 2002
***************
*** 35,37 ****
            NR::frenel(x,cs);
!           cout << setw(15) << x << setw(15) << xs << setw(15) << imag(cs);;
            cout << setw(15) << xc << setw(15) << real(cs)<< endl;
--- 35,37 ----
            NR::frenel(x,cs);
!           cout << setw(15) << x << setw(15) << xs << setw(15) << imag(cs);
            cout << setw(15) << xc << setw(15) << real(cs)<< endl;
diff -r -b -B -c1 cpp/examples/xlubksb.cpp ../cpp/examples/xlubksb.cpp
*** cpp/examples/xlubksb.cpp	Sun Nov 25 15:27:42 2001
--- ../cpp/examples/xlubksb.cpp	Sun Oct  6 17:43:46 2002
***************
*** 26,28 ****
            Vec_DP x(n);
!           Mat_DP a(n,n),b(n,n),c(n,n);
            getline(fp,txt);
--- 26,28 ----
            Vec_DP x(n);
!           Mat_DP a(n,n),b(n,m),c(n,n);
            getline(fp,txt);
diff -r -b -B -c1 cpp/examples/xpowell.cpp ../cpp/examples/xpowell.cpp
*** cpp/examples/xpowell.cpp	Sun Nov 25 15:27:42 2001
--- ../cpp/examples/xpowell.cpp	Sun Oct  6 17:43:28 2002
***************
*** 27,29 ****
          NR::powell(p,xi,FTOL,iter,fret,func);
!         cout << "Iterations: " << iter << endl << endl;;
          cout << "Minimum found at: " << endl;
--- 27,29 ----
          NR::powell(p,xi,FTOL,iter,fret,func);
!         cout << "Iterations: " << iter << endl << endl;
          cout << "Minimum found at: " << endl;
diff -r -b -B -c1 cpp/examples/xrealft.cpp ../cpp/examples/xrealft.cpp
*** cpp/examples/xrealft.cpp	Sun Nov 25 15:27:44 2001
--- ../cpp/examples/xrealft.cpp	Sun Oct  6 17:43:32 2002
***************
*** 41,43 ****
            cout << "press RETURN to continue ..." << endl;
!           cin.get();;
            NR::realft(data,-1);
--- 41,43 ----
            cout << "press RETURN to continue ..." << endl;
!           cin.get();
            NR::realft(data,-1);
diff -r -b -B -c1 cpp/examples/xsphbes.cpp ../cpp/examples/xsphbes.cpp
*** cpp/examples/xsphbes.cpp	Sun Nov 25 15:27:44 2001
--- ../cpp/examples/xsphbes.cpp	Sun Oct  6 17:43:38 2002
***************
*** 7,9 ****
  
! // Driver for routine rj
  
--- 7,9 ----
  
! // Driver for routine sphbes
  
diff -r -b -B -c1 cpp/other/nrtypes_lib.h ../cpp/other/nrtypes_lib.h
*** cpp/other/nrtypes_lib.h	Sun Nov 25 15:29:56 2001
--- ../cpp/other/nrtypes_lib.h	Sun Oct  6 17:39:28 2002
***************
*** 65,67 ****
  
! typedef const NRVec<unsigned long> Mat_I_ULNG;
  typedef const NRMat<unsigned long> Mat_ULNG, Mat_O_ULNG, Mat_IO_ULNG;
--- 65,67 ----
  
! typedef const NRMat<unsigned long> Mat_I_ULNG;
  typedef const NRMat<unsigned long> Mat_ULNG, Mat_O_ULNG, Mat_IO_ULNG;
diff -r -b -B -c1 cpp/other/nrtypes_nr.h ../cpp/other/nrtypes_nr.h
*** cpp/other/nrtypes_nr.h	Sun Nov 25 15:29:56 2001
--- ../cpp/other/nrtypes_nr.h	Sun Oct  6 17:39:38 2002
***************
*** 65,67 ****
  
! typedef const NRVec<unsigned long> Mat_I_ULNG;
  typedef NRMat<unsigned long> Mat_ULNG, Mat_O_ULNG, Mat_IO_ULNG;
--- 65,67 ----
  
! typedef const NRMat<unsigned long> Mat_I_ULNG;
  typedef NRMat<unsigned long> Mat_ULNG, Mat_O_ULNG, Mat_IO_ULNG;
diff -r -b -B -c1 cpp/recipes/eulsum.cpp ../cpp/recipes/eulsum.cpp
*** cpp/recipes/eulsum.cpp	Sun Nov 25 15:26:46 2001
--- ../cpp/recipes/eulsum.cpp	Sun Oct  6 17:31:42 2002
***************
*** 14,16 ****
  	} else {
! 		if (nterm+1 > wksp.size()) nrerror("wksp too small in euler");
  		tmp=wksp[0];
--- 14,16 ----
  	} else {
! 		if (nterm+1 > wksp.size()) nrerror("wksp too small in eulsum");
  		tmp=wksp[0];
diff -r -b -B -c1 cpp/recipes/hpsort.cpp ../cpp/recipes/hpsort.cpp
*** cpp/recipes/hpsort.cpp	Sun Nov 25 15:26:48 2001
--- ../cpp/recipes/hpsort.cpp	Sun Oct  6 17:31:58 2002
***************
*** 10,12 ****
  		jold=l;
! 		j=l+1;
  		while (j <= r) {
--- 10,12 ----
  		jold=l;
! 		j=2*l+1;
  		while (j <= r) {
diff -r -b -B -c1 cpp/recipes/mnbrak.cpp ../cpp/recipes/mnbrak.cpp
*** cpp/recipes/mnbrak.cpp	Sun Nov 25 15:26:48 2001
--- ../cpp/recipes/mnbrak.cpp	Sun Oct  6 17:32:04 2002
***************
*** 51,53 ****
  			if (fu < fc) {
! 				shft3(bx,cx,u,cx+GOLD*(cx-bx));
  				shft3(fb,fc,fu,func(u));
--- 51,53 ----
  			if (fu < fc) {
! 				shft3(bx,cx,u,u+GOLD*(u-cx));
  				shft3(fb,fc,fu,func(u));
diff -r -b -B -c1 cpp/recipes/mrqmin.cpp ../cpp/recipes/mrqmin.cpp
*** cpp/recipes/mrqmin.cpp	Sun Nov 25 15:26:48 2001
--- ../cpp/recipes/mrqmin.cpp	Sun Oct  6 17:31:32 2002
***************
*** 11,15 ****
  	int ma=a.size();
! 	static Mat_DP oneda(ma,1);
! 	static Vec_DP atry(ma),beta(ma),da(ma);
  	if (alamda < 0.0) {
  		mfit=0;
--- 11,18 ----
  	int ma=a.size();
! 	static Mat_DP *oneda_p;
! 	static Vec_DP *atry_p,*beta_p,*da_p;
  	if (alamda < 0.0) {
+ 		atry_p = new Vec_DP(ma);
+ 		beta_p = new Vec_DP(ma);
+ 		da_p = new Vec_DP(ma);
  		mfit=0;
***************
*** 17,23 ****
  			if (ia[j]) mfit++;
  		alamda=0.001;
! 		mrqcof(x,y,sig,a,ia,alpha,beta,chisq,funcs);
  		ochisq=chisq;
! 		for (j=0;j<ma;j++) atry[j]=a[j];
  	}
  	Mat_DP temp(mfit,mfit);
--- 20,29 ----
  			if (ia[j]) mfit++;
+ 		oneda_p = new Mat_DP(mfit,1);
  		alamda=0.001;
! 		mrqcof(x,y,sig,a,ia,alpha,*beta_p,chisq,funcs);
  		ochisq=chisq;
! 		for (j=0;j<ma;j++) (*atry_p)[j]=a[j];
  	}
+ 	Mat_DP &oneda=*oneda_p;
+ 	Vec_DP &atry=*atry_p,&beta=*beta_p,&da=*da_p;
  	Mat_DP temp(mfit,mfit);
***************
*** 37,38 ****
--- 43,45 ----
  		covsrt(alpha,ia,mfit);
+ 		delete oneda_p; delete da_p; delete beta_p; delete atry_p;
  		return;
diff -r -b -B -c1 cpp/recipes/simplx.cpp ../cpp/recipes/simplx.cpp
*** cpp/recipes/simplx.cpp	Sun Nov 25 15:26:50 2001
--- ../cpp/recipes/simplx.cpp	Sun Oct  6 17:32:10 2002
***************
*** 61,63 ****
  				kh=iposv[ip-1]-m1-n+1;
! 				if (kh >= 1 && l3[kh-1]) {
  					l3[kh-1]=0;
--- 61,63 ----
  				kh=iposv[ip-1]-m1-n+1;
! 				if (kh >= 1 && l3[kh-1] != 0) {
  					l3[kh-1]=0;
diff -r -b -B -c1 cpp/recipes/sprstp.cpp ../cpp/recipes/sprstp.cpp
*** cpp/recipes/sprstp.cpp	Sun Nov 25 15:26:50 2001
--- ../cpp/recipes/sprstp.cpp	Sun Oct  6 17:31:52 2002
***************
*** 9,12 ****
  	for (j=0;j<n2-1;j++) sb[j]=sa[j];
! 	Vec_INT ija_vec((int *) &ija[n2],ija[n2-1]-ija[0]);
! 	Vec_INT ijb_vec(&ijb[n2],ija[n2-1]-ija[0]);
  	indexx(ija_vec,ijb_vec);
--- 9,12 ----
  	for (j=0;j<n2-1;j++) sb[j]=sa[j];
! 	Vec_INT ija_vec(&ija[n2],ija[n2-1]-ija[0]);
! 	Vec_INT ijb_vec(ija[n2-1]-ija[0]);
  	indexx(ija_vec,ijb_vec);
diff -r -b -B -c1 cpp/recipes/svdcmp.cpp ../cpp/recipes/svdcmp.cpp
*** cpp/recipes/svdcmp.cpp	Sun Nov 25 15:26:50 2001
--- ../cpp/recipes/svdcmp.cpp	Sun Oct  6 17:31:38 2002
***************
*** 39,41 ****
  		g=s=scale=0.0;
! 		if (i+1 <= m && i != n) {
  			for (k=l-1;k<n;k++) scale += fabs(a[i][k]);
--- 39,41 ----
  		g=s=scale=0.0;
! 		if (i+1 <= m && i+1 != n) {
  			for (k=l-1;k<n;k++) scale += fabs(a[i][k]);
***************
*** 105,107 ****
  				s=1.0;
! 				for (i=l-1;i<k+1;i++) {
  					f=s*rv1[i];
--- 105,107 ----
  				s=1.0;
! 				for (i=l;i<k+1;i++) {
  					f=s*rv1[i];
diff -r -b -B -c1 cpp/recipes/wwghts.cpp ../cpp/recipes/wwghts.cpp
*** cpp/recipes/wwghts.cpp	Sun Nov 25 15:26:50 2001
--- ../cpp/recipes/wwghts.cpp	Sun Oct  6 17:32:14 2002
***************
*** 6,8 ****
  	int j,k;
- 	Vec_DP wold(4);
  	DP hh,hi,c,fac,a,b;
--- 6,7 ----