===================================================================
RCS file: /home/cvs/OpenXM_contrib2/asir2000/builtin/numerical.c,v
retrieving revision 1.6
retrieving revision 1.7
diff -u -p -r1.6 -r1.7
--- OpenXM_contrib2/asir2000/builtin/numerical.c	2003/03/07 06:39:55	1.6
+++ OpenXM_contrib2/asir2000/builtin/numerical.c	2018/03/29 01:32:50	1.7
@@ -45,7 +45,7 @@
  * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE,
  * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE.
  *
- * $OpenXM: OpenXM_contrib2/asir2000/builtin/numerical.c,v 1.5 2003/03/07 03:12:24 noro Exp $ 
+ * $OpenXM: OpenXM_contrib2/asir2000/builtin/numerical.c,v 1.6 2003/03/07 06:39:55 noro Exp $ 
 */
 #include "ca.h"
 #include "parse.h"
@@ -56,52 +56,52 @@
 void Pflinsolve();
 
 struct ftab numerical_tab[] = {
-	{"flinsolve",Pflinsolve,2},
-	{0,0,0},
+  {"flinsolve",Pflinsolve,2},
+  {0,0,0},
 };
 
 void Pflinsolve(arg,rp)
 NODE arg;
 VECT *rp;
 {
-	MAT mat;
-	VECT b,r;
-	int i,j;
-	integer n,nrhs,lda,ldb,info;
-	integer *ipiv;
-	doublereal *m_array,*b_array;
-	Num **mb;
-	Num *vb;
-	Real *rb;
-	char trans;
+  MAT mat;
+  VECT b,r;
+  int i,j;
+  integer n,nrhs,lda,ldb,info;
+  integer *ipiv;
+  doublereal *m_array,*b_array;
+  Num **mb;
+  Num *vb;
+  Real *rb;
+  char trans;
 
-	mat=(MAT)ARG0(arg); b=(VECT)ARG1(arg);
-	n=mat->row;
-	if ( n != mat->col )
-		error("flinsolve : non-square matrix");
-	if ( n != b->len )
-		error("flinsolve : inconsistent rhs");
-	m_array = (doublereal *)MALLOC_ATOMIC(n*n*sizeof(double));
-	b_array = (doublereal *)MALLOC_ATOMIC(n*sizeof(double));
-	ipiv = (integer *)MALLOC_ATOMIC(n*sizeof(integer));
-	for ( i = 0, mb = (Num **)mat->body; i < n; i++ )
-		for ( j = 0; j < n; j++ )
-			m_array[i*n+j] = ToReal(mb[j][i]);
-	for ( i = 0, vb = (Num *)b->body; i < n; i++ )
-			b_array[i] = ToReal(vb[i]);
-	trans = 'N';
-	nrhs = 1;
-	lda = n;
-	ldb = n;
-	dgetrf_(&n,&n,m_array,&lda,ipiv,&info);
-	if ( info )
-		error("flinsolve : LU factorization failed");
-	dgetrs_(&trans,&n,&nrhs,m_array,&lda,ipiv,b_array,&ldb,&info);
-	if ( info )
-		error("flinsolve : solving failed");
-	MKVECT(r,n);
-	for ( i = 0, rb = (Real *)r->body; i < n; i++ )
-		MKReal(b_array[i],rb[i]);
-	*rp = r;		
+  mat=(MAT)ARG0(arg); b=(VECT)ARG1(arg);
+  n=mat->row;
+  if ( n != mat->col )
+    error("flinsolve : non-square matrix");
+  if ( n != b->len )
+    error("flinsolve : inconsistent rhs");
+  m_array = (doublereal *)MALLOC_ATOMIC(n*n*sizeof(double));
+  b_array = (doublereal *)MALLOC_ATOMIC(n*sizeof(double));
+  ipiv = (integer *)MALLOC_ATOMIC(n*sizeof(integer));
+  for ( i = 0, mb = (Num **)mat->body; i < n; i++ )
+    for ( j = 0; j < n; j++ )
+      m_array[i*n+j] = ToReal(mb[j][i]);
+  for ( i = 0, vb = (Num *)b->body; i < n; i++ )
+      b_array[i] = ToReal(vb[i]);
+  trans = 'N';
+  nrhs = 1;
+  lda = n;
+  ldb = n;
+  dgetrf_(&n,&n,m_array,&lda,ipiv,&info);
+  if ( info )
+    error("flinsolve : LU factorization failed");
+  dgetrs_(&trans,&n,&nrhs,m_array,&lda,ipiv,b_array,&ldb,&info);
+  if ( info )
+    error("flinsolve : solving failed");
+  MKVECT(r,n);
+  for ( i = 0, rb = (Real *)r->body; i < n; i++ )
+    MKReal(b_array[i],rb[i]);
+  *rp = r;    
 }
 #endif