-HEAD = ../includes/logm.h \

- logm_internal.h

- zlogma.c\

- wbdiag.c \

- cbal.c \

- corth.c \

- cortr.c \

- comqr3.c

-

- $(top_builddir)/matrixOperations/inversion/libMatrixInversion.la\

- $(top_builddir)/elementaryFunctions/sqrt/libSqrt.la \

- $(top_builddir)/auxiliaryFunctions/abs/libAbs.la \

- $(top_builddir)/auxiliaryFunctions/sign/libSign.la \

- $(top_builddir)/auxiliaryFunctions/conj/libConj.la \

- $(top_builddir)/auxiliaryFunctions/pythag/libPythag.la \

- $(top_builddir)/operations/addition/libAddition.la \

- $(top_builddir)/operations/multiplication/libMultiplication.la \

- $(top_builddir)/operations/division/libDivision.la \

- $(top_builddir)/operations/subtraction/libSubtraction.la \

- $(top_builddir)/matrixOperations/transpose/libMatrixTranspose.la\

- $(top_builddir)/matrixOperations/multiplication/libMatrixMultiplication.la\

- $(top_builddir)/elementaryFunctions/pow/libPow.la \

- $(top_builddir)/elementaryFunctions/exp/libExp.la \

- $(top_builddir)/elementaryFunctions/sin/libSin.la \

- $(top_builddir)/elementaryFunctions/cos/libCos.la \

- $(top_builddir)/elementaryFunctions/sinh/libSinh.la \

- $(top_builddir)/elementaryFunctions/cosh/libCosh.la \

-

- libLogm_la-wbdiag.lo libLogm_la-cbal.lo libLogm_la-corth.lo \

- libLogm_la-cortr.lo libLogm_la-comqr3.lo

-HEAD = ../includes/logm.h \

- logm_internal.h

-

- zlogma.c\

- wbdiag.c \

- cbal.c \

- corth.c \

- cortr.c \

- comqr3.c

- $(top_builddir)/matrixOperations/inversion/libMatrixInversion.la\

- $(top_builddir)/elementaryFunctions/sqrt/libSqrt.la \

- $(top_builddir)/auxiliaryFunctions/abs/libAbs.la \

- $(top_builddir)/auxiliaryFunctions/sign/libSign.la \

- $(top_builddir)/auxiliaryFunctions/conj/libConj.la \

- $(top_builddir)/auxiliaryFunctions/pythag/libPythag.la \

- $(top_builddir)/operations/addition/libAddition.la \

- $(top_builddir)/operations/multiplication/libMultiplication.la \

- $(top_builddir)/operations/division/libDivision.la \

- $(top_builddir)/operations/subtraction/libSubtraction.la \

- $(top_builddir)/matrixOperations/transpose/libMatrixTranspose.la\

- $(top_builddir)/matrixOperations/multiplication/libMatrixMultiplication.la\

- $(top_builddir)/elementaryFunctions/pow/libPow.la \

- $(top_builddir)/elementaryFunctions/exp/libExp.la \

- $(top_builddir)/elementaryFunctions/sin/libSin.la \

- $(top_builddir)/elementaryFunctions/cos/libCos.la \

- $(top_builddir)/elementaryFunctions/sinh/libSinh.la \

- $(top_builddir)/elementaryFunctions/cosh/libCosh.la \

-@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libLogm_la-cbal.Plo@am__quote@

-@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libLogm_la-comqr3.Plo@am__quote@

-@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libLogm_la-corth.Plo@am__quote@

-@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libLogm_la-cortr.Plo@am__quote@

-@AMDEP_TRUE@@am__include@ @am__quote@./$(DEPDIR)/libLogm_la-wbdiag.Plo@am__quote@

-libLogm_la-wbdiag.lo: wbdiag.c

-@am__fastdepCC_TRUE@ $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libLogm_la_CFLAGS) $(CFLAGS) -MT libLogm_la-wbdiag.lo -MD -MP -MF $(DEPDIR)/libLogm_la-wbdiag.Tpo -c -o libLogm_la-wbdiag.lo `test -f 'wbdiag.c' || echo '$(srcdir)/'`wbdiag.c

-@am__fastdepCC_TRUE@ mv -f $(DEPDIR)/libLogm_la-wbdiag.Tpo $(DEPDIR)/libLogm_la-wbdiag.Plo

-@AMDEP_TRUE@@am__fastdepCC_FALSE@ source='wbdiag.c' object='libLogm_la-wbdiag.lo' libtool=yes @AMDEPBACKSLASH@

-@AMDEP_TRUE@@am__fastdepCC_FALSE@ DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@

-@am__fastdepCC_FALSE@ $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libLogm_la_CFLAGS) $(CFLAGS) -c -o libLogm_la-wbdiag.lo `test -f 'wbdiag.c' || echo '$(srcdir)/'`wbdiag.c

-

-libLogm_la-cbal.lo: cbal.c

