My Documents/Codes/MD++.svn/trunk/src/stringmethod.cpp File Reference

#include "stringmethod.h"

Go to the source code of this file.

Defines
#define	get_Xi0(y0, y1, y2, y3)
#define	get_XiN(y_nm3, y_nm2, y_nm1, y_n)
#define	get_cubicspline(h_i, w, wbar, sigma_i, sigma_im1, y_i, y_im1)

---

Define Documentation

#define get_cubicspline	(	h_i,
		w,
		wbar,
		sigma_i,
		sigma_im1,
		y_i,
		y_im1		)

Value:

protect(\
        y_csp = w*y_i + wbar*y_im1 + SQR(h_i)*((pow(w,3)-w)*sigma_i + (pow(wbar,3)-wbar)*sigma_im1); \
        )

#define get_Xi0	(	y0,
		y1,
		y2,
		y3		)

Value:

protect(\
        Xi0 = DL(0)*( y0/(_nebinterp[0]-_nebinterp[1])/(_nebinterp[0]-_nebinterp[2])/(_nebinterp[0]-_nebinterp[3]) \
                    + y1/(_nebinterp[1]-_nebinterp[0])/(_nebinterp[1]-_nebinterp[2])/(_nebinterp[1]-_nebinterp[3]) \
                    + y2/(_nebinterp[2]-_nebinterp[0])/(_nebinterp[2]-_nebinterp[1])/(_nebinterp[2]-_nebinterp[3]) \
                    + y3/(_nebinterp[3]-_nebinterp[0])/(_nebinterp[3]-_nebinterp[1])/(_nebinterp[3]-_nebinterp[2])); \
        )

#define get_XiN	(	y_nm3,
		y_nm2,
		y_nm1,
		y_n		)

Value:

protect(\
        XiN = DU(_CHAINLENGTH-1)*( y_nm3/(_nebinterp[_CHAINLENGTH-3]-_nebinterp[_CHAINLENGTH-2])/(_nebinterp[_CHAINLENGTH-3]-_nebinterp[_CHAINLENGTH-1])/(_nebinterp[_CHAINLENGTH-3]-_nebinterp[_CHAINLENGTH]) \
                    + y_nm2/(_nebinterp[_CHAINLENGTH-2]-_nebinterp[_CHAINLENGTH-3])/(_nebinterp[_CHAINLENGTH-2]-_nebinterp[_CHAINLENGTH-1])/(_nebinterp[_CHAINLENGTH-2]-_nebinterp[_CHAINLENGTH]) \
                    + y_nm1/(_nebinterp[_CHAINLENGTH-1]-_nebinterp[_CHAINLENGTH-3])/(_nebinterp[_CHAINLENGTH-1]-_nebinterp[_CHAINLENGTH-2])/(_nebinterp[_CHAINLENGTH-1]-_nebinterp[_CHAINLENGTH]) \
                    + y_n/(_nebinterp[_CHAINLENGTH]-_nebinterp[_CHAINLENGTH-3])/(_nebinterp[_CHAINLENGTH]-_nebinterp[_CHAINLENGTH-2])/(_nebinterp[_CHAINLENGTH]-_nebinterp[_CHAINLENGTH-1])); \
        )