/[escript]/branches/trilinos_from_5897/paso/src/FCT_Solver.h
ViewVC logotype

Contents of /branches/trilinos_from_5897/paso/src/FCT_Solver.h

Parent Directory Parent Directory | Revision Log Revision Log


Revision 5933 - (show annotations)
Wed Feb 17 23:53:30 2016 UTC (21 months, 4 weeks ago) by caltinay
File MIME type: text/plain
File size: 2028 byte(s)
sync with trunk.

1
2 /*****************************************************************************
3 *
4 * Copyright (c) 2003-2016 by The University of Queensland
5 * http://www.uq.edu.au
6 *
7 * Primary Business: Queensland, Australia
8 * Licensed under the Open Software License version 3.0
9 * http://www.opensource.org/licenses/osl-3.0.php
10 *
11 * Development until 2012 by Earth Systems Science Computational Center (ESSCC)
12 * Development 2012-2013 by School of Earth Sciences
13 * Development from 2014 by Centre for Geoscience Computing (GeoComp)
14 *
15 *****************************************************************************/
16
17 #ifndef __PASO_FCTSOLVER_H__
18 #define __PASO_FCTSOLVER_H__
19
20 #include "Transport.h"
21 #include "FluxLimiter.h"
22 #include "Solver.h"
23
24 namespace paso {
25
26 PASO_DLL_API
27 struct FCT_Solver
28 {
29 FCT_Solver(const_TransportProblem_ptr tp, Options* options);
30
31 ~FCT_Solver();
32
33 err_t update(double* u, double* u_old, Options* options, Performance* pp);
34
35 err_t updateNL(double* u, double* u_old, Options* options, Performance* pp);
36
37 err_t updateLCN(double* u, double* u_old, Options* options, Performance* pp);
38
39 void initialize(double dt, Options* options, Performance* pp);
40
41 static double getSafeTimeStepSize(const_TransportProblem_ptr tp);
42
43 static void setLowOrderOperator(TransportProblem_ptr tp);
44
45 void setAntiDiffusionFlux_linearCN(SystemMatrix_ptr flux_matrix);
46
47 void setAntiDiffusionFlux_BE(SystemMatrix_ptr flux_matrix);
48
49 void setAntiDiffusionFlux_CN(SystemMatrix_ptr flux_matrix);
50
51 void setMuPaLu(double* out, const_Coupler_ptr coupler, double a);
52
53 inline double getTheta()
54 {
55 return method == PASO_BACKWARD_EULER ? 1. : 0.5;
56 }
57
58 const_TransportProblem_ptr transportproblem;
59 esysUtils::JMPI mpi_info;
60 FCT_FluxLimiter* flux_limiter;
61 index_t method;
62 double omega;
63 double dt;
64 double* b;
65 double* z;
66 double* du;
67 Coupler_ptr u_coupler;
68 Coupler_ptr u_old_coupler; /* last time step */
69 };
70
71
72 } // namespace paso
73
74 #endif // __PASO_FCTSOLVER_H__
75

  ViewVC Help
Powered by ViewVC 1.1.26