summaryrefslogtreecommitdiff
path: root/src/c/linearAlgebra/sva
diff options
context:
space:
mode:
authorSunil Shetye2019-05-17 17:30:07 +0530
committerSunil Shetye2019-05-29 11:08:06 +0530
commit716c7a821a76e65bff30d0292095c585e37c3ac0 (patch)
treee56be10d77d70990c749b139d063b6bda69523ec /src/c/linearAlgebra/sva
parent953447ba5cd043470d0f35c8a780e1b4c4d0aebf (diff)
downloadscilab2c-716c7a821a76e65bff30d0292095c585e37c3ac0.tar.gz
scilab2c-716c7a821a76e65bff30d0292095c585e37c3ac0.tar.bz2
scilab2c-716c7a821a76e65bff30d0292095c585e37c3ac0.zip
fix compilation warningsHEADmaster
Diffstat (limited to 'src/c/linearAlgebra/sva')
-rw-r--r--src/c/linearAlgebra/sva/dsvaa.c17
1 files changed, 8 insertions, 9 deletions
diff --git a/src/c/linearAlgebra/sva/dsvaa.c b/src/c/linearAlgebra/sva/dsvaa.c
index 691694e4..6c9d95e1 100644
--- a/src/c/linearAlgebra/sva/dsvaa.c
+++ b/src/c/linearAlgebra/sva/dsvaa.c
@@ -17,6 +17,8 @@
#include <stdlib.h>
#include "string.h"
#include "matrixTranspose.h"
+#include "max.h"
+#include "min.h"
#define eps 2.22044604925e-16
@@ -30,18 +32,17 @@ void dsvaa(int ninp,double *in1,int row,int col,double in2,double *out1, \
if(row == 0 && col == 0) return;
int i,j;
int arow; /*Actual row of given matrix*/
- int acol; /*Actual col of given matrix*/
/* Calculation of svd of a given matrix */
double *U,*S,*V;
- U = (double *)malloc((double)row*Min(row,col)*sizeof(double));
- S = (double *)malloc((double)Min(row,col)*Min(row,col)*sizeof(double));
- V = (double *)malloc((double)col*Min(row,col)*sizeof(double));
+ U = (double *)malloc((double)row*min(row,col)*sizeof(double));
+ S = (double *)malloc((double)min(row,col)*min(row,col)*sizeof(double));
+ V = (double *)malloc((double)col*min(row,col)*sizeof(double));
dsvda(0,in1,M,N,1,3,U,S,V);
if (ninp == 1){ /* [u,s,v] = sva(A) when input is only matrix */
- tol = Max(row,col)*S[0]*eps;
+ tol = max(row,col)*S[0]*eps;
rk = 0;
for(i=0;i<col;i++){
if(S[i+i*row] > tol){
@@ -53,7 +54,7 @@ void dsvaa(int ninp,double *in1,int row,int col,double in2,double *out1, \
tol = in2;
if(tol > 1){
rk = tol;
- if(rk > Min(row,col)){
+ if(rk > min(row,col)){
printf("ERROR: Wrong value for input argument !");
out1 = NULL;
out2 = NULL;
@@ -71,20 +72,18 @@ void dsvaa(int ninp,double *in1,int row,int col,double in2,double *out1, \
}
}
arow = M;
- acol = Min(M,N); /* Copying, the output in required format */
for(i=0;i<arow;i++){
for(j=0;j<rk;j++){
out1[i+j*row]=U[i+j*arow];
}
}
- arow = Min(M,N);
+ arow = min(M,N);
for(i=0;i<rk;i++){ /* Copying, the output in required format */
for(j=0;j<rk;j++){
out2[i+j*(int)rk] = S[i+j*arow];
}
}
arow = N;
- acol = Min(M,N);
for(i=0;i<arow;i++){ /* Copying, the output in required format */
for(j=0;j<rk;j++){
out3[i+j*arow] = V[i+j*arow];