diff options
| -rw-r--r-- | ORDERINGS/ccn.f | 8 |
1 files changed, 5 insertions, 3 deletions
diff --git a/ORDERINGS/ccn.f b/ORDERINGS/ccn.f index eee5836..037cfa9 100644 --- a/ORDERINGS/ccn.f +++ b/ORDERINGS/ccn.f @@ -90,7 +90,7 @@ C Laura C. Dutto - email: dutto@cerca.umontreal.ca c July 1992 - Update: March 1994 C----------------------------------------------------------------------- integer izs(nw), lpw(n), nsbloc(0:nblcmx), ia(n+1), ja(*) - real*8 amat(*) + real*8 amat(*), iziama(1) logical impr character*6 chsubr C----------------------------------------------------------------------- @@ -147,10 +147,12 @@ c..........We copy IA and JA on IAT and JAT respectively * ja, ia, izs(ilpw), izs(ilpw), job) ipos = 1 c..........We sort columns inside JA. - call csrcsc(n, job, ipos, amat, ja, ia, izs(iamat), + iziama(1) = izs(iamat) + call csrcsc(n, job, ipos, amat, ja, ia, iziama, * izs(ijat), izs(iiat)) - call csrcsc(n, job, ipos, izs(iamat), izs(ijat), izs(iiat), + call csrcsc(n, job, ipos, iziama, izs(ijat), izs(iiat), * amat, ja, ia) + izs(iamat) = iziama(1) endif c.....We modify the ordering of unknowns in LPW call compos(n, lpw, izs(ilpw)) |