qfunc: disentangle the API from the global solver API

1 files changed, 3 insertions, 1 deletions

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++)