Have to be more careful with the contiguous keyword given the coeffs …

…are now pointers
stevendargaville · Nov 8, 2024 · 15c04c8 · 15c04c8
1 parent 736c999
commit 15c04c8
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/Approx_Inverse_Setup.F90 b/src/Approx_Inverse_Setup.F90
@@ -264,7 +264,7 @@ subroutine finish_approximate_inverse(matrix, inverse_type, &
       integer, intent(in)                               :: inverse_type, poly_order
       integer, intent(in)                               :: inverse_sparsity_order
       type(tsqr_buffers), intent(inout)                 :: buffers      
-      real, dimension(:, :), pointer, intent(inout)     :: coefficients
+      real, dimension(:, :), pointer, contiguous, intent(inout) :: coefficients
       logical, intent(in)                               :: matrix_free
       type(tMat), intent(inout)                         :: reuse_mat, inv_matrix
 

diff --git a/src/Gmres_Poly_Data_Type.F90 b/src/Gmres_Poly_Data_Type.F90
@@ -19,7 +19,7 @@ module gmres_poly_data_type
       ! Coefficients for the gmres polynomial
       ! If using the newton basis this has two columns with the real 
       ! and imaginary roots
-      real, pointer, dimension(:, :) :: coefficients => null()
+      real, pointer, dimension(:, :), contiguous :: coefficients => null()
 
       ! Asynchronous buffers for the TSQR
       type(tsqr_buffers) :: buffers