diff options
| author | siddhu8990 | 2017-02-07 16:16:31 +0530 |
|---|---|---|
| committer | siddhu8990 | 2017-02-07 16:16:31 +0530 |
| commit | e7dba33a909e158dffc6ddb3361df10d427db631 (patch) | |
| tree | d71a27ff842d6088798ce8a0db16decb1fe10853 /src/c/linearAlgebra/schur/dgschura.c | |
| parent | 765d9c44f94634406eeff50e20e8cdfcf1b7699c (diff) | |
| download | scilab2c-e7dba33a909e158dffc6ddb3361df10d427db631.tar.gz scilab2c-e7dba33a909e158dffc6ddb3361df10d427db631.tar.bz2 scilab2c-e7dba33a909e158dffc6ddb3361df10d427db631.zip | |
Support for 'lqr' and 'lqe' added
Diffstat (limited to 'src/c/linearAlgebra/schur/dgschura.c')
| -rw-r--r-- | src/c/linearAlgebra/schur/dgschura.c | 17 |
1 files changed, 13 insertions, 4 deletions
diff --git a/src/c/linearAlgebra/schur/dgschura.c b/src/c/linearAlgebra/schur/dgschura.c index 56d48455..f17da8c8 100644 --- a/src/c/linearAlgebra/schur/dgschura.c +++ b/src/c/linearAlgebra/schur/dgschura.c @@ -15,14 +15,15 @@ #include "lapack.h" #include "stdlib.h" #include "string.h" - +#include <math.h> /*flag --> 0: nothing --> 1: continuous --> 2: discrete */ -lapack_logical selctg2( double* in1, double* in2, double* in3); +lapack_logical selctg21( double* in1, double* in2, double* in3); +lapack_logical selctg22( double* in1, double* in2, double* in3); double dgschura(double* in1, int size, double* in2, int flag, int nout, \ double* out1, double* out2, double* out3, double* out4) @@ -40,7 +41,7 @@ double dgschura(double* in1, int size, double* in2, int flag, int nout, \ modifies the input variable*/ /*Used incase of flag > 0*/ - LAPACK_D_SELECT3 selctg = &selctg2; + LAPACK_D_SELECT3 selctg = &selctg21; if(nout >= 2){ JOBVSL = 'V'; @@ -143,7 +144,15 @@ double dgschura(double* in1, int size, double* in2, int flag, int nout, \ return ret; } -lapack_logical selctg2(double* in1, double* in2, double* in3) +lapack_logical selctg21(double* in1, double* in2, double* in3) +{ + if((sqrt(*in1**in1+*in2**in2)/ *in3) < 1) + return 1; + else + return 0; +} + +lapack_logical selctg22(double* in1, double* in2, double* in3) { if((sqrt(*in1**in1+*in2**in2)/ *in3) < 1) return 1; |