Add missing SUNMATRIX module initialization in the DENSE solver case …

…(solverType==3).
AtChem · Aug 15, 2018 · def0415 · def0415
1 parent 7f0f9ec
commit def0415
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/atchem2.f90 b/src/atchem2.f90
@@ -351,6 +351,7 @@ end subroutine FCVFUN
   else if ( solverType == 3 ) then
     call FSUNDENSEMATINIT( 1, neq,  neq, ier )
     ! call FCVDENSE( neq, ier )
+    call FSUNDENSELINSOLINIT( 1, ier )
     call FCVDLSINIT( ier )
     ! UNEXPECTED SOLVER TYPE
   else