20#include "hdrl_cat_polynm.h"
22#include "hdrl_cat_solve.h"
55 double xdat[],
double xcor[], cpl_size n,
double polycf[], cpl_size m, cpl_size ilim)
59 double a[sizeAray][sizeAray];
61 for (cpl_size i = 0; i < sizeAray; i++) {
65 for (cpl_size j = 0; j < sizeAray; j++) {
71 for (cpl_size i = 0; i < n; i++) {
73 for (cpl_size k = 0; k < m; k++) {
78 temp = pow(xcor[i], (
double)(k + ilim));
81 b[k] += xdat[i] * temp;
83 for (cpl_size j = 0; j <= k; j++) {
87 if (k+j+2*ilim != 0) {
88 temp = pow(xcor[i], (
double)(k + j + 2 * ilim));
96 for (cpl_size k = 1; k < m; k++) {
97 for (cpl_size j = 0; j < k; j++) {
105 for (cpl_size i = 0; i < m; i++) {
109 return CPL_ERROR_NONE;
cpl_error_code hdrl_polynm(double xdat[], double xcor[], cpl_size n, double polycf[], cpl_size m, cpl_size ilim)
Work out the median seeing.
cpl_error_code hdrl_solve(double a[25][25], double b[25], cpl_size m)
Use Gauss-Jordan elimination to solve ax=b.