2.1.18.12 ocmath_modified_shepard_eval

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

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;
}