Changeset 2215 for trunk/psLib/src/math/psMinimize.c
- Timestamp:
- Oct 27, 2004, 10:29:12 AM (22 years ago)
- File:
-
- 1 edited
-
trunk/psLib/src/math/psMinimize.c (modified) (5 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/math/psMinimize.c
r2214 r2215 9 9 * @author GLG, MHPCC 10 10 * 11 * @version $Revision: 1.6 8$ $Name: not supported by cvs2svn $12 * @date $Date: 2004-10-27 20:2 0:11$11 * @version $Revision: 1.69 $ $Name: not supported by cvs2svn $ 12 * @date $Date: 2004-10-27 20:29:12 $ 13 13 * 14 14 * Copyright 2004 Maui High Performance Computing Center, University of Hawaii … … 255 255 PS_PTR_CHECK_NULL(y, NULL); 256 256 PS_VECTOR_CHECK_TYPE_F32_OR_F64(y, NULL); 257 CONVERT_VECTOR_F64_TO_F32(y, y32, y32Static);257 PS_VECTOR_CONVERT_F64_TO_F32_STATIC(y, y32, y32Static); 258 258 259 259 // If yErr==NULL, set all errors equal. 260 260 if (yErr == NULL) { 261 GEN_PERSIST_YERR_VEC_F32(yErr32Static, y->n);261 PS_VECTOR_GEN_YERR_STATIC_F32(yErr32Static, y->n); 262 262 yErr32 = yErr32Static; 263 263 } else { 264 264 PS_VECTOR_CHECK_TYPE_F32_OR_F64(yErr, NULL); 265 CONVERT_VECTOR_F64_TO_F32(yErr, yErr32, yErr32Static);265 PS_VECTOR_CONVERT_F64_TO_F32_STATIC(yErr, yErr32, yErr32Static); 266 266 } 267 267 268 268 // If x==NULL, create an x32 vector with x values set to (0:n). 269 269 if (x == NULL) { 270 GEN_PERSIST_X_INDEX_VEC_F32(x32Static, y->n);270 PS_VECTOR_GEN_X_INDEX_STATIC_F32(x32Static, y->n); 271 271 x32 = x32Static; 272 272 } else { 273 273 PS_VECTOR_CHECK_TYPE_F32_OR_F64(x, NULL); 274 CONVERT_VECTOR_F64_TO_F32(x, x32, x32Static);274 PS_VECTOR_CONVERT_F64_TO_F32_STATIC(x, x32, x32Static); 275 275 } 276 276 PS_VECTOR_CHECK_SIZE_EQUAL(x32, y32, NULL); … … 668 668 double fac; 669 669 double sum; 670 GEN_STATIC_RECYCLED_VECTOR(f, n, PS_TYPE_F64);670 PS_VECTOR_GEN_STATIC_RECYCLED(f, n, PS_TYPE_F64); 671 671 psScalar *fScalar; 672 672 psScalar tmpScalar; … … 866 866 867 867 PS_VECTOR_CHECK_TYPE_F32_OR_F64(y, NULL); 868 CONVERT_VECTOR_F32_TO_F64(y, y64, y64Static);868 PS_VECTOR_CONVERT_F32_TO_F64_STATIC(y, y64, y64Static); 869 869 870 870 // If yErr==NULL, set all errors equal. 871 871 if (yErr == NULL) { 872 GEN_PERSIST_YERR_VEC_F64(yErr64Static, y->n);872 PS_VECTOR_GEN_YERR_STATIC_F64(yErr64Static, y->n); 873 873 yErr64 = yErr64Static; 874 874 } else { 875 875 PS_VECTOR_CHECK_TYPE_F32_OR_F64(yErr, NULL); 876 CONVERT_VECTOR_F32_TO_F64(yErr, yErr64, yErr64Static);876 PS_VECTOR_CONVERT_F32_TO_F64_STATIC(yErr, yErr64, yErr64Static); 877 877 } 878 878 879 879 // If x==NULL, create an x64 vector with x values set to (0:n). 880 880 if (x == NULL) { 881 GEN_PERSIST_X_INDEX_VEC_F64(x64Static, y->n);881 PS_VECTOR_GEN_X_INDEX_STATIC_F64(x64Static, y->n); 882 882 if (myPoly->type == PS_POLYNOMIAL_CHEB) { 883 883 p_psNormalizeVector(x64Static); … … 886 886 } else { 887 887 PS_VECTOR_CHECK_TYPE_F32_OR_F64(x, NULL); 888 CONVERT_VECTOR_F32_TO_F64(x, x64, x64Static);888 PS_VECTOR_CONVERT_F32_TO_F64_STATIC(x, x64, x64Static); 889 889 } 890 890 PS_VECTOR_CHECK_SIZE_EQUAL(x64, y64, NULL);
Note:
See TracChangeset
for help on using the changeset viewer.
