5 void nrerror(error_text)
10 fprintf(stderr,"Numerical Recipes run-time error...\n");
11 fprintf(stderr,"%s\n",error_text);
12 fprintf(stderr,"...now exiting to system...\n");
23 v=(double *)malloc((unsigned) (nh-nl+1)*sizeof(double));
24 if (!v) nrerror("allocation failure in vector()");
33 v=(int *)malloc((unsigned) (nh-nl+1)*sizeof(int));
34 if (!v) nrerror("allocation failure in ivector()");
38 double *dvector(nl,nh)
43 v=(double *)malloc((unsigned) (nh-nl+1)*sizeof(double));
44 if (!v) nrerror("allocation failure in dvector()");
50 double **matrix(nrl,nrh,ncl,nch)
56 m=(double **) malloc((unsigned) (nrh-nrl+1)*sizeof(double*));
57 if (!m) nrerror("allocation failure 1 in matrix()");
60 for(i=nrl;i<=nrh;i++) {
61 m[i]=(double *) malloc((unsigned) (nch-ncl+1)*sizeof(double));
62 if (!m[i]) nrerror("allocation failure 2 in matrix()");
68 double **dmatrix(nrl,nrh,ncl,nch)
74 m=(double **) malloc((unsigned) (nrh-nrl+1)*sizeof(double*));
75 if (!m) nrerror("allocation failure 1 in dmatrix()");
78 for(i=nrl;i<=nrh;i++) {
79 m[i]=(double *) malloc((unsigned) (nch-ncl+1)*sizeof(double));
80 if (!m[i]) nrerror("allocation failure 2 in dmatrix()");
86 int **imatrix(nrl,nrh,ncl,nch)
91 m=(int **)malloc((unsigned) (nrh-nrl+1)*sizeof(int*));
92 if (!m) nrerror("allocation failure 1 in imatrix()");
95 for(i=nrl;i<=nrh;i++) {
96 m[i]=(int *)malloc((unsigned) (nch-ncl+1)*sizeof(int));
97 if (!m[i]) nrerror("allocation failure 2 in imatrix()");
105 double **submatrix(a,oldrl,oldrh,oldcl,oldch,newrl,newcl)
107 int oldrl,oldrh,oldcl,oldch,newrl,newcl;
112 m=(double **) malloc((unsigned) (oldrh-oldrl+1)*sizeof(double*));
113 if (!m) nrerror("allocation failure in submatrix()");
116 for(i=oldrl,j=newrl;i<=oldrh;i++,j++) m[j]=a[i]+oldcl-newcl;
123 void free_vector(v,nl,nh)
127 free((char*) (v+nl));
130 void free_ivector(v,nl,nh)
133 free((char*) (v+nl));
136 void free_dvector(v,nl,nh)
140 free((char*) (v+nl));
145 void free_matrix(m,nrl,nrh,ncl,nch)
151 for(i=nrh;i>=nrl;i--) free((char*) (m[i]+ncl));
152 free((char*) (m+nrl));
155 void free_dmatrix(m,nrl,nrh,ncl,nch)
161 for(i=nrh;i>=nrl;i--) free((char*) (m[i]+ncl));
162 free((char*) (m+nrl));
165 void free_imatrix(m,nrl,nrh,ncl,nch)
171 for(i=nrh;i>=nrl;i--) free((char*) (m[i]+ncl));
172 free((char*) (m+nrl));
177 void free_submatrix(b,nrl,nrh,ncl,nch)
181 free((char*) (b+nrl));
186 double **convert_matrix(a,nrl,nrh,ncl,nch)
195 m = (double **) malloc((unsigned) (nrow)*sizeof(double*));
196 if (!m) nrerror("allocation failure in convert_matrix()");
198 for(i=0,j=nrl;i<=nrow-1;i++,j++) m[j]=a+ncol*i-ncl;
204 void free_convert_matrix(b,nrl,nrh,ncl,nch)
208 free((char*) (b+nrl));