IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Ignore:
Timestamp:
Sep 21, 2005, 4:47:16 PM (21 years ago)
Author:
gusciora
Message:

The nX terms for the polynomials now refer to the order, not the number
of terms.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/psLib/test/math/tst_psFunc10.c

    r4731 r5090  
    44*  ORD and CHEB type polynomials.
    55*
    6 *  @version  $Revision: 1.3 $  $Name: not supported by cvs2svn $
    7 *  @date  $Date: 2005-08-08 21:44:09 $
     6*  @version  $Revision: 1.4 $  $Name: not supported by cvs2svn $
     7*  @date  $Date: 2005-09-22 02:47:16 $
    88*
    99* Copyright 2004-2005 Maui High Performance Computing Center, Univ. of Hawaii
    … …  
    139139 
    140140    // Allocate polynomial structure
    141     psPolynomial3D*  polyOrd = psPolynomial3DAlloc(TERMS, TERMS, TERMS, PS_POLYNOMIAL_ORD);
    142     psPolynomial3D*  polyCheb = psPolynomial3DAlloc(TERMS, TERMS, TERMS, PS_POLYNOMIAL_CHEB);
     141    psPolynomial3D*  polyOrd = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_ORD);
     142    psPolynomial3D*  polyCheb = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1, PS_POLYNOMIAL_CHEB);
    143143    // Set polynomial members
    144144    for(psS32 i = 0; i < TERMS; i++) {
    … …  
    173173 
    174174    // Allocate polynomial with invalid type
    175     polyOrd = psPolynomial3DAlloc(TERMS, TERMS, TERMS, 99);
     175    polyOrd = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1, 99);
    176176    // Attempt to evaluation invalid polynomial type
    177177    psLogMsg(__func__,PS_LOG_INFO,"Following should generate error message invalid type");
    … …  
    193193
    194194    // Allocate polynomial structure
    195     psPolynomial3D*  polyOrd = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_ORD);
    196     psPolynomial3D*  polyCheb = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_CHEB);
     195    psPolynomial3D*  polyOrd = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_ORD);
     196    psPolynomial3D*  polyCheb = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_CHEB);
    197197    // Set polynomial members
    198198    for(psS32 i = 0; i < TERMS; i++) {
    … …  
    227227
    228228    // Allocate polynomial with invalid type
    229     polyOrd = psPolynomial3DAlloc(TERMS, TERMS, TERMS, 99);
     229    polyOrd = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1, 99);
    230230    // Attempt to evaluation invalid polynomial type
    231231    psLogMsg(__func__,PS_LOG_INFO,"Following should generate error message invalid type");
    … …  
    243243{
    244244    // Allocate polynomial
    245     psPolynomial3D* polyOrd  = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_ORD);
    246     psPolynomial3D* polyCheb = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_CHEB);
     245    psPolynomial3D* polyOrd  = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_ORD);
     246    psPolynomial3D* polyCheb = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_CHEB);
    247247 
    248248    // Set polynomial members
    … …  
    378378{
    379379    // Allocate polynomial
    380     psPolynomial3D* polyOrd = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_ORD);
    381     psPolynomial3D* polyCheb = psPolynomial3DAlloc(TERMS,TERMS,TERMS,PS_POLYNOMIAL_CHEB);
     380    psPolynomial3D* polyOrd = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_ORD);
     381    psPolynomial3D* polyCheb = psPolynomial3DAlloc(TERMS-1, TERMS-1, TERMS-1,PS_POLYNOMIAL_CHEB);
    382382
    383383    // Set polynomial members
Note: See TracChangeset for help on using the changeset viewer.