/[escript]/trunk/paso/src/SystemMatrix_getSystemMatrixTypeId.c
ViewVC logotype

Annotation of /trunk/paso/src/SystemMatrix_getSystemMatrixTypeId.c

Parent Directory Parent Directory | Revision Log Revision Log


Revision 633 - (hide annotations)
Thu Mar 23 05:37:00 2006 UTC (13 years, 6 months ago) by dhawcroft
File MIME type: text/plain
File size: 2302 byte(s)


1 jgs 150 /* $Id$ */
2    
3 dhawcroft 631 /*
4     ********************************************************************************
5 dhawcroft 633 * Copyright 2006 by ACcESS MNRF *
6 dhawcroft 631 * *
7     * http://www.access.edu.au *
8     * Primary Business: Queensland, Australia *
9     * Licensed under the Open Software License version 3.0 *
10     * http://www.opensource.org/licenses/osl-3.0.php *
11     ********************************************************************************
12     */
13 jgs 150 /**************************************************************/
14    
15     /* Paso: returns the matrix format requested by a particular linear solver */
16    
17     /**************************************************************/
18    
19     /* Copyrights by ACcESS Australia 2004,2005 */
20     /* Author: gross@access.edu.au */
21    
22     /**************************************************************/
23    
24     #include "Paso.h"
25     #include "Options.h"
26     #include "SystemMatrix.h"
27    
28     /**************************************************************/
29    
30     index_t Paso_SystemMatrix_getSystemMatrixTypeId(index_t solver,index_t package, bool_t symmetry) {
31 gross 415 index_t out=MATRIX_FORMAT_DEFAULT;
32 jgs 150 package=Paso_Options_getPackage(solver,package,symmetry);
33    
34     switch(package) {
35    
36     case PASO_PASO:
37 gross 415 out=MATRIX_FORMAT_DEFAULT;
38 jgs 150 break;
39    
40     case PASO_SCSL:
41 gross 415 out=MATRIX_FORMAT_CSC + MATRIX_FORMAT_BLK1;
42 gross 425 /* if (solver == PASO_CHOLEVSKY) out+=MATRIX_FORMAT_SYM */
43 jgs 150 break;
44 gross 415
45 jgs 150 case PASO_MKL:
46 gross 425 out=MATRIX_FORMAT_BLK1 + MATRIX_FORMAT_OFFSET1;
47     /* if (solver == PASO_CHOLEVSKY) out+=MATRIX_FORMAT_SYM */
48 jgs 150 break;
49    
50     case PASO_UMFPACK:
51 gross 415 out=MATRIX_FORMAT_BLK1;
52 jgs 150 break;
53    
54     default:
55     Paso_setError(VALUE_ERROR,"unknown package code");
56     }
57     return out;
58     }
59     /*
60     * $Log$
61     * Revision 1.2 2005/09/15 03:44:39 jgs
62     * Merge of development branch dev-02 back to main trunk on 2005-09-15
63     *
64     * Revision 1.1.2.1 2005/09/05 06:29:48 gross
65     * These files have been extracted from finley to define a stand alone libray for iterative
66     * linear solvers on the ALTIX. main entry through Paso_solve. this version compiles but
67     * has not been tested yet.
68     *
69     *
70     */

Properties

Name Value
svn:eol-style native
svn:keywords Author Date Id Revision

  ViewVC Help
Powered by ViewVC 1.1.26