diff options
| author | Brijeshcr | 2017-07-07 17:09:19 +0530 |
|---|---|---|
| committer | Brijeshcr | 2017-07-07 17:09:19 +0530 |
| commit | aa250753b439916e4df235839e9266f058632f80 (patch) | |
| tree | 75f7c141ed31ba0c2291d744beeb8df655b1288d /src/c/linearAlgebra/hess | |
| parent | 124ef41ebab756797803b30c4c07ce37213a6477 (diff) | |
| parent | 45bb4c81ba70d5f711d8a5ea207ac18b141d0cb2 (diff) | |
| download | scilab2c-aa250753b439916e4df235839e9266f058632f80.tar.gz scilab2c-aa250753b439916e4df235839e9266f058632f80.tar.bz2 scilab2c-aa250753b439916e4df235839e9266f058632f80.zip | |
Linear Algebra changes by Sandeep
Diffstat (limited to 'src/c/linearAlgebra/hess')
| -rw-r--r-- | src/c/linearAlgebra/hess/dhessa.c | 8 |
1 files changed, 5 insertions, 3 deletions
diff --git a/src/c/linearAlgebra/hess/dhessa.c b/src/c/linearAlgebra/hess/dhessa.c index 57f81b35..e1f2e2d1 100644 --- a/src/c/linearAlgebra/hess/dhessa.c +++ b/src/c/linearAlgebra/hess/dhessa.c @@ -20,11 +20,13 @@ #include "matrixTranspose.h" #include "matrixMultiplication.h" +/* Lapack subroutines - which are used*/ extern int dgehrd_(int *, int *,int *,double *,int *,double *,double *,int *,int *); extern int dorghr_(int *, int *,int *,double *,int *,double *,double *,int *,int *); - +/* All the vairbale names are given exactly the same name as scilab source code */ void dhessa(double *in1,int size,int nout,double *out1, double *out2){ +/* Variables names are done through, Lapack library. */ int i,j,k; int N = size; int ILO=1; @@ -41,11 +43,11 @@ void dhessa(double *in1,int size,int nout,double *out1, double *out2){ WORK = (double *)malloc((double)LWORK*sizeof(double)); dgehrd_(&N,&ILO,&IHI,A,&LDA,TAU,WORK,&N,&INFO); - for(i=0;i<N;i++) + for(i=0;i<N;i++) /* copying it in output */ for(j=0;j<N;j++) out2[i+j*N] = A[i+j*N]; - for(j=1;j<=N-2;j++){ + for(j=1;j<=N-2;j++){ /* copying it in output */ for(i=j+2;i<=N;i++){ out2[(i-1)+(j-1)*N] = 0; } |