RH/bezier.h at master · DKISTDC/RH · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
/* ------- file: -------------------------- bezier.h ----------------

   Cubic DELO-Bezier (polarized) and cubic short-char Bezier solvers.

   References: de la Cruz Rodriguez & Piskunov (2013), Auer (2003)
               (Derivatives) Fritsch & Butland (1984),

   Coded by J. de la Cruz Rodriguez (ISP-SU 2017)

       --------------------------                      ----------RH-- */


#ifndef BEZIER_H
#define BEZIER_H

/* --- Macros --                                       -------------- */

#define SWAP(a, b, c) (c=a,a=b,b=c)


/* --- Prototypes auxiliary functions -                -------------- */

static inline double signFortran(const double val)
{
  return ((val >= 0.0)? 1.0 : -1.0);
}

/* --- Matrix inversion Cramer method SSE instructions -- ----------- */

void SIMD_MatInv(float* src);

/* --- Matrix inversion Shipley-Coleman (not used) --  -------------- */

void m4inv(double MI[4][4]);


void m4v(float a[4][4], double b[4], double c[4]);
void m4m(double a[4][4], double b[4][4], double c[4][4]);

void cent_deriv_vec(double wprime[4], double dsup, double dsdn,
		    double chiup[4], double chic[4], double chidn[4]);

void cent_deriv_mat(double wprime[4][4], double dsup, double dsdn,
		    double chiup[4][4], double chic[4][4],
		    double chidn[4][4]);

double cent_deriv(double dsup,double dsdn, double chiup,
		  double chic, double chidn);
void Svec(int k, double **S, double *Sf);

void Bezier3_coeffs(double dt, double *alpha, double *beta,
		    double *gamma, double *theta, double *eps);

void solveLinearFast(double A[4][4], double B[4]);

#endif /* !__BEZIER_H__ */

/* ---------------------------------------- bezier.h ---------------- */