Fix examples BAND _nbdr2 returning error info = -9

Both examples access the M matrix, so we initialize M to be the identity matrix. See #443 for discussion.
2023-12-20 09:41:27 +01:00
parent 31b50762a5
commit 0eaf38bd51
2 changed files with 2 additions and 0 deletions
@@ -206,6 +206,7 @@ c
      idiag = kl+ku+1
      do 30 j = 1, n
         a(idiag,j) = (4.0E+0, 0.0E+0)  / h2
+         m(idiag,j) = (1.0E+0, 0.0E+0)
  30  continue
 c
 c     %-------------------------------------%
@@ -206,6 +206,7 @@ c
      idiag = kl+ku+1
      do 30 j = 1, n
         a(idiag,j) = (4.0D+0, 0.0D+0)  / h2
+         m(idiag,j) = (1.0D+0, 0.0D+0)
  30  continue
 c
 c     %-------------------------------------%