kallemov / fast-spatial-interpolation-for-irregular-grid Goto Github PK

#define qqa (360.)
#define zzi (7.)
double sbf[4] = {100.,102,103.,104.};
double sbxa[4] = {qqa-120.-zzi,qqa-120.-zzi,360. -119.5,360. -119.5};
double sbya[4] = {38.+zzi,38.25+zzi,38.25,38. };
double sbx = qqa-119.85;
double sby = 38.18;
double sbv = 0;
// sb_init(sbf, sbxa, sbya);
// sbv = sb_interp(sbx, sby);
// when qqa = 0
// zzi of 7 returns +99
// zzi of 5 returns +135
// zzi of 2 returns -442
// zzi of 1 returns -4208
// zzi of 0 returns +810576
// when qqa = 360
// zzi of 7 returns -786
// zzi of 5 returns +2900
// zzi of 2 returns +3605
// zzi of 1 returns +225578
// zzi of 0 returns -93498082
``

Can the interpolation routines be formulated slightly differently so as to provide more numerical stability?