diff options
Diffstat (limited to 'solve.c')
-rw-r--r-- | solve.c | 4 |
1 files changed, 3 insertions, 1 deletions
@@ -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++) |