-@am__fastdepCC_TRUE@ $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libLogm_la_CFLAGS) $(CFLAGS) -MT libLogm_la-cbal.lo -MD -MP -MF $(DEPDIR)/libLogm_la-cbal.Tpo -c -o libLogm_la-cbal.lo `test -f 'cbal.c' || echo '$(srcdir)/'`cbal.c

-@am__fastdepCC_TRUE@ mv -f $(DEPDIR)/libLogm_la-cbal.Tpo $(DEPDIR)/libLogm_la-cbal.Plo

-@AMDEP_TRUE@@am__fastdepCC_FALSE@ source='cbal.c' object='libLogm_la-cbal.lo' libtool=yes @AMDEPBACKSLASH@

-@AMDEP_TRUE@@am__fastdepCC_FALSE@ DEPDIR=$(DEPDIR) $(CCDEPMODE) $(depcomp) @AMDEPBACKSLASH@

-@am__fastdepCC_FALSE@ $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libLogm_la_CFLAGS) $(CFLAGS) -c -o libLogm_la-cbal.lo `test -f 'cbal.c' || echo '$(srcdir)/'`cbal.c

-

-libLogm_la-corth.lo: corth.c

-@am__fastdepCC_TRUE@ $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libLogm_la_CFLAGS) $(CFLAGS) -MT libLogm_la-corth.lo -MD -MP -MF $(DEPDIR)/libLogm_la-corth.Tpo -c -o libLogm_la-corth.lo `test -f 'corth.c' || echo '$(srcdir)/'`corth.c

-@am__fastdepCC_TRUE@ mv -f $(DEPDIR)/libLogm_la-corth.Tpo $(DEPDIR)/libLogm_la-corth.Plo

-@AMDEP_TRUE@@am__fastdepCC_FALSE@ source='corth.c' object='libLogm_la-corth.lo' libtool=yes @AMDEPBACKSLASH@

-@am__fastdepCC_FALSE@ $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libLogm_la_CFLAGS) $(CFLAGS) -c -o libLogm_la-corth.lo `test -f 'corth.c' || echo '$(srcdir)/'`corth.c

-libLogm_la-cortr.lo: cortr.c

-@am__fastdepCC_TRUE@ $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libLogm_la_CFLAGS) $(CFLAGS) -MT libLogm_la-cortr.lo -MD -MP -MF $(DEPDIR)/libLogm_la-cortr.Tpo -c -o libLogm_la-cortr.lo `test -f 'cortr.c' || echo '$(srcdir)/'`cortr.c

-@am__fastdepCC_TRUE@ mv -f $(DEPDIR)/libLogm_la-cortr.Tpo $(DEPDIR)/libLogm_la-cortr.Plo

-@AMDEP_TRUE@@am__fastdepCC_FALSE@ source='cortr.c' object='libLogm_la-cortr.lo' libtool=yes @AMDEPBACKSLASH@

-@am__fastdepCC_FALSE@ $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libLogm_la_CFLAGS) $(CFLAGS) -c -o libLogm_la-cortr.lo `test -f 'cortr.c' || echo '$(srcdir)/'`cortr.c

-libLogm_la-comqr3.lo: comqr3.c

-@am__fastdepCC_TRUE@ $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libLogm_la_CFLAGS) $(CFLAGS) -MT libLogm_la-comqr3.lo -MD -MP -MF $(DEPDIR)/libLogm_la-comqr3.Tpo -c -o libLogm_la-comqr3.lo `test -f 'comqr3.c' || echo '$(srcdir)/'`comqr3.c

-@am__fastdepCC_TRUE@ mv -f $(DEPDIR)/libLogm_la-comqr3.Tpo $(DEPDIR)/libLogm_la-comqr3.Plo

-@AMDEP_TRUE@@am__fastdepCC_FALSE@ source='comqr3.c' object='libLogm_la-comqr3.lo' libtool=yes @AMDEPBACKSLASH@

-@am__fastdepCC_FALSE@ $(LIBTOOL) --tag=CC $(AM_LIBTOOLFLAGS) $(LIBTOOLFLAGS) --mode=compile $(CC) $(DEFS) $(DEFAULT_INCLUDES) $(INCLUDES) $(AM_CPPFLAGS) $(CPPFLAGS) $(libLogm_la_CFLAGS) $(CFLAGS) -c -o libLogm_la-comqr3.lo `test -f 'comqr3.c' || echo '$(srcdir)/'`comqr3.c

-#include "schur.h"

-#include "log.h"

-#include "multiplication.h"

-#include "tranpose.h"

-#include "inverse.h"

-

-static void bdiag(floatComplex *in, int size, double n, floatComplex *out1, floatComplex *out2){

-

-}

-

- doubleComplex *out1, *out2, *dtmp;

- int i=0, egaux=0;

- float eps = 2.220e-16;

- int select, sdim, lwork, bwork, info = 0;

- doubleComplex *eigenvalues, *work, *rwork;

-

-

- out1=malloc((uint)(size*size)*sizeof(doubleComplex));

- out2=malloc((uint)(size*size)*sizeof(doubleComplex));

- dtmp=malloc((uint)(size*size)*sizeof(doubleComplex));

- eigenvalues=malloc((uint)size*sizeof(doubleComplex));

- work=malloc((uint)(pow(2,size))*sizeof(doubleComplex));

- rwork=malloc((uint)size*sizeof(doubleComplex));

-

- ctransposea(in,size,size,tmp);/*pb entre float et double*/

-

- for (i=0;i<size*size;i++){

- if (in[i]!=tmp[i]) break;

- }

- if (i==size*size) egaux =1;

-

-

- if (egaux){

- /*Hermitian matrix*/

- for (i=0;i<size*size,i++)

- out1[i]=DoubleComplex((double)creals(in[i]),(double)cimags(in[i]));;

-

- C2F(zgees)("V","S", &select, &size, out1, &size, &sdim, eigenvalues, out2, &size, work, &lwork, rwork, &bwork, &info);

-

- for (i=0;i<size;i++){

- for (j=0;j<size;j++){

- tmp[i*size+j]=0;

- if (i==j) dtmp[i*size+j]=zlogs(out2[i*size+j]);

- }

- }

-

- zmulma(out1,size,size,dtmp,size,size,out2);

- ztransposea(out1,size,size,dtmp);

- cmula(out2,size,size,tmp,size,size,out1);

- for (i=0;i<size*size,i++)

- out[i]=FloatComplex((float)zreals(in[i]),(float)zimags(in[i]));;

- }

- else

- {

- /*General Matrix*/

- bdiag(in, size, eps, out1, out2);

- for (i=0;i<size;i++){

- for (j=0;j<size;j++){

- tmp[i*size+j]=0;

- if (i==j) tmp[i*size+j]=clogs(out2[i*size+j]);

- }

- }

- out2 = cmulma(out1,tmp);

- tmp = cinversea(out1,size,size,out1);

- out= cmula(out2,tmp);

- }

- floatComplex *inCpx;

- inCpx=malloc((uint)(size*size)*sizeof(floatComplex));

- for (i=0;i<size*size) inCpx[i] = FloatComplex(in[i],0);

- clogma(inCpx, size, out);

- wbdiag(in,3,out);

-

-#include "pow.h"

-#include "logm_internal.h"

- doubleComplex *out1, *out2, *tmp, *inCopy;

- int i=0, j=0, egaux=0;

- int mon_select, sdim, lwork, bwork, info = 0;

- doubleComplex *eigenvalues, *work, *rwork;

- inCopy=malloc((uint)(size*size)*sizeof(doubleComplex));

- out1=malloc((uint)(size*size)*sizeof(doubleComplex));

- out2=malloc((uint)(size*size)*sizeof(doubleComplex));

- tmp=malloc((uint)(size*size)*sizeof(doubleComplex));

- eigenvalues=malloc((uint)size*sizeof(doubleComplex));

- work=malloc((uint)(dpows(2,size))*sizeof(doubleComplex));

- rwork=malloc((uint)size*sizeof(doubleComplex));

- for (i=0;i<size*size;i++) inCopy[i]=in[i];

- ztransposea(in,size,size,tmp);

- /* Check if in and transpose(in) are equals */

- for (i=0;i<size*size;i++){

- if ( (zreals(in[i])!=zreals(tmp[i])) || (zimags(in[i])!=zimags(tmp[i])) ) break;

- }

- if (i==size*size) egaux =1;

- for (i=0;i<size*size;i++) out1[i]=in[i];

- if (egaux){

- /*Hermitian matrix*/

-

-

- C2F(zgees)("V","S", &mon_select, &size, out1, &size, &sdim, eigenvalues, out2, &size, work, &lwork, rwork, &bwork, &info);

- for (i=0;i<size;i++){

- for (j=0;j<size;j++){

- tmp[i*size+j]=DoubleComplex(0,0);

- if (i==j) tmp[i*size+j]=zlogs(out2[i*size+j]);

- }

- }

- zmulma(out1,size,size,tmp,size,size,out2);

- ztransposea(out1,size,size,tmp);

-

- else

- {

- /*General Matrix*/

- /* A faire

- wbdiag (appel à bdiag dans le script mais le bdiag de scilab

- appelle sci_bdiag qui lui appelle wbdiag car input complexe)

- */

- wbdiag(out1,size,out2);

- for (i=0;i<size;i++){

- for (j=0;j<size;j++){

- tmp[i*size+j]=DoubleComplex(0,0);

- if (i==j) tmp[i*size+j]=zlogs(out2[i*size+j]);

- }

- zmulma(out1,size,size,tmp,size,size,out2);

- zinverma(out1,tmp,size);

-

- zmulma(out2,size,size,tmp,size,size,out);