#ifndef _Eigen_h_
#define _Eigen_h_
/* Eigen.h
 *
 * Copyright (C) 1993-2002 David Weenink
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or (at
 * your option) any later version.
 *
 * This program is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
 */

/*
 djmw 20020813 GPL header
*/

#ifndef _Data_h_
	#include "Data.h"
#endif
#ifndef _Graphics_h_
	#include "Graphics.h"
#endif

#define Eigen_createFromSymmetricMatrix Eigen_createFromSymmetricMatrix_d
#define Eigen_initFromSymmetricMatrix Eigen_initFromSymmetricMatrix_d

#define Eigen_members Data_members \
	long numberOfEigenvalues, dimension;	\
	double *eigenvalues;	\
	double **eigenvectors; /* eigenvectors stored in row */
#define Eigen_methods Data_methods
class_create (Eigen, Data)

Eigen Eigen_create (long numberOfEigenvalues, long dimension);

Eigen Eigen_createFromSymmetricMatrix_f (float **a, long n);
Eigen Eigen_createFromSymmetricMatrix_d (double **a, long n);

int Eigen_init (I, long numberOfEigenvalues, long dimension);

int Eigen_initFromSymmetricMatrix_f (I, float **a, long n);
int Eigen_initFromSymmetricMatrix_d (I, double **a, long n);

int Eigen_initFromSquareRoot (I, double **a, long numberOfRows,
	long numberOfColumns);
/*
	Calculate eigenstructure for symmetric matrix A'A (e.g. covariance matrix),
	when only A is given.
	Precondition: numberOfRows > 1
	Method: SVD.
*/

int Eigen_initFromSquareRootPair (I, double **a, long numberOfRows,
	long numberOfColumns, double **b, long numberOfRows_b);
/*
	Calculate eigenstructure for A'Ax - lambda B'Bx = 0
	Preconditions: numberOfRows >= numberOfColumns &&
		numberOfRows_b >= numberOfColumns
	Method: Generalized SVD.
*/

long Eigen_getNumberOfEigenvectors (I);

long Eigen_getDimensionOfComponents (I);

double Eigen_getCumulativeContributionOfComponents (I, long from, long to);

long Eigen_getDimensionOfFraction (I, double fraction);

double Eigen_getSumOfEigenvalues (I);

void Eigen_sort (I);
/*
	Sort eigenvalues and corresponding eigenvectors in decreasing order.
*/

void Eigen_invertEigenvector (I, long ivec);

void Eigen_drawEigenvalues_scree (I, Graphics g, long first, long last, double minimum,
	double maximum, double size_mm, const char *mark, int garnish);
	
void Eigen_drawEigenvector (I, Graphics g, long ivec, long first, long last,
	double minimum, double maximum, int weigh, double size_mm, const char *mark,
	int connect, int garnish);

#endif /* _Eigen_h_ */

/* End of file Eigen.h */