diff options
| -rw-r--r-- | src/matrixOperations/includes/spec.h | 48 |
1 files changed, 48 insertions, 0 deletions
diff --git a/src/matrixOperations/includes/spec.h b/src/matrixOperations/includes/spec.h new file mode 100644 index 00000000..0ab5e4b6 --- /dev/null +++ b/src/matrixOperations/includes/spec.h @@ -0,0 +1,48 @@ +/* + * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab + * Copyright (C) 2008 - INRIA - Arnaud TORSET + * + * This file must be used under the terms of the CeCILL. + * This source file is licensed as described in the file COPYING, which + * you should have received as part of this distribution. The terms + * are also available at + * http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt + * + */ + + +#ifndef __SPEC_H__ +#define __SPEC_H__ + +#include "doubleComplex.h" +#include "floatComplex.h" + +/* spec gives only the eigenvalues + If you want the eigenvalues and the eigenvectors, use spec2 */ + +/* spec */ +#define sspecs(in) in +#define dspecs(in) in +#define cspecs(in) in +#define zspecs(in) in + +void sspeca(float* in, int rows, floatComplex* out); +void dspeca(double* in, int rows, doubleComplex* out); +void cspeca(floatComplex* in, int rows, floatComplex* out); +void zspeca(doubleComplex* in, int rows,doubleComplex* out); + + + +/* spec2 */ +#define sspec2s(in,out) {out=in;return 1;} +#define dspec2s(in,out) {out=in;return 1;} +#define cspec2s(in,out) {out=in;return FloatComplex(1,0);} +#define zspec2s(in,out) {out=in;return DoubleCompelx(1,0);} + +void sspec2a(float* in, int rows, floatComplex* eigenvalues,floatComplex* eigenvectors); +void dspec2a(double* in, int rows, doubleComplex* eigenvalues,doubleComplex* eigenvectors); +void cspec2a(floatComplex* in, int rows, floatComplex* eigenvalues,floatComplex* eigenvectors); +void zspec2a(doubleComplex* in, int rows,doubleComplex* eigenvalues,doubleComplex* eigenvectors); + + +#endif /* __SPEC_H__ */ |