--- a/LinPack.mod.m2pp
+++ b/LinPack.mod.m2pp
@@ -7,9 +7,9 @@
(* *)
(* m2pp -D __{Parameter}__ < LinPack.mod.m2pp > LibPack.mod *)
(* *)
- (* mit Parameter = {inline|BLAS|BLASF77} erzeugt werden koennen. *)
+ (* mit Parameter = {INLINE|BLAS|BLASF77} erzeugt werden koennen. *)
(* *)
- (* inline : Alle Schleifen werden expliziet ausgefuehrt *)
+ (* INLINE : Alle Schleifen werden expliziet ausgefuehrt *)
(* BLAS : Schleifen werden durch Routinen aus LibDBlas ausgefuehrt *)
(* BLASF77: Schleifen werden durch Routinen aus der Fortran 77 *)
(* Bibiliothek BLAS ausgefuehrt - diese muessen mit gelinkt *)
@@ -51,7 +51,7 @@
(* $Id: LinPack.mod.m2pp,v 1.3 2018/01/16 09:19:51 mriedl Exp mriedl $ *)
(*------------------------------------------------------------------------*)
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
(* HINT: This version uses inlined code. *)
(* This is most probably slower than using BLAS level 1 routines. *)
<* END *>
@@ -82,7 +82,7 @@
VAR t : LONGREAL;
j,k,kp1,l,nm1 : INTEGER;
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
dmax : LONGREAL;
ii : INTEGER;
<* END *>
@@ -107,7 +107,7 @@
(* find l := pivot index *)
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
dmax:=ABS(A[k,k]);
l:=k;
FOR ii:=k TO n-1 DO
@@ -142,7 +142,7 @@
(* compute multipliers *)
t := -1.0 / A[k,k];
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
FOR ii:=k+1 TO n-1 DO A[k,ii]:=A[k,ii]*t; END;
<* END *>
<* IF (__BLAS__) THEN *>
@@ -160,7 +160,7 @@
A[j,l] := A[j,k];
A[j,k] := t;
END;
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
FOR ii:=k+1 TO n-1 DO A[j,ii]:=A[j,ii] + A[k,ii]*t ; END;
<* END *>
<* IF (__BLAS__) THEN *>
@@ -196,7 +196,7 @@
VAR t : LONGREAL;
k,l,nm1 : INTEGER;
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
dot : LONGREAL;
ii : INTEGER;
<* END *>
@@ -219,7 +219,7 @@
(* first solve trans(u)*y = b *)
FOR k:=0 TO n-1 DO
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
(* t := LibDBlasM2.ddot(k,A[k],1,B,1); *)
t:=0.0; FOR ii:=0 TO k-1 DO t:=t + A[k,ii]*B[ii]; END;
<* END *>
@@ -241,7 +241,7 @@
k := n - kb - 1;
*)
FOR k:=nm1-1 TO 0 BY -1 DO (* TO 0 statt TO 1, MRi 26.12.16 *)
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
dot:=0.0; FOR ii:=k+1 TO n-1 DO dot:=dot + A[k,ii]*B[ii]; END;
B[k]:=B[k] + dot;
<* END *>
@@ -271,7 +271,7 @@
B[l] := B[k];
B[k] := t;
END;
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
FOR ii:=k+1 TO n-1 DO B[ii]:=B[ii] + t*A[k,ii]; END;
<* END *>
<* IF (__BLAS__) THEN *>
@@ -292,7 +292,7 @@
FOR k:=n-1 TO 0 BY -1 DO
B[k] := B[k] / A[k,k];
t := -B[k];
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
(* LibDBlasM2.daxpy(k,t,A[k],1,B,1); *)
FOR ii:=0 TO k-1 DO B[ii]:=B[ii] + t*A[k,ii]; END;
<* END *>
@@ -321,7 +321,7 @@
<* IF (__BLASF77__) THEN *>
ione,num : INTEGER;
<* END *>
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
ii : INTEGER;
<* END *>
BEGIN
@@ -366,7 +366,7 @@
FOR k:=0 TO n-1 DO
A[k,k] := 1.0 / A[k,k];
t := -A[k,k];
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
FOR ii:=0 TO k-1 DO A[k,ii]:=A[k,ii]*t; END;
(* LibDBlasM2.dscal(k,t,A[k],1); *)
(* LibDBlas.dscal(k,t,A[k,0],1); *)
@@ -381,7 +381,7 @@
FOR j:=k+1 TO n-1 DO
t := A[j,k];
A[j,k] := 0.0;
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
FOR ii:=0 TO k DO A[j,ii]:=A[j,ii] + t*A[k,ii]; END;
(* LibDBlasM2.daxpy(k+1,t,A[k],1,A[j],1); *)
(* LibDBlas.daxpy(k+1,t,A[k,0],1,A[j,0],1); *)
@@ -405,7 +405,7 @@
END;
FOR j:=k+1 TO n-1 DO
t := work[j];
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
FOR ii:=0 TO n-1 DO A[k,ii]:=A[k,ii] + t*A[j,ii]; END;
(* LibDBlasM2.daxpy(n,t,A[j],1,A[k],1); *)
(* LibDBlas.daxpy(n,t,A[j,0],1,A[k,0],1); *)
@@ -420,7 +420,7 @@
END;
l := ipvt[k];
IF (l # k) THEN
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
FOR ii:=0 TO n-1 DO
t := A[k,ii];
A[k,ii] := A[l,ii];
@@ -449,7 +449,7 @@
VAR j,jj,k,kj,kk : INTEGER;
s,t : LONGREAL;
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
dot : LONGREAL;
ij,jjk : INTEGER;
<* END *>
@@ -469,7 +469,7 @@
kj := jj;
kk := 0;
FOR k:=1 TO j-1 DO
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
dot:=0.0; jjk:=jj;
FOR ij:=kk TO kk+k-2 DO dot:=dot + AP[ij]*AP[jjk]; INC(jjk); END;
t := AP[kj] - dot;
@@ -506,7 +506,7 @@
VAR i,ii,j,j1,jj,k,k1,kj,kk : INTEGER;
t : LONGREAL;
exit : BOOLEAN;
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
ij,j1j,k1j : INTEGER;
<* END *>
<* IF (__BLASF77__) THEN *>
@@ -553,7 +553,7 @@
INC(kk,k);
AP[kk-1] := 1.0 / AP[kk-1];
t := -AP[kk-1];
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
FOR ij:=k1 TO k1+k-2 DO AP[ij]:=AP[ij]*t; END;
(* dscal(k-1,t,AP[k1],1); *)
<* END *>
@@ -569,7 +569,7 @@
FOR j:=k+1 TO n DO
t := AP[kj];
AP[kj] := 0.0;
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
k1j:=k1;
FOR ij:=j1 TO j1+k-1 DO AP[ij]:=AP[ij] + t*AP[k1j]; INC(k1j); END;
(* daxpy(k,t,AP[k1],1,AP[j1],1); *)
@@ -595,7 +595,7 @@
k1 := 0;
FOR k:=1 TO j-1 DO
t := AP[kj];
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
j1j:=j1;
FOR ij:=k1 TO k1+k-1 DO AP[ij]:=AP[ij] + t*AP[j1j]; INC(j1j); END;
(* daxpy(k,t,AP[j1],1,AP[k1],1); *)
@@ -612,7 +612,7 @@
END;
INC(jj,j);
t := AP[jj-1];
-<* IF (__inline__) THEN *>
+<* IF (__INLINE__) THEN *>
FOR ij:=j1 TO jj-1 DO AP[ij]:=AP[ij]*t; END;
(* dscal(j,t,AP[j1],1); *)
<* END *>