From 63b25eca76d8481a0fdad339464d1eb617826583 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Sun, 28 Aug 2016 10:07:20 +0200 Subject: qfunc: disentangle the API from the global solver API --- solve.c | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) (limited to 'solve.c') diff --git a/solve.c b/solve.c index 64dcb55..5803d87 100644 --- a/solve.c +++ b/solve.c @@ -32,6 +32,7 @@ #include "brill_data.h" #include "internal.h" +#include "qfunc.h" static int brill_construct_matrix(const BDContext *bd, double *mat, double *rhs) @@ -76,7 +77,8 @@ static int brill_construct_matrix(const BDContext *bd, double *mat, for (idx_grid_rho = 0; idx_grid_rho < s->nb_colloc_points[0]; idx_grid_rho++) { double x_val = s->basis[0]->colloc_point(s->nb_coeffs[0], idx_grid_rho, sf[0]); - double d2q = s->q_func->d2q_rho(bd, x_val, z_val) + s->q_func->d2q_z(bd, x_val, z_val); + double d2q = bdi_qfunc_eval(s->qfunc, x_val, z_val, 2) + + bdi_qfunc_eval(s->qfunc, x_val, z_val, 6); int idx_grid = idx_grid_z * s->nb_colloc_points[0] + idx_grid_rho; for (idx_coeff_z = 0; idx_coeff_z < s->nb_coeffs[1]; idx_coeff_z++) -- cgit v1.2.3