/*
 *  Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
 *  Copyright (C) 2008-2008 - INRIA - Allan SIMON
 *
 *  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
 *
 */

#include <assert.h>
#include <stdio.h>
#include "matrixEyes.h" 
#include <math.h>

#define LEADDIM 10

static void deyesaTest ( void )
{
    int i = 0 ;
    
    double result = 0 ;
    
    double in [LEADDIM*LEADDIM] ;
    
    deyesa ( in , LEADDIM , LEADDIM ) ;
    
    for ( i = 0 ; i <  LEADDIM*LEADDIM ; i++ )
        {
            if ( i%(LEADDIM+1) == 0 )
                result = 1 ;
            else
                result = 0 ;
            
         printf ( "\t\t %d in : %e\tresult : %e\tassert : %e \n" , i, in[i] , result , fabs( in[i] - result) / fabs( in[i]) ) ;	
		 
         if ( in[i] < 1e-14 && result < 1e-14 )
            assert(1);
         else
            assert ( fabs ( in[i] - result) / fabs( in[i]) < 1e-14 ) ;
        }
}


static void seyesaTest ( void )
{
    int i = 0 ;
    
    float result = 0 ;
    
    float in[LEADDIM*LEADDIM] ;
    
    seyesa ( in , LEADDIM , LEADDIM ) ;
    
    for ( i = 0 ; i <  LEADDIM*LEADDIM ; i++)
        {
            if ( i%(LEADDIM+1) == 0 )
                result = 1.0f ;
            else
                result = 0 ;
            
         printf ( "\t\t %d in : %e\tresult : %e\tassert : %e \n" , i, in[i] , result , fabs( in[i] - result) / fabs( in[i]) ) ;	
         if ( in[i] < 1e-6 && result < 1e-6 )
            assert(1);
         else
		    assert ( fabs ( in[i] - result) / fabs( in[i]) < 1e-6 ) ;
        }
}


static void zeyesaTest ( void )
{
    int i = 0 ;
    
    doubleComplex result = DoubleComplex ( 0 , 0) ;
    
    doubleComplex in[LEADDIM*LEADDIM] ;
    
    zeyesa ( in , LEADDIM , LEADDIM ) ;
    
    for ( i = 0 ; i <  LEADDIM*LEADDIM ; i++ )
        {
            if ( i%(LEADDIM+1) == 0 )
                result = DoubleComplex ( 1 , 0 ) ;
            else
                result = DoubleComplex ( 0 , 0 ) ;
            
         printf ( "\t\t %d in : %e\tresult : %e \n" , i, zreals( in[i]) , zreals ( result ) ) ;
            
         if ( zreals( in[i]) < 1e-14 && zreals( result) < 1e-14 )
            assert(1);
         else   
		    assert ( fabs ( zreals(in[i]) - zreals(result)) / fabs( zreals(in[i])) < 1e-14 ) ;
         
        if ( zimags ( in[i]) < 1e-14 )
            assert (1);
        else 
            assert (0);
         
        }
}





static void ceyesaTest ( void )
{
    int i = 0 ;
    
    floatComplex result = FloatComplex ( 0 , 0  ) ;
    
    floatComplex in[LEADDIM*LEADDIM] ;
    
    ceyesa ( in , LEADDIM , LEADDIM ) ;
    
    for ( i = 0 ; i <  LEADDIM*LEADDIM ; i++)
        {
            if ( i%(LEADDIM+1) == 0 )
                result = FloatComplex ( 1.0f , 0 ) ;
            else
                result = FloatComplex ( 0    , 0 ) ;
            
         printf ( "\t\t %d in : %e\tresult : %e \n" , i, creals( in[i]) , creals ( result ) ) ;
         if ( creals( in[i]) < 1e-6 && creals( result) < 1e-6 )
            assert(1);
         else              
		    assert ( fabs ( creals(in[i]) - creals(result)) / fabs( creals(in[i])) < 1e-6 ) ;
         
        if ( cimags ( in[i]) < 1e-6 )
            assert (1);
        else 
            assert (0);
         
        }
}



static int testExponential(void) {

  printf("\n>>>> Matrix Exponential Tests\n");
  printf("\t>>>> Matrix Double Realt Tests\n");  
  deyesaTest();

  printf("\n\n\t>>>> Matrix Float Realt Tests\n");     
  seyesaTest();
  
  printf("\n\n\n");
  printf("\t>>>> Matrix Float Complex Tests\n");      
  ceyesaTest();  
    
  printf("\n\n\n");
  printf("\t>>>> Matrix Double Complex Tests\n");  
  zeyesaTest();

  return 0;
}



int main(void) {
  assert(testExponential() == 0);
  return 0;
}