ocmath_renka_cline_eval

1 Description
2 Syntax
3 Parameters
4 Return
5 Examples
6 Remark
7 See Also
8 header to Included
9 Reference

Description

This subroutine computes the value at given points using the method of renka and cline.

Syntax

int ocmath_renka_cline_eval( ocmath_RenkaCline_Struct * comm, int n, double * px, double * py, double * pf )

Parameters

comm: [input] the structure computed from ocmath_renka_cline_interpolation()
n: [input] number of nodes and data values defining q.
px: [input] x cartesian coordinates of the points at which q are to be evaluated.
py: [input] y cartesian coordinates of the points at which q are to be evaluated.
pf: [output] value of q at (PX,PY)

Return

Err_Intp_No_Error: No error occur;

Err_Intp_Mem_Fail: Fails to allocate memory;

OR error from ocmath_renka_cline_interpolation().

Examples

EX1

#include <wks2mat.h>
void ocmath_renka_cline_eval_ex1()
{
        vector vX = { 1.7, 2.1, 3.9, 7.2, 8.6, 8.5, 7.3, 5.1, 2.8, 1.8, 1.7 };
    vector vY = { 3.2, 3.9, 4.9, 5.3, 5.5, 6.2, 6.5, 6.9, 7.5, 8.3, 9.4 }; 
    vector vZ = { 3.2, 3.9, 4.9, 5.3, 5.5, 6.2, 6.5, 6.9, 7.5, 8.3, 9.4 }; 
 
    double dXMin, dXMax, dYMin, dYMax;
    vX.GetMinMax(dXMin, dXMax);
    vY.GetMinMax(dYMin, dYMax);
 
    int rows = 20, cols = 20;
        MatrixPage mp;
        mp.Create("Origin");
        MatrixLayer ml = mp.Layers();
        MatrixObject mo = ml.MatrixObjects(0);
 
    Matrix& mat = mo.GetDataObject();
        mat.SetSize(rows, cols);
    mo.SetXY(dXMin, dYMin, dXMax, dYMax);
 
    vector vXGrid(rows*cols), vYGrid(rows*cols);
 
    int iRet = ocmath_mat_to_regular_xyz(NULL, rows, cols, dXMin, dXMax, dYMin, dYMax, vXGrid, vYGrid, NULL);
 
    ocmath_RenkaCline_Struct comm;
    ocmath_renka_cline_interpolation(vX.GetSize(), vX, vY, vZ, &comm);//Get comm from ocmath_renka_cline_interpolation
    iRet = ocmath_renka_cline_eval(&comm, rows*cols, vXGrid, vYGrid, mat);    
    out_int("iRet=", iRet);
    ocmath_renka_cline_struct_free(&comm);    
}