/[escript]/branches/RW_WIN32/finley/src/CPPAdapter/SystemMatrixAdapter.cpp
ViewVC logotype

Annotation of /branches/RW_WIN32/finley/src/CPPAdapter/SystemMatrixAdapter.cpp

Parent Directory Parent Directory | Revision Log Revision Log


Revision 100 - (hide annotations)
Wed Dec 15 03:48:48 2004 UTC (16 years, 2 months ago) by jgs
Original Path: trunk/esys2/finley/src/CPPAdapter/SystemMatrixAdapter.cpp
File size: 3933 byte(s)
*** empty log message ***

1 jgs 82 /*
2     ******************************************************************************
3     * *
4     * COPYRIGHT ACcESS 2004 - All Rights Reserved *
5     * *
6     * This software is the property of ACcESS. No part of this code *
7     * may be copied in any form or by any means without the expressed written *
8     * consent of ACcESS. Copying, use or modification of this software *
9     * by any unauthorised person is illegal unless that person has a software *
10     * license agreement with ACcESS. *
11     * *
12     ******************************************************************************
13     */
14     extern "C" {
15     #include "finley/finleyC/System.h"
16     }
17     #include "escript/Data/Data.h"
18     #include "finley/CPPAdapter/SystemMatrixAdapter.h"
19     #include "finley/CPPAdapter/FinleyAdapterException.h"
20     #include "finley/CPPAdapter/FinleyError.h"
21     #include <boost/python/extract.hpp>
22    
23     using namespace std;
24    
25     namespace finley {
26    
27     SystemMatrixAdapter::SystemMatrixAdapter()
28     {
29     throw FinleyAdapterException("Error - Illegal to generate default SystemMatrixAdapter.");
30     }
31    
32 jgs 100 SystemMatrixAdapter::SystemMatrixAdapter(const Finley_SystemMatrix* system_matrix,
33 jgs 82 const int row_blocksize,
34     const escript::FunctionSpace& row_functionspace,
35     const int column_blocksize,
36     const escript::FunctionSpace& column_functionspace):
37 jgs 100 AbstractSystemMatrix(row_blocksize,row_functionspace,column_blocksize,column_functionspace),
38     m_system_matrix(system_matrix)
39 jgs 82 {
40     }
41    
42     SystemMatrixAdapter::~SystemMatrixAdapter()
43     {
44     if (m_system_matrix.unique()) {
45     Finley_SystemMatrix* mat=m_system_matrix.get();
46     Finley_SystemMatrix_dealloc(mat);
47     }
48     }
49    
50     Finley_SystemMatrix* SystemMatrixAdapter::getFinley_SystemMatrix() const
51     {
52     return m_system_matrix.get();
53     }
54    
55     void SystemMatrixAdapter::ypAx(escript::Data& y, const escript::Data& x) const
56     {
57     Finley_SystemMatrix* mat=getFinley_SystemMatrix();
58     Finley_SystemMatrixVector(&(y.getDataC()),mat,&(x.getDataC()));
59     checkFinleyError();
60     }
61    
62     void SystemMatrixAdapter::setToSolution(escript::Data& out, const escript::Data& in, const boost::python::dict& options) const
63     {
64     Finley_SolverOptions finley_options;
65     Finley_SystemMatrix_setDefaults(&finley_options);
66     // extract options
67     #define EXTRACT(__key__,__val__,__type__) if ( options.has_key(__key__)) finley_options.__val__=boost::python::extract<__type__>(options.get(__key__))
68     EXTRACT("verbose",verbose,int);
69     EXTRACT("reordering",reordering,int);
70 jgs 100 EXTRACT("tolerance",tolerance,double);
71     EXTRACT("iterative_method",iterative_method,int);
72 jgs 82 EXTRACT("preconditioner",preconditioner,int);
73     EXTRACT("iter_max",iter_max,int);
74     EXTRACT("drop_tolerance",drop_tolerance,double);
75     EXTRACT("drop_storage",drop_storage,double);
76 jgs 100 EXTRACT("iterative",iterative,int);
77 jgs 82 #undef EXTRACT
78 jgs 100 if (finley_options.iterative) {
79     Finley_SystemMatrix_iterative(getFinley_SystemMatrix(),&(out.getDataC()),&(in.getDataC()),&finley_options);
80     } else {
81     Finley_SystemMatrix_solve(getFinley_SystemMatrix(),&(out.getDataC()),&(in.getDataC()),&finley_options);
82     }
83 jgs 82 checkFinleyError();
84     }
85    
86     void SystemMatrixAdapter::nullifyRowsAndCols(const escript::Data& row_q, const escript::Data& col_q, const double mdv) const
87     {
88     Finley_SystemMatrix* system_matrix_ptr = getFinley_SystemMatrix();
89     escriptDataC row_qC = row_q.getDataC();
90     escriptDataC col_qC = col_q.getDataC();
91    
92     Finley_SystemMatrix_nullifyRowsAndCols(system_matrix_ptr, &row_qC, &col_qC, mdv);
93     }
94    
95     } // end of namespace

Properties

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

  ViewVC Help
Powered by ViewVC 1.1.26