diff options
Diffstat (limited to 'src/elliptic/dense_Jacobian.cc')
-rw-r--r-- | src/elliptic/dense_Jacobian.cc | 11 |
1 files changed, 7 insertions, 4 deletions
diff --git a/src/elliptic/dense_Jacobian.cc b/src/elliptic/dense_Jacobian.cc index 3e4e582..4fd9e66 100644 --- a/src/elliptic/dense_Jacobian.cc +++ b/src/elliptic/dense_Jacobian.cc @@ -165,8 +165,9 @@ void dense_Jacobian::zero_matrix() // // This function constructs a dense_Jacobian__LAPACK object. // -dense_Jacobian__LAPACK::dense_Jacobian__LAPACK(patch_system& ps, - bool print_msg_flag /* = false */) +dense_Jacobian__LAPACK::dense_Jacobian__LAPACK + (patch_system& ps, + bool print_msg_flag /* = false */) : dense_Jacobian(ps, print_msg_flag), pivot_(new integer[ N_rows_]), iwork_(new integer[ N_rows_]), @@ -201,8 +202,10 @@ delete[] pivot_; // It returns the estimated infinity-norm condition number of the linear // system, or 0.0 if the matrix is numerically singular // -fp dense_Jacobian__LAPACK::solve_linear_system(int rhs_gfn, int x_gfn, - bool print_msg_flag) +fp dense_Jacobian__LAPACK::solve_linear_system + (int rhs_gfn, int x_gfn, + const struct linear_solver_pars& pars, + bool print_msg_flag) { const fp *rhs = ps_.gridfn_data(rhs_gfn); fp *x = ps_.gridfn_data(x_gfn); |