1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
|
/* This file is part of Charon.
Charon is free software: you can redistribute it and/or modify
it under the terms of the GNU Lesser General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
Charon is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public License
along with Charon. If not, see <http://www.gnu.org/licenses/>.
*/
/** \file petsc4_PetscSharedLib.cpp
* Library to perform simple calls to petsc.
* \author <a href="Jens-Malte.Gottfried@iwr.uni-heidelberg.de">
* Jens-Malte Gottfried </a>
* \date 25.10.2010
*/
#include <petscksp.h>
#include <iostream>
/// do something arbitrary with petsc
extern "C" {
int run(int& argc, char**& argv) {
// testing Petsc Initialization
std::cout << "entering PetscTestDLL" << std::endl;
PetscErrorCode ierr = MPI_Init(&argc,&argv);
CHKERRQ(ierr);
ierr = PetscInitialize(&argc,&argv,PETSC_NULL,PETSC_NULL);
CHKERRQ(ierr);
Vec x;
ierr = VecCreate(PETSC_COMM_WORLD,&x);
CHKERRQ(ierr);
ierr = PetscFinalize();
CHKERRQ(ierr);
int initialized = 0;
MPI_Initialized(&initialized);
if (initialized) {
ierr = MPI_Finalize();
CHKERRQ(ierr);
}
std::cout << "leaving PetscTestDLL" << std::endl;
return 0;
}
}
|