Index: trunk/psLib/test/math/tst_psFunc00.c
===================================================================
--- trunk/psLib/test/math/tst_psFunc00.c	(revision 5070)
+++ trunk/psLib/test/math/tst_psFunc00.c	(revision 5090)
@@ -14,6 +14,6 @@
 *    orders are created.
 * 
-*    @version $Revision: 1.3 $  $Name: not supported by cvs2svn $
-*    @date $Date: 2005-08-08 21:44:09 $
+*    @version $Revision: 1.4 $  $Name: not supported by cvs2svn $
+*    @date $Date: 2005-09-22 02:47:16 $
 *
 *  Copyright 2004-2005 Maui High Performance Computing Center, Univ. of Hawaii
@@ -63,7 +63,7 @@
     }
     // Verify polynomial structure members set properly
-    if(my1DPoly->n != ORDER) {
-        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
-                my1DPoly->n, ORDER);
+    if(my1DPoly->COOL_1D_n != ORDER) {
+        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
+                my1DPoly->COOL_1D_n, ORDER);
         return 2;
     }
@@ -73,5 +73,5 @@
         return 3;
     }
-    for(psS32 i = 0; i < ORDER; i++) {
+    for(psS32 i = 0; i < ORDER+1; i++) {
         if(my1DPoly->coeff[i] != 0.0) {
             psError(PS_ERR_UNKNOWN,true,"Coeff[%d] %lg not as expected %lg",
@@ -115,13 +115,13 @@
     }
     // Verify polynomial structure members set properly
-    if(my2DPoly->nX != ORDER) {
-        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
-                my2DPoly->nX, ORDER);
+    if(my2DPoly->COOL_2D_nX != ORDER) {
+        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
+                my2DPoly->COOL_2D_nX, ORDER);
         return 2;
     }
     // Verify polynomial structure members set properly
-    if(my2DPoly->nY != ORDER+1) {
-        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
-                my2DPoly->nY, ORDER+1);
+    if(my2DPoly->COOL_2D_nY != ORDER+1) {
+        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
+                my2DPoly->COOL_2D_nY, ORDER+1);
         return 3;
     }
@@ -131,6 +131,6 @@
         return 4;
     }
-    for(psS32 i = 0; i < ORDER; i++) {
-        for(psS32 j = 0; j < ORDER+1; j++) {
+    for(psS32 i = 0; i < ORDER+1; i++) {
+        for(psS32 j = 0; j < ORDER+2; j++) {
             if(my2DPoly->coeff[i][j] != 0.0) {
                 psError(PS_ERR_UNKNOWN,true,"Coeff[%d][%d] %lg not as expected %lg",
@@ -181,19 +181,19 @@
     }
     // Verify polynomial structure members set properly
-    if(my3DPoly->nX != ORDER) {
-        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
-                my3DPoly->nX, ORDER);
+    if(my3DPoly->COOL_3D_nX != ORDER) {
+        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
+                my3DPoly->COOL_3D_nX, ORDER);
         return 2;
     }
     // Verify polynomial structure members set properly
-    if(my3DPoly->nY != ORDER+1) {
-        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
-                my3DPoly->nY, ORDER+1);
+    if(my3DPoly->COOL_3D_nY != ORDER+1) {
+        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
+                my3DPoly->COOL_3D_nY, ORDER+1);
         return 3;
     }
     // Verify polynomial structure members set properly
-    if(my3DPoly->nZ != ORDER+2) {
-        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
-                my3DPoly->nZ, ORDER+2);
+    if(my3DPoly->COOL_3D_nZ != ORDER+2) {
+        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
+                my3DPoly->COOL_3D_nZ, ORDER+2);
         return 4;
     }
@@ -203,7 +203,7 @@
         return 5;
     }
-    for(psS32 i = 0; i < ORDER; i++) {
-        for(psS32 j = 0; j < ORDER+1; j++) {
-            for(psS32 k = 0; k < ORDER+2; k++) {
+    for(psS32 i = 0; i < ORDER+1; i++) {
+        for(psS32 j = 0; j < ORDER+2; j++) {
+            for(psS32 k = 0; k < ORDER+3; k++) {
                 if(my3DPoly->coeff[i][j][k] != 0.0) {
                     psError(PS_ERR_UNKNOWN,true,"Coeff[%d][%d][%d] %lg not as expected %lg",
@@ -254,5 +254,5 @@
 
     // Allocate polynomial
-    my4DPoly = psPolynomial4DAlloc(ORDER+3,ORDER,ORDER+1,ORDER+2,PS_POLYNOMIAL_ORD);
+    my4DPoly = psPolynomial4DAlloc(ORDER,ORDER+1,ORDER+2,ORDER+3,PS_POLYNOMIAL_ORD);
     // Verify structure allocated
     if(my4DPoly == NULL) {
@@ -261,25 +261,25 @@
     }
     // Verify polynomial structure members set properly
-    if(my4DPoly->nY != ORDER) {
-        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
-                my4DPoly->nY, ORDER);
+    if(my4DPoly->COOL_4D_nX != ORDER) {
+        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
+                my4DPoly->COOL_4D_nX, ORDER);
         return 2;
     }
     // Verify polynomial structure members set properly
-    if(my4DPoly->nZ != ORDER+1) {
-        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
-                my4DPoly->nZ, ORDER+1);
+    if(my4DPoly->COOL_4D_nY != ORDER+1) {
+        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
+                my4DPoly->COOL_4D_nX, ORDER+1);
         return 3;
     }
     // Verify polynomial structure members set properly
-    if(my4DPoly->nT != ORDER+2) {
-        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
-                my4DPoly->nT, ORDER+2);
+    if(my4DPoly->COOL_4D_nZ != ORDER+2) {
+        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
+                my4DPoly->COOL_4D_nZ, ORDER+2);
         return 4;
     }
     // Verify polynomial structure members set properly
-    if(my4DPoly->nX != ORDER+3) {
-        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
-                my4DPoly->nX, ORDER+3);
+    if(my4DPoly->COOL_4D_nT != ORDER+3) {
+        psError(PS_ERR_UNKNOWN,true,"Number of terms %d not as expected %d",
+                my4DPoly->COOL_4D_nT, ORDER+3);
         return 5;
     }
@@ -289,8 +289,8 @@
         return 6;
     }
-    for(psS32 i = 0; i < ORDER+3; i++) {
-        for(psS32 j = 0; j < ORDER; j++) {
-            for(psS32 k = 0; k < ORDER+1; k++) {
-                for(psS32 l = 0; l < ORDER+2; l++) {
+    for(psS32 i = 0; i < ORDER+1; i++) {
+        for(psS32 j = 0; j < ORDER+2; j++) {
+            for(psS32 k = 0; k < ORDER+3; k++) {
+                for(psS32 l = 0; l < ORDER+4; l++) {
                     if(my4DPoly->coeff[i][j][k][l] != 0.0) {
                         psError(PS_ERR_UNKNOWN,true,"Coeff[%d][%d][%d][%d] %lg not as expected %lg",
Index: trunk/psLib/test/math/tst_psFunc08.c
===================================================================
--- trunk/psLib/test/math/tst_psFunc08.c	(revision 5070)
+++ trunk/psLib/test/math/tst_psFunc08.c	(revision 5090)
@@ -4,6 +4,6 @@
 *  ORD and CHEB type polynomials.
 *
-*  @version  $Revision: 1.3 $  $Name: not supported by cvs2svn $
-*  @date  $Date: 2005-08-08 21:44:09 $
+*  @version  $Revision: 1.4 $  $Name: not supported by cvs2svn $
+*  @date  $Date: 2005-09-22 02:47:16 $
 *
 *  XXX: Probably should test single- and multi-dimensional polynomials in
@@ -59,6 +59,6 @@
 
     // Allocate polynomial structure
-    psPolynomial1D*  polyOrd = psPolynomial1DAlloc(TERMS, PS_POLYNOMIAL_ORD);
-    psPolynomial1D*  polyCheb = psPolynomial1DAlloc(TERMS, PS_POLYNOMIAL_CHEB);
+    psPolynomial1D*  polyOrd = psPolynomial1DAlloc(TERMS-1, PS_POLYNOMIAL_ORD);
+    psPolynomial1D*  polyCheb = psPolynomial1DAlloc(TERMS-1, PS_POLYNOMIAL_CHEB);
     // Set polynomial members
     for(psS32 i = 0; i < TERMS; i++) {
@@ -87,5 +87,5 @@
 
     // Allocate polynomial with invalid type
-    polyOrd = psPolynomial1DAlloc(TERMS, 99);
+    polyOrd = psPolynomial1DAlloc(TERMS-1, 99);
     // Attempt to evaluation invalid polynomial type
     psLogMsg(__func__,PS_LOG_INFO,"Following should generate error message invalid type");
@@ -104,6 +104,6 @@
 {
     // Allocate polynomial
-    psPolynomial1D* polyOrd = psPolynomial1DAlloc(TERMS,PS_POLYNOMIAL_ORD);
-    psPolynomial1D* polyCheb = psPolynomial1DAlloc(TERMS,PS_POLYNOMIAL_CHEB);
+    psPolynomial1D* polyOrd = psPolynomial1DAlloc(TERMS-1,PS_POLYNOMIAL_ORD);
+    psPolynomial1D* polyCheb = psPolynomial1DAlloc(TERMS-1,PS_POLYNOMIAL_CHEB);
 
     // Set polynomial members
Index: trunk/psLib/test/math/tst_psFunc09.c
===================================================================
--- trunk/psLib/test/math/tst_psFunc09.c	(revision 5070)
+++ trunk/psLib/test/math/tst_psFunc09.c	(revision 5090)
@@ -4,6 +4,6 @@
 *  ORD and CHEB type polynomials.
 *
-*  @version  $Revision: 1.3 $  $Name: not supported by cvs2svn $
-*  @date  $Date: 2005-08-08 21:44:09 $
+*  @version  $Revision: 1.4 $  $Name: not supported by cvs2svn $
+*  @date  $Date: 2005-09-22 02:47:16 $
 *
 * Copyright 2004-2005 Maui High Performance Computing Center, Univ. of Hawaii
@@ -88,6 +88,6 @@
  
     // Allocate polynomial structure
-    psPolynomial2D*  polyOrd = psPolynomial2DAlloc(TERMS, TERMS, PS_POLYNOMIAL_ORD);
-    psPolynomial2D*  polyCheb = psPolynomial2DAlloc(TERMS, TERMS, PS_POLYNOMIAL_CHEB);
+    psPolynomial2D*  polyOrd = psPolynomial2DAlloc(TERMS-1, TERMS-1, PS_POLYNOMIAL_ORD);
+    psPolynomial2D*  polyCheb = psPolynomial2DAlloc(TERMS-1, TERMS-1, PS_POLYNOMIAL_CHEB);
     // Set polynomial members
     for(psS32 i = 0; i < TERMS; i++) {
@@ -118,5 +118,5 @@
  
     // Allocate polynomial with invalid type
-    polyOrd = psPolynomial2DAlloc(TERMS, TERMS, 99);
+    polyOrd = psPolynomial2DAlloc(TERMS-1, TERMS-1, 99);
     // Attempt to evaluation invalid polynomial type
     psLogMsg(__func__,PS_LOG_INFO,"Following should generate error message invalid type");
@@ -138,6 +138,6 @@
 
     // Allocate polynomial structure
-    psPolynomial2D*  polyOrd = psPolynomial2DAlloc(TERMS, TERMS, PS_POLYNOMIAL_ORD);
-    psPolynomial2D*  polyCheb = psPolynomial2DAlloc(TERMS, TERMS, PS_POLYNOMIAL_CHEB);
+    psPolynomial2D*  polyOrd = psPolynomial2DAlloc(TERMS-1, TERMS-1, PS_POLYNOMIAL_ORD);
+    psPolynomial2D*  polyCheb = psPolynomial2DAlloc(TERMS-1, TERMS-1, PS_POLYNOMIAL_CHEB);
     // Set polynomial members
     for(psS32 i = 0; i < TERMS; i++) {
@@ -168,5 +168,5 @@
 
     // Allocate polynomial with invalid type
-    polyOrd = psPolynomial2DAlloc(TERMS, TERMS, 99);
+    polyOrd = psPolynomial2DAlloc(TERMS-1, TERMS-1, 99);
     // Attempt to evaluation invalid polynomial type
     psLogMsg(__func__,PS_LOG_INFO,"Following should generate error message invalid type");
@@ -184,6 +184,6 @@
 {
     // Allocate polynomial
-    psPolynomial2D* polyOrd = psPolynomial2DAlloc(TERMS,TERMS,PS_POLYNOMIAL_ORD);
-    psPolynomial2D* polyCheb = psPolynomial2DAlloc(TERMS,TERMS,PS_POLYNOMIAL_CHEB);
+    psPolynomial2D* polyOrd = psPolynomial2DAlloc(TERMS-1,TERMS-1,PS_POLYNOMIAL_ORD);
+    psPolynomial2D* polyCheb = psPolynomial2DAlloc(TERMS-1,TERMS-1,PS_POLYNOMIAL_CHEB);
  
     // Set polynomial members
@@ -297,6 +297,6 @@
 {
     // Allocate polynomial
-    psPolynomial2D* polyOrd = psPolynomial2DAlloc(TERMS,TERMS,PS_POLYNOMIAL_ORD);
-    psPolynomial2D* polyCheb = psPolynomial2DAlloc(TERMS,TERMS,PS_POLYNOMIAL_CHEB);
+    psPolynomial2D* polyOrd = psPolynomial2DAlloc(TERMS-1,TERMS-1,PS_POLYNOMIAL_ORD);
+    psPolynomial2D* polyCheb = psPolynomial2DAlloc(TERMS-1,TERMS-1,PS_POLYNOMIAL_CHEB);
 
     // Set polynomial members
Index: trunk/psLib/test/math/tst_psFunc10.c
===================================================================
--- trunk/psLib/test/math/tst_psFunc10.c	(revision 5070)
+++ trunk/psLib/test/math/tst_psFunc10.c	(revision 5090)
@@ -4,6 +4,6 @@
 *  ORD and CHEB type polynomials.
 *
-*  @version  $Revision: 1.3 $  $Name: not supported by cvs2svn $
-*  @date  $Date: 2005-08-08 21:44:09 $
+*  @version  $Revision: 1.4 $  $Name: not supported by cvs2svn $
+*  @date  $Date: 2005-09-22 02:47:16 $
 *
 * Copyright 2004-2005 Maui High Performance Computing Center, Univ. of Hawaii
@@ -139,6 +139,6 @@
  
     // Allocate polynomial structure
-    psPolynomial3D*  polyOrd = psPolynomial3DAlloc(TERMS, TERMS, TERMS, PS_POLYNOMIAL_ORD);
-    psPolynomial3D*  polyCheb = psPolynomial3DAlloc(TERMS, TERMS, TERMS, PS_POLYNOMIAL_CHEB);
+    psPolynomial3D*  polyOrd = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_ORD);
+    psPolynomial3D*  polyCheb = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_CHEB);
     // Set polynomial members
     for(psS32 i = 0; i < TERMS; i++) {
@@ -173,5 +173,5 @@
  
     // Allocate polynomial with invalid type
-    polyOrd = psPolynomial3DAlloc(TERMS, TERMS, TERMS, 99);
+    polyOrd = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1, 99);
     // Attempt to evaluation invalid polynomial type
     psLogMsg(__func__,PS_LOG_INFO,"Following should generate error message invalid type");
@@ -193,6 +193,6 @@
 
     // Allocate polynomial structure
-    psPolynomial3D*  polyOrd = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_ORD);
-    psPolynomial3D*  polyCheb = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_CHEB);
+    psPolynomial3D*  polyOrd = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_ORD);
+    psPolynomial3D*  polyCheb = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_CHEB);
     // Set polynomial members
     for(psS32 i = 0; i < TERMS; i++) {
@@ -227,5 +227,5 @@
 
     // Allocate polynomial with invalid type
-    polyOrd = psPolynomial3DAlloc(TERMS, TERMS, TERMS, 99);
+    polyOrd = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1, 99);
     // Attempt to evaluation invalid polynomial type
     psLogMsg(__func__,PS_LOG_INFO,"Following should generate error message invalid type");
@@ -243,6 +243,6 @@
 {
     // Allocate polynomial
-    psPolynomial3D* polyOrd  = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_ORD);
-    psPolynomial3D* polyCheb = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_CHEB);
+    psPolynomial3D* polyOrd  = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_ORD);
+    psPolynomial3D* polyCheb = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_CHEB);
  
     // Set polynomial members
@@ -378,6 +378,6 @@
 {
     // Allocate polynomial
-    psPolynomial3D* polyOrd = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_ORD);
-    psPolynomial3D* polyCheb = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_CHEB);
+    psPolynomial3D* polyOrd = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_ORD);
+    psPolynomial3D* polyCheb = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_CHEB);
 
     // Set polynomial members
Index: trunk/psLib/test/math/tst_psFunc11.c
===================================================================
--- trunk/psLib/test/math/tst_psFunc11.c	(revision 5070)
+++ trunk/psLib/test/math/tst_psFunc11.c	(revision 5090)
@@ -4,6 +4,6 @@
 *  ORD and CHEB type polynomials.
 *
-*  @version  $Revision: 1.3 $  $Name: not supported by cvs2svn $
-*  @date  $Date: 2005-08-08 21:44:09 $
+*  @version  $Revision: 1.4 $  $Name: not supported by cvs2svn $
+*  @date  $Date: 2005-09-22 02:47:16 $
 *
 * Copyright 2004-2005 Maui High Performance Computing Center, Univ. of Hawaii
@@ -381,6 +381,6 @@
  
     // Allocate polynomial structure
-    psPolynomial4D*  polyOrd  = psPolynomial4DAlloc(TERMS,TERMS,TERMS,TERMS,PS_POLYNOMIAL_ORD);
-    psPolynomial4D*  polyCheb = psPolynomial4DAlloc(TERMS,TERMS,TERMS,TERMS,PS_POLYNOMIAL_CHEB);
+    psPolynomial4D*  polyOrd  = psPolynomial4DAlloc(TERMS-1, TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_ORD);
+    psPolynomial4D*  polyCheb = psPolynomial4DAlloc(TERMS-1, TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_CHEB);
     // Set polynomial members
     for(psS32 i = 0; i < TERMS; i++) {
@@ -417,5 +417,5 @@
  
     // Allocate polynomial with invalid type
-    polyOrd = psPolynomial4DAlloc(TERMS, TERMS, TERMS, TERMS, 99);
+    polyOrd = psPolynomial4DAlloc(TERMS-1, TERMS-1, TERMS-1, TERMS-1, 99);
     // Attempt to evaluation invalid polynomial type
     psLogMsg(__func__,PS_LOG_INFO,"Following should generate error message invalid type");
@@ -437,6 +437,6 @@
 
     // Allocate polynomial structure
-    psPolynomial4D*  polyOrd = psPolynomial4DAlloc(TERMS,TERMS,TERMS,TERMS,PS_POLYNOMIAL_ORD);
-    psPolynomial4D*  polyCheb = psPolynomial4DAlloc(TERMS,TERMS,TERMS,TERMS,PS_POLYNOMIAL_CHEB);
+    psPolynomial4D*  polyOrd = psPolynomial4DAlloc(TERMS-1, TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_ORD);
+    psPolynomial4D*  polyCheb = psPolynomial4DAlloc(TERMS-1, TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_CHEB);
     // Set polynomial members
     for(psS32 i = 0; i < TERMS; i++) {
@@ -473,5 +473,5 @@
 
     // Allocate polynomial with invalid type
-    polyOrd = psPolynomial4DAlloc(TERMS, TERMS, TERMS, TERMS, 99);
+    polyOrd = psPolynomial4DAlloc(TERMS-1, TERMS-1, TERMS-1, TERMS-1, 99);
     // Attempt to evaluation invalid polynomial type
     psLogMsg(__func__,PS_LOG_INFO,"Following should generate error message invalid type");
@@ -489,6 +489,6 @@
 {
     // Allocate polynomial
-    psPolynomial4D* polyOrd  = psPolynomial4DAlloc(TERMS,TERMS,TERMS,TERMS,PS_POLYNOMIAL_ORD);
-    psPolynomial4D* polyCheb = psPolynomial4DAlloc(TERMS,TERMS,TERMS,TERMS,PS_POLYNOMIAL_CHEB);
+    psPolynomial4D* polyOrd  = psPolynomial4DAlloc(TERMS-1, TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_ORD);
+    psPolynomial4D* polyCheb = psPolynomial4DAlloc(TERMS-1, TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_CHEB);
  
     // Set polynomial members
@@ -646,6 +646,6 @@
 {
     // Allocate polynomial
-    psPolynomial4D* polyOrd = psPolynomial4DAlloc(TERMS,TERMS,TERMS,TERMS,PS_POLYNOMIAL_ORD);
-    psPolynomial4D* polyCheb = psPolynomial4DAlloc(TERMS,TERMS,TERMS,TERMS,PS_POLYNOMIAL_CHEB);
+    psPolynomial4D* polyOrd = psPolynomial4DAlloc(TERMS-1, TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_ORD);
+    psPolynomial4D* polyCheb = psPolynomial4DAlloc(TERMS-1, TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_CHEB);
 
     // Set polynomial members
Index: trunk/psLib/test/math/tst_psMinimize04.c
===================================================================
--- trunk/psLib/test/math/tst_psMinimize04.c	(revision 5070)
+++ trunk/psLib/test/math/tst_psMinimize04.c	(revision 5090)
@@ -45,5 +45,5 @@
     double actualData;
 
-    psPolynomial2D *myPoly = psPolynomial2DAlloc(POLY_ORDER+1, POLY_ORDER+1, PS_POLYNOMIAL_ORD);
+    psPolynomial2D *myPoly = psPolynomial2DAlloc(POLY_ORDER, POLY_ORDER, PS_POLYNOMIAL_ORD);
     psVector *x = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
     psVector *y = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
Index: trunk/psLib/test/math/tst_psMinimize04_F32.c
===================================================================
--- trunk/psLib/test/math/tst_psMinimize04_F32.c	(revision 5070)
+++ trunk/psLib/test/math/tst_psMinimize04_F32.c	(revision 5090)
@@ -45,5 +45,5 @@
     float actualData;
 
-    myPoly = psPolynomial1DAlloc(POLY_ORDER+1, PS_POLYNOMIAL_ORD);
+    myPoly = psPolynomial1DAlloc(POLY_ORDER, PS_POLYNOMIAL_ORD);
     x = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
     y = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
@@ -114,5 +114,5 @@
     float actualData;
 
-    myPoly = psPolynomial1DAlloc(POLY_ORDER+1, PS_POLYNOMIAL_ORD);
+    myPoly = psPolynomial1DAlloc(POLY_ORDER, PS_POLYNOMIAL_ORD);
     x = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
     y = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
@@ -179,5 +179,5 @@
     float actualData;
 
-    myPoly = psPolynomial1DAlloc(POLY_ORDER+1, PS_POLYNOMIAL_ORD);
+    myPoly = psPolynomial1DAlloc(POLY_ORDER, PS_POLYNOMIAL_ORD);
     y = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
 
Index: trunk/psLib/test/math/tst_psMinimize04b.c
===================================================================
--- trunk/psLib/test/math/tst_psMinimize04b.c	(revision 5070)
+++ trunk/psLib/test/math/tst_psMinimize04b.c	(revision 5090)
@@ -38,5 +38,5 @@
     psS32 memLeaks = 0;
 
-    myPoly = psPolynomial1DAlloc(POLY_ORDER+1, PS_POLYNOMIAL_CHEB);
+    myPoly = psPolynomial1DAlloc(POLY_ORDER, PS_POLYNOMIAL_CHEB);
     x = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
     y = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
@@ -109,5 +109,5 @@
     psS32 memLeaks = 0;
 
-    myPoly = psPolynomial1DAlloc(POLY_ORDER+1, PS_POLYNOMIAL_CHEB);
+    myPoly = psPolynomial1DAlloc(POLY_ORDER, PS_POLYNOMIAL_CHEB);
     x = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
     y = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
@@ -172,5 +172,5 @@
     psS32 memLeaks = 0;
 
-    myPoly = psPolynomial1DAlloc(POLY_ORDER+1, PS_POLYNOMIAL_CHEB);
+    myPoly = psPolynomial1DAlloc(POLY_ORDER, PS_POLYNOMIAL_CHEB);
     x = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
     y = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
Index: trunk/psLib/test/math/tst_psMinimize04b_F32.c
===================================================================
--- trunk/psLib/test/math/tst_psMinimize04b_F32.c	(revision 5070)
+++ trunk/psLib/test/math/tst_psMinimize04b_F32.c	(revision 5090)
@@ -38,5 +38,5 @@
     psS32 memLeaks = 0;
 
-    myPoly = psPolynomial1DAlloc(POLY_ORDER+1, PS_POLYNOMIAL_CHEB);
+    myPoly = psPolynomial1DAlloc(POLY_ORDER, PS_POLYNOMIAL_CHEB);
     x = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
     y = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
@@ -109,5 +109,5 @@
     psS32 memLeaks = 0;
 
-    myPoly = psPolynomial1DAlloc(POLY_ORDER+1, PS_POLYNOMIAL_CHEB);
+    myPoly = psPolynomial1DAlloc(POLY_ORDER, PS_POLYNOMIAL_CHEB);
     x = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
     y = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
@@ -172,5 +172,5 @@
     psS32 memLeaks = 0;
 
-    myPoly = psPolynomial1DAlloc(POLY_ORDER+1, PS_POLYNOMIAL_CHEB);
+    myPoly = psPolynomial1DAlloc(POLY_ORDER, PS_POLYNOMIAL_CHEB);
     x = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
     y = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
Index: trunk/psLib/test/math/tst_psMinimizeVector2D_F32.c
===================================================================
--- trunk/psLib/test/math/tst_psMinimizeVector2D_F32.c	(revision 5070)
+++ trunk/psLib/test/math/tst_psMinimizeVector2D_F32.c	(revision 5090)
@@ -45,5 +45,5 @@
     double actualData;
 
-    psPolynomial2D *myPoly = psPolynomial2DAlloc(POLY_ORDER+1, POLY_ORDER+1, PS_POLYNOMIAL_ORD);
+    psPolynomial2D *myPoly = psPolynomial2DAlloc(POLY_ORDER, POLY_ORDER, PS_POLYNOMIAL_ORD);
     psVector *x = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
     psVector *y = psVectorAlloc(NUM_DATA, PS_TYPE_F64);
Index: trunk/psLib/test/math/tst_psMinimizeVector2D_F64.c
===================================================================
--- trunk/psLib/test/math/tst_psMinimizeVector2D_F64.c	(revision 5070)
+++ trunk/psLib/test/math/tst_psMinimizeVector2D_F64.c	(revision 5090)
@@ -45,5 +45,5 @@
     double actualData;
 
-    psPolynomial2D *myPoly = psPolynomial2DAlloc(POLY_ORDER+1, POLY_ORDER+1, PS_POLYNOMIAL_ORD);
+    psPolynomial2D *myPoly = psPolynomial2DAlloc(POLY_ORDER, POLY_ORDER, PS_POLYNOMIAL_ORD);
     psVector *x = psVectorAlloc(NUM_DATA, PS_TYPE_F32);
     psVector *y = psVectorAlloc(NUM_DATA, PS_TYPE_F32);