2.1.18.12 ocmath_modified_shepard_eval

Description

This subroutine computes the value at given points using the modified shepard's method.

Syntax

int ocmath_modified_shepard_eval( ocmath_ModShepard_Struct * comm, int n, double * px, double * py, double * pf )

Parameters

comm

[input] the structure computed from ocmath_modified_shepard_interpolant()

n

[input] number of nodes and 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>
int ocmath_Modified_Shepard_ex1()
{
	double x1[10], y1[10], z1[10] ;
    int n = 10 ;
	for(int ii = 0; ii < n; ii++)
	{
		x1[ii] = ii;
		y1[ii] = 1;
		z1[ii] = ii+1;
	}
		
	ocmath_ModShepard_Opt optional;
	optional.nq = 5;
	optional.nw = 7;
	ocmath_ModShepard_Struct comm;
	int nRet = ocmath_modified_shepard_interpolant(n, x1, y1, z1,&comm, NULL);
	if(nRet < 0)
		return 0;

	double px[5] = { .1,.3,.5,.7,.9 };
	double py[5] = { .2,.4,.6,.8, .4 };
    int	m = 5;

	double pf[5];
    nRet = ocmath_modified_shepard_eval(&comm, m, px, py, pf);
	if(nRet <0)
		return 0;
	
	return 1;
}

Remark

See Also

Header to Include

wks2mat.h

Reference