removed tabs and formatted c and h files with emacs.

Author: pseudotensor

advance.c

@@ -1,42 +1,42 @@
 #include "decs.h"
 
 void bcucof(FTYPE y[], FTYPE y1[], FTYPE y2[], FTYPE y12[], FTYPE d1, FTYPE d2,
-	FTYPE **c)
+            FTYPE **c)
 {
-	static int wt[16][16]=
-		{ 1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
-		0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,
-		-3,0,0,3,0,0,0,0,-2,0,0,-1,0,0,0,0,
-		2,0,0,-2,0,0,0,0,1,0,0,1,0,0,0,0,
-		0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,
-		0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,
-		0,0,0,0,-3,0,0,3,0,0,0,0,-2,0,0,-1,
-		0,0,0,0,2,0,0,-2,0,0,0,0,1,0,0,1,
-		-3,3,0,0,-2,-1,0,0,0,0,0,0,0,0,0,0,
-		0,0,0,0,0,0,0,0,-3,3,0,0,-2,-1,0,0,
-		9,-9,9,-9,6,3,-3,-6,6,-6,-3,3,4,2,1,2,
-		-6,6,-6,6,-4,-2,2,4,-3,3,3,-3,-2,-1,-1,-2,
-		2,-2,0,0,1,1,0,0,0,0,0,0,0,0,0,0,
-		0,0,0,0,0,0,0,0,2,-2,0,0,1,1,0,0,
-		-6,6,-6,6,-3,-3,3,3,-4,4,2,-2,-2,-2,-1,-1,
-		4,-4,4,-4,2,2,-2,-2,2,-2,-2,2,1,1,1,1};
-	int l,k,j,i;
-	FTYPE xx,d1d2,cl[16],x[16];
+  static int wt[16][16]=
+    { 1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
+      0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,
+      -3,0,0,3,0,0,0,0,-2,0,0,-1,0,0,0,0,
+      2,0,0,-2,0,0,0,0,1,0,0,1,0,0,0,0,
+      0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,
+      0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,
+      0,0,0,0,-3,0,0,3,0,0,0,0,-2,0,0,-1,
+      0,0,0,0,2,0,0,-2,0,0,0,0,1,0,0,1,
+      -3,3,0,0,-2,-1,0,0,0,0,0,0,0,0,0,0,
+      0,0,0,0,0,0,0,0,-3,3,0,0,-2,-1,0,0,
+      9,-9,9,-9,6,3,-3,-6,6,-6,-3,3,4,2,1,2,
+      -6,6,-6,6,-4,-2,2,4,-3,3,3,-3,-2,-1,-1,-2,
+      2,-2,0,0,1,1,0,0,0,0,0,0,0,0,0,0,
+      0,0,0,0,0,0,0,0,2,-2,0,0,1,1,0,0,
+      -6,6,-6,6,-3,-3,3,3,-4,4,2,-2,-2,-2,-1,-1,
+      4,-4,4,-4,2,2,-2,-2,2,-2,-2,2,1,1,1,1};
+  int l,k,j,i;
+  FTYPE xx,d1d2,cl[16],x[16];
 
-	d1d2=d1*d2;
-	for (i=1;i<=4;i++) {
-		x[i-1]=y[i];
-		x[i+3]=y1[i]*d1;
-		x[i+7]=y2[i]*d2;
-		x[i+11]=y12[i]*d1d2;
-	}
-	for (i=0;i<=15;i++) {
-		xx=0.0;
-		for (k=0;k<=15;k++) xx += wt[i][k]*x[k];
-		cl[i]=xx;
-	}
-	l=0;
-	for (i=1;i<=4;i++)
-		for (j=1;j<=4;j++) c[i][j]=cl[l++];
+  d1d2=d1*d2;
+  for (i=1;i<=4;i++) {
+    x[i-1]=y[i];
+    x[i+3]=y1[i]*d1;
+    x[i+7]=y2[i]*d2;
+    x[i+11]=y12[i]*d1d2;
+  }
+  for (i=0;i<=15;i++) {
+    xx=0.0;
+    for (k=0;k<=15;k++) xx += wt[i][k]*x[k];
+    cl[i]=xx;
+  }
+  l=0;
+  for (i=1;i<=4;i++)
+    for (j=1;j<=4;j++) c[i][j]=cl[l++];
 }
 /* (C) Copr. 1986-92 Numerical Recipes Software *1.@Q.. */

bcucof.c

@@ -2,30 +2,30 @@
 #include "decs.h"
 
 void bcuint(FTYPE y[], FTYPE y1[], FTYPE y2[], FTYPE y12[], FTYPE x1l,
-	FTYPE x1u, FTYPE x2l, FTYPE x2u, FTYPE x1, FTYPE x2, FTYPE *ansy,
-	FTYPE *ansy1, FTYPE *ansy2)
+            FTYPE x1u, FTYPE x2l, FTYPE x2u, FTYPE x1, FTYPE x2, FTYPE *ansy,
+            FTYPE *ansy1, FTYPE *ansy2)
 {
-	void bcucof(FTYPE y[], FTYPE y1[], FTYPE y2[], FTYPE y12[], FTYPE d1,
-		FTYPE d2, FTYPE **c);
-	int i;
-	FTYPE t,u,d1,d2,**c;
+  void bcucof(FTYPE y[], FTYPE y1[], FTYPE y2[], FTYPE y12[], FTYPE d1,
+              FTYPE d2, FTYPE **c);
+  int i;
+  FTYPE t,u,d1,d2,**c;
 
-	c=fmatrix(1,4,1,4);
-	d1=x1u-x1l;
-	d2=x2u-x2l;
-	bcucof(y,y1,y2,y12,d1,d2,c);
-	if (x1u == x1l || x2u == x2l) nrerror("Bad input in routine bcuint");
-	t=(x1-x1l)/d1;
-	u=(x2-x2l)/d2;
-	*ansy=(*ansy2)=(*ansy1)=0.0;
-	for (i=4;i>=1;i--) {
-		*ansy=t*(*ansy)+((c[i][4]*u+c[i][3])*u+c[i][2])*u+c[i][1];
-		*ansy2=t*(*ansy2)+(3.0*c[i][4]*u+2.0*c[i][3])*u+c[i][2];
-		*ansy1=u*(*ansy1)+(3.0*c[4][i]*t+2.0*c[3][i])*t+c[2][i];
-	}
-	*ansy1 /= d1;
-	*ansy2 /= d2;
-	free_fmatrix(c,1,4,1,4);
+  c=fmatrix(1,4,1,4);
+  d1=x1u-x1l;
+  d2=x2u-x2l;
+  bcucof(y,y1,y2,y12,d1,d2,c);
+  if (x1u == x1l || x2u == x2l) nrerror("Bad input in routine bcuint");
+  t=(x1-x1l)/d1;
+  u=(x2-x2l)/d2;
+  *ansy=(*ansy2)=(*ansy1)=0.0;
+  for (i=4;i>=1;i--) {
+    *ansy=t*(*ansy)+((c[i][4]*u+c[i][3])*u+c[i][2])*u+c[i][1];
+    *ansy2=t*(*ansy2)+(3.0*c[i][4]*u+2.0*c[i][3])*u+c[i][2];
+    *ansy1=u*(*ansy1)+(3.0*c[4][i]*t+2.0*c[3][i])*t+c[2][i];
+  }
+  *ansy1 /= d1;
+  *ansy2 /= d2;
+  free_fmatrix(c,1,4,1,4);
 }
 #undef NRANSI
 /* (C) Copr. 1986-92 Numerical Recipes Software *1.@Q.. */