version 1.1.1.1, 2000/01/09 17:00:49 |
version 1.1.1.3, 2003/09/15 07:09:23 |
Line 212 register int nl, nh; |
|
Line 212 register int nl, nh; |
|
{ |
{ |
register float GPFAR *vec; |
register float GPFAR *vec; |
|
|
if (!(vec = (float GPFAR *) gp_alloc((unsigned long) (nh - nl + 1) * sizeof(float), NULL))) { |
if (!(vec = (float GPFAR *) gp_alloc((nh - nl + 1) * sizeof(float), NULL))) { |
int_error("not enough memory to create vector", NO_CARET); |
int_error("not enough memory to create vector", NO_CARET); |
return NULL; /* Not reached */ |
return NULL; /* Not reached */ |
} |
} |
Line 286 register int nrl, nrh, ncl, nch; |
|
Line 286 register int nrl, nrh, ncl, nch; |
|
register int i; |
register int i; |
register float GPFAR *GPFAR * m; |
register float GPFAR *GPFAR * m; |
|
|
m = (float GPFAR * GPFAR *) gp_alloc((unsigned long) (nrh - nrl + 1) * sizeof(float GPFAR *), "matrix"); |
m = (float GPFAR * GPFAR *) gp_alloc((nrh - nrl + 1) * sizeof(float GPFAR *), "matrix"); |
m -= nrl; |
m -= nrl; |
|
|
for (i = nrl; i <= nrh; i++) { |
for (i = nrl; i <= nrh; i++) { |
if (!(m[i] = (float GPFAR *) gp_alloc((unsigned long) (nch - ncl + 1) * sizeof(float), NULL))) { |
if (!(m[i] = (float GPFAR *) gp_alloc((nch - ncl + 1) * sizeof(float), NULL))) { |
free_matrix(m, nrl, i - 1, ncl, nch); |
free_matrix(m, nrl, i - 1, ncl, nch); |
int_error("not enough memory to create matrix", NO_CARET); |
int_error("not enough memory to create matrix", NO_CARET); |
return NULL; |
return NULL; |
Line 344 register int srh, sch; |
|
Line 344 register int srh, sch; |
|
} |
} |
} |
} |
for (i = nrh + 1; i <= srh; i++) { |
for (i = nrh + 1; i <= srh; i++) { |
if (!(m[i] = (float GPFAR *) gp_alloc((unsigned long) (nch - ncl + 1) * sizeof(float), NULL))) { |
if (!(m[i] = (float GPFAR *) gp_alloc((nch - ncl + 1) * sizeof(float), NULL))) { |
free_matrix(m, nrl, i - 1, nrl, sch); |
free_matrix(m, nrl, i - 1, nrl, sch); |
int_error("not enough memory to extend matrix", NO_CARET); |
int_error("not enough memory to extend matrix", NO_CARET); |
return NULL; |
return NULL; |
Line 402 register int nrl, nrh, ncl, nch; |
|
Line 402 register int nrl, nrh, ncl, nch; |
|
|
|
nrow = nrh - nrl + 1; |
nrow = nrh - nrl + 1; |
ncol = nch - ncl + 1; |
ncol = nch - ncl + 1; |
m = (float GPFAR * GPFAR *) gp_alloc((unsigned long) (nrh - nrl + 1) * sizeof(float GPFAR *), "convert_matrix"); |
m = (float GPFAR * GPFAR *) gp_alloc((nrh - nrl + 1) * sizeof(float GPFAR *), "convert_matrix"); |
m -= nrl; |
m -= nrl; |
|
|
m[nrl] = a - ncl; |
m[nrl] = a - ncl; |