Add missing SUNMATRIX module initialization in the DENSE solver case (solverType==3).

Sam Cox · Sam Cox · commit 8869c1f26bdb · 2018-08-22T10:39:27.000+01:00
diff --git a/src/atchem2.f90 b/src/atchem2.f90
@@ -346,6 +346,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
