/[escript]/trunk/paso/src/SolverFCT.h
ViewVC logotype

Contents of /trunk/paso/src/SolverFCT.h

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1366 - (show annotations)
Tue Dec 18 05:49:17 2007 UTC (12 years, 9 months ago) by gross
File MIME type: text/plain
File size: 1950 byte(s)
finley can create a TransportProblem now
1 /* $Id: $ */
2
3 /*******************************************************
4 *
5 * Copyright 2007 by University of Queensland
6 *
7 * http://esscc.uq.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
14 #ifndef INC_SOLVERFCT
15 #define INC_SOLVERFCT
16
17 #include "SystemMatrix.h"
18
19 typedef struct Paso_FCTransportProblem {
20
21 double theta;
22 bool_t valid_matrices;
23
24 double* u;
25
26 Paso_SystemMatrix * transport_matrix;
27 Paso_SystemMatrix * flux_matrix;
28 double* lumped_mass_matrix;
29 double* row_sum_flux_matrix;
30
31 dim_t num_colors;
32 index_t *colorOf;
33 index_t *main_iptr;
34
35 Paso_MPIInfo *mpi_info;
36 dim_t reference_counter;
37
38 } Paso_FCTransportProblem;
39
40 void Paso_FCTransportProblem_free(Paso_FCTransportProblem* in);
41 Paso_FCTransportProblem* Paso_FCTransportProblem_getReference(Paso_FCTransportProblem* in);
42 Paso_SystemMatrix* Paso_FCTransportProblem_borrowTransportMatrix(Paso_FCTransportProblem* in);
43 Paso_SystemMatrix* Paso_FCTransportProblem_borrowFluxMatrix(Paso_FCTransportProblem* in);
44 double* Paso_FCTransportProblem_borrowLumpedMassMatrix(Paso_FCTransportProblem* in);
45 dim_t Paso_FCTransportProblem_getTotalNumRows(Paso_FCTransportProblem* in);
46 Paso_FCTransportProblem* Paso_FCTransportProblem_alloc(double theta, double dt_max, Paso_SystemMatrixPattern *pattern, int block_size);
47 void Paso_FCTransportProblem_setAntiDiffusiveFlux(Paso_FCTransportProblem * fc, double * u, double* fa);
48 void Paso_FCTransportProblem_addAdvectivePart(Paso_FCTransportProblem * fc, double alpha);
49 void Paso_SolverFCT_solve(Paso_FCTransportProblem* fctp, double* u, double dt, double* source, Paso_Options* options);
50 void Paso_FCTransportProblem_checkinSolution(Paso_FCTransportProblem* in, double* u) ;
51
52
53 #endif /* #ifndef INC_SOLVERFCT */

  ViewVC Help
Powered by ViewVC 1.1.26