Changeset 38939

Timestamp:

Sep 1, 2009, 2:00:36 AM (15 years ago)

Author:

hamish

Message:

G_math_solvps() should return something

File:

• grass/branches/develbranch_6/lib/gmath/ccmath_grass_wrapper.c (modified) (15 diffs)

grass/branches/develbranch_6/lib/gmath/ccmath_grass_wrapper.c

@@ -165 +165 @@
     return solv(a[0],b, n);
 }
+
+
 /**
      \brief Solve a symmetric positive definite linear system S*x = b.
@@ -171 +173 @@
      \param  b = array containing system vector b as input and solution vector x as output
      \param  n = dimension of system
-     \return: 0 -> normal exit 1 -> input matrix not positive definite
+     \return: 0 -> normal exit; -1 -> input matrix not positive definite
  */
  int G_math_solvps(double **a,double *b,int n)
 {
-    solvps(a[0], b,n);
-}
+    return solvps(a[0], b,n);
+}
+
+
 /**
      \brief Solve a tridiagonal linear system M*x = y.
@@ -192 +196 @@
     return;
 }
+
+
 /*
      \brief Solve an upper right triangular linear system T*x = b.
@@ -200 +206 @@
      \return value: f = status flag, with 0 -> normal exit, -1 -> system singular
 */
-
 int G_math_solvru(double **a,double *b,int n)
 {
     return solvru(a[0], b, n);
 }
+
+
 /**
      \brief Invert (in place) a general real matrix A -> Inv(A).
-
 
      \param  a = array containing the input matrix A. This is converted to the inverse matrix.
@@ -217 +223 @@
     return minv(a[0], n);
 }
+
+
 /**
      \brief Invert (in place) a symmetric real matrix, V -> Inv(V).
@@ -230 +238 @@
 }
 
-/**
-
+
+/**
      \brief Invert an upper right triangular matrix T -> Inv(T).
-
 
      \param  a = pointer to array of upper right triangular matrix, This is replaced by the inverse matrix.
@@ -243 +250 @@
     return ruinv(a[0], n);
 }
+
+
 /*
 -----------------------------------------------------------------------------
@@ -252 +261 @@
 
      \brief Compute the eigenvalues of a real symmetric matrix A.
-
 
      \param  a = pointer to array of symmetric n by n input matrix A. The computation alters these values.
@@ -263 +271 @@
     return;
 }
+
+
 /**
      \brief Compute the eigenvalues and eigenvectors of a real symmetric matrix A.
@@ -282 +292 @@
     return;
 }
+
+
 /*
      \brief Compute the maximum (absolute) eigenvalue and corresponding eigenvector of a real symmetric matrix A.
@@ -359 +371 @@
 -------------------------------------------------------------------------------
 */
+
+
 /**
      \brief Compute the singular values of a real m by n matrix A.
@@ -374 +388 @@
     return svdval(d, a[0], m, n);
 }
+
+
 /**
 
      \brief Compute singular values when m >> n.
 
-     
      \param  d = pointer to double array of dimension n (output = singular values of A)
      \param  a = pointer to store of the m by n input matrix A (A is altered by the computation)
@@ -389 +404 @@
     return sv2val(d, a[0], m, n);
 }
+
+
 /*
      \brief Compute the singular value transformation S = U~*A*V.
-
      
      \param  d = pointer to double array of dimension n (output = singular values of A)
@@ -406 +422 @@
 }
 
+
 /**
      \brief Compute the singular value transformation when m >> n.
-
      
      \param  d = pointer to double array of dimension n (output = singular values of A)
@@ -422 +438 @@
     return sv2uv(d, a[0], u[0], m, v[0], n);
 }
+
+
 /**
 
      \brief Compute the singular value transformation with A overloaded by the partial U-matrix.
-
      
      \param  d = pointer to double array of dimension n