1 SUBROUTINE QIMUDS
(Q
,F
,MAXA
,NN
,LENAA
)
3 C computes f := [Q**(-1)] * f .
7 C Q is the preconditioning matrix, and contains an approximate
10 C f is the right hand side vector, Q z = f .
12 C MAXA, NN, LENAA describe Q in packed skyline storage format.
16 C Q, MAXA, NN, LENAA are unchanged.
18 C f contains the solution z of Q z = f .
23 INTEGER LENAA
,LENQ
,NN
,MAXA
(NN
+2),NQ
24 DOUBLE PRECISION Q
(LENAA
+NN
+1),F
(NN
+1)
29 CALL SOLVDS
(NQ
,Q
,LENQ
,MAXA
,F
)