icrout_cholesky_mex.h

00001 /*@HEADER
00002 // ***********************************************************************
00003 // 
00004 //       Ifpack: Object-Oriented Algebraic Preconditioner Package
00005 //                 Copyright (2002) Sandia Corporation
00006 // 
00007 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
00008 // license for use of this work by or on behalf of the U.S. Government.
00009 // 
00010 // This library is free software; you can redistribute it and/or modify
00011 // it under the terms of the GNU Lesser General Public License as
00012 // published by the Free Software Foundation; either version 2.1 of the
00013 // License, or (at your option) any later version.
00014 //  
00015 // This library is distributed in the hope that it will be useful, but
00016 // WITHOUT ANY WARRANTY; without even the implied warranty of
00017 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00018 // Lesser General Public License for more details.
00019 //  
00020 // You should have received a copy of the GNU Lesser General Public
00021 // License along with this library; if not, write to the Free Software
00022 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
00023 // USA
00024 // Questions? Contact Michael A. Heroux (maherou@sandia.gov) 
00025 // 
00026 // ***********************************************************************
00027 //@HEADER
00028 */
00029 
00030 #undef IFPACK
00031 
00032 #ifdef __cplusplus
00033 extern "C" {
00034 #endif
00035 
00036 /*
00037  * Data structure for sparse matrices is CSR, 0-based indexing.
00038  */
00039 typedef struct {
00040     double *val;  /* also known as A  */
00041     int    *col;  /* also known as JA; first column is column 0 */
00042     int    *ptr;  /* also known as IA; with ptr[0] = 0 */
00043 } Matrix;
00044 
00045 void quicksort (int *const pbase, double *const daux, size_t total_elems);
00046 
00047 void Matrix_dealloc(Matrix *a);
00048 void crout_ict(
00049     int n,
00050 #ifdef IFPACK
00051     void * A,
00052     int maxentries,
00053     int (*getcol)( void * A, int col, int ** nentries, double * val, int * ind),
00054     int (*getdiag)( void *A, double * diag),
00055 #else
00056     const Matrix *AL,
00057     const double *Adiag,
00058 #endif
00059     double droptol,
00060     int lfil,
00061     Matrix *L,
00062     double **pdiag);
00063 
00064 #ifdef __cplusplus
00065 }
00066 #endif
00067 

Generated on Tue Oct 20 12:48:53 2009 for IFPACK by doxygen 1.4.7