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_psFunc09.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
    … …  
    8888 
    8989    // Allocate polynomial structure
    90     psPolynomial2D*  polyOrd = psPolynomial2DAlloc(TERMS, TERMS, PS_POLYNOMIAL_ORD);
    91     psPolynomial2D*  polyCheb = psPolynomial2DAlloc(TERMS, TERMS, PS_POLYNOMIAL_CHEB);
     90    psPolynomial2D*  polyOrd = psPolynomial2DAlloc(TERMS-1, TERMS-1, PS_POLYNOMIAL_ORD);
     91    psPolynomial2D*  polyCheb = psPolynomial2DAlloc(TERMS-1, TERMS-1, PS_POLYNOMIAL_CHEB);
    9292    // Set polynomial members
    9393    for(psS32 i = 0; i < TERMS; i++) {
    … …  
    118118 
    119119    // Allocate polynomial with invalid type
    120     polyOrd = psPolynomial2DAlloc(TERMS, TERMS, 99);
     120    polyOrd = psPolynomial2DAlloc(TERMS-1, TERMS-1, 99);
    121121    // Attempt to evaluation invalid polynomial type
    122122    psLogMsg(__func__,PS_LOG_INFO,"Following should generate error message invalid type");
    … …  
    138138
    139139    // Allocate polynomial structure
    140     psPolynomial2D*  polyOrd = psPolynomial2DAlloc(TERMS, TERMS, PS_POLYNOMIAL_ORD);
    141     psPolynomial2D*  polyCheb = psPolynomial2DAlloc(TERMS, TERMS, PS_POLYNOMIAL_CHEB);
     140    psPolynomial2D*  polyOrd = psPolynomial2DAlloc(TERMS-1, TERMS-1, PS_POLYNOMIAL_ORD);
     141    psPolynomial2D*  polyCheb = psPolynomial2DAlloc(TERMS-1, TERMS-1, PS_POLYNOMIAL_CHEB);
    142142    // Set polynomial members
    143143    for(psS32 i = 0; i < TERMS; i++) {
    … …  
    168168
    169169    // Allocate polynomial with invalid type
    170     polyOrd = psPolynomial2DAlloc(TERMS, TERMS, 99);
     170    polyOrd = psPolynomial2DAlloc(TERMS-1, TERMS-1, 99);
    171171    // Attempt to evaluation invalid polynomial type
    172172    psLogMsg(__func__,PS_LOG_INFO,"Following should generate error message invalid type");
    … …  
    184184{
    185185    // Allocate polynomial
    186     psPolynomial2D* polyOrd = psPolynomial2DAlloc(TERMS,TERMS,PS_POLYNOMIAL_ORD);
    187     psPolynomial2D* polyCheb = psPolynomial2DAlloc(TERMS,TERMS,PS_POLYNOMIAL_CHEB);
     186    psPolynomial2D* polyOrd = psPolynomial2DAlloc(TERMS-1,TERMS-1,PS_POLYNOMIAL_ORD);
     187    psPolynomial2D* polyCheb = psPolynomial2DAlloc(TERMS-1,TERMS-1,PS_POLYNOMIAL_CHEB);
    188188 
    189189    // Set polynomial members
    … …  
    297297{
    298298    // Allocate polynomial
    299     psPolynomial2D* polyOrd = psPolynomial2DAlloc(TERMS,TERMS,PS_POLYNOMIAL_ORD);
    300     psPolynomial2D* polyCheb = psPolynomial2DAlloc(TERMS,TERMS,PS_POLYNOMIAL_CHEB);
     299    psPolynomial2D* polyOrd = psPolynomial2DAlloc(TERMS-1,TERMS-1,PS_POLYNOMIAL_ORD);
     300    psPolynomial2D* polyCheb = psPolynomial2DAlloc(TERMS-1,TERMS-1,PS_POLYNOMIAL_CHEB);
    301301
    302302    // Set polynomial members
Note: See TracChangeset for help on using the changeset viewer.