cmatrix.h Source File

00001 /*
00002  * cmatrix.h
00003  *
00004  * Copyright (C) 1996 Limit Point Systems, Inc.
00005  *
00006  * Author: Curtis Janssen <cljanss@ca.sandia.gov>
00007  * Maintainer: LPS
00008  *
00009  * This file is part of the SC Toolkit.
00010  *
00011  * The SC Toolkit is free software; you can redistribute it and/or modify
00012  * it under the terms of the GNU Library General Public License as published by
00013  * the Free Software Foundation; either version 2, or (at your option)
00014  * any later version.
00015  *
00016  * The SC Toolkit is distributed in the hope that it will be useful,
00017  * but WITHOUT ANY WARRANTY; without even the implied warranty of
00018  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00019  * GNU Library General Public License for more details.
00020  *
00021  * You should have received a copy of the GNU Library General Public License
00022  * along with the SC Toolkit; see the file COPYING.LIB.  If not, write to
00023  * the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
00024  *
00025  * The U.S. Government is granted a limited license as per AL 91-7.
00026  */
00027 
00028 #ifndef _math_scmat_cmatrix_h
00029 #define _math_scmat_cmatrix_h
00030 
00031 #ifdef __cplusplus
00032 extern "C" {
00033 #endif
00034     
00035 /* These routines work in terms of simple arrays.  No testing is done
00036  * on the input.
00037  */
00038     double cmat_determ(double**matrix,int sym,int dim);
00039     double cmat_invert(double**matrix,int sym,int dim);
00040     double cmat_solve_lin(double**,int sym,double*,int dim);
00041     void cmat_mxm(double**a,int transpose_a,
00042                   double**b,int transpose_b,
00043                   double**c,int transpose_c,
00044                   int nrow, int nlink, int ncol,
00045                   int add);
00046     void cmat_symmetric_mxm(double**a,int na, /* a is (na,na) */
00047                             double**b,int nb, /* b is (na,nb) */
00048                             int add);
00049     void cmat_transform_symmetric_matrix(double**a,int na, /* a is (na,na) */
00050                                          double**b,int nb, /* b is (nb,nb) */
00051                                          double**c,        /* c is (na,nb) */
00052                                          int add);
00053     void cmat_transform_diagonal_matrix(double**a,int na, /* a is (na,na) */
00054                                         double*b,int nb,  /* b is (nb,nb) */
00055                                         double**c,        /* c is (na,nb) */
00056                                         int add);
00057     double** cmat_new_square_matrix(int n);
00058     double** cmat_new_rect_matrix(int n,int m);
00059     void cmat_delete_matrix(double**matrix);
00060     void cmat_transpose_square_matrix(double**matrix,int n);
00061     void cmat_transpose_matrix(double**a,int nrow,int ncol);
00062     void cmat_matrix_pointers(double**ptrs,double*matrix,int nrow, int ncol);
00063     void cmat_diag(double**symm_a, double*evals, double**evecs, int n,
00064                    int matz, double tol);
00065     void cmat_schmidt(double **rows, double *S, int nrow, int nc);
00066     int cmat_schmidt_tol(double **C, double *S, int nrow, int ncol,
00067                          double tolerance, double *res);
00068 
00069 
00070 #ifdef __cplusplus
00071 }
00072 #endif
00073 
00074 #endif