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
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
|
<center><a href="ex1.c">Actual source code: ex1.c</a></center><br>
<html>
<head> <link rel="canonical" href="http://www.mcs.anl.gov/petsc/petsc-current/tutorials/python/ex1.c.html" />
<title></title>
<meta name="generator" content="c2html 0.9.4">
<meta name="date" content="2017-01-01T16:37:10+00:00">
</head>
<body bgcolor="#FFFFFF">
<div id="version" align=right><b>petsc-3.7.5 2017-01-01</b></div>
<div id="bugreport" align=right><a href="mailto:petsc-maint@mcs.anl.gov?subject=Typo or Error in Documentation &body=Please describe the typo or error in the documentation: petsc-3.7.5 v3.7.5 tutorials/python/ex1.c.html "><small>Report Typos and Errors</small></a></div>
<pre width="80"><a name="line1"> 1: </a><font color="#A020F0">#include <petscksp.h></font>
<a name="line3"> 3: </a><font color="#B22222">/* ------------------------------------------------------- */</font>
<a name="line7"> 7: </a><strong><font color="#4169E1"><a name="RunTest"></a><a href="tutorials/../../docs/manualpages/Sys/PetscErrorCode.html#PetscErrorCode">PetscErrorCode</a> RunTest(void)</font></strong>
<a name="line8"> 8: </a>{
<a name="line9"> 9: </a> <a href="tutorials/../../docs/manualpages/Sys/PetscInt.html#PetscInt">PetscInt</a> N = 100;
<a name="line10"> 10: </a> <a href="tutorials/../../docs/manualpages/Sys/PetscBool.html#PetscBool">PetscBool</a> draw = <a href="tutorials/../../docs/manualpages/Sys/PETSC_FALSE.html#PETSC_FALSE">PETSC_FALSE</a>;
<a name="line11"> 11: </a> <a href="tutorials/../../docs/manualpages/Sys/PetscReal.html#PetscReal">PetscReal</a> rnorm;
<a name="line12"> 12: </a> <a href="tutorials/../../docs/manualpages/Mat/Mat.html#Mat">Mat</a> A;
<a name="line13"> 13: </a> <a href="tutorials/../../docs/manualpages/Vec/Vec.html#Vec">Vec</a> b,x,r;
<a name="line14"> 14: </a> <a href="tutorials/../../docs/manualpages/KSP/KSP.html#KSP">KSP</a> ksp;
<a name="line15"> 15: </a> <a href="tutorials/../../docs/manualpages/PC/PC.html#PC">PC</a> pc;
<a name="line20"> 20: </a> <a href="tutorials/../../docs/manualpages/Sys/PetscOptionsGetInt.html#PetscOptionsGetInt">PetscOptionsGetInt</a>(NULL,NULL,<font color="#666666">"-N"</font>,&N,NULL);
<a name="line21"> 21: </a> <a href="tutorials/../../docs/manualpages/Sys/PetscOptionsGetBool.html#PetscOptionsGetBool">PetscOptionsGetBool</a>(NULL,NULL,<font color="#666666">"-draw"</font>,&draw,NULL);
<a name="line23"> 23: </a> <a href="tutorials/../../docs/manualpages/Mat/MatCreate.html#MatCreate">MatCreate</a>(<a href="tutorials/../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</a>,&A);
<a name="line24"> 24: </a> <a href="tutorials/../../docs/manualpages/Mat/MatSetSizes.html#MatSetSizes">MatSetSizes</a>(A,<a href="tutorials/../../docs/manualpages/Sys/PETSC_DECIDE.html#PETSC_DECIDE">PETSC_DECIDE</a>,<a href="tutorials/../../docs/manualpages/Sys/PETSC_DECIDE.html#PETSC_DECIDE">PETSC_DECIDE</a>,N,N);
<a name="line25"> 25: </a> <a href="tutorials/../../docs/manualpages/Mat/MatSetType.html#MatSetType">MatSetType</a>(A,MATPYTHON);
<a name="line26"> 26: </a> <a href="tutorials/../../docs/manualpages/Mat/MatPythonSetType.html#MatPythonSetType">MatPythonSetType</a>(A,<font color="#666666">"example1.py:Laplace1D"</font>);
<a name="line27"> 27: </a> <a href="tutorials/../../docs/manualpages/Mat/MatSetUp.html#MatSetUp">MatSetUp</a>(A);
<a name="line29"> 29: </a> <a href="tutorials/../../docs/manualpages/Mat/MatCreateVecs.html#MatCreateVecs">MatCreateVecs</a>(A,&x,&b);
<a name="line30"> 30: </a> <a href="tutorials/../../docs/manualpages/Vec/VecSet.html#VecSet">VecSet</a>(b,1);
<a name="line32"> 32: </a> <a href="tutorials/../../docs/manualpages/KSP/KSPCreate.html#KSPCreate">KSPCreate</a>(<a href="tutorials/../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</a>,&ksp);
<a name="line33"> 33: </a> <a href="tutorials/../../docs/manualpages/KSP/KSPSetType.html#KSPSetType">KSPSetType</a>(ksp,KSPPYTHON);
<a name="line34"> 34: </a> <a href="tutorials/../../docs/manualpages/KSP/KSPPythonSetType.html#KSPPythonSetType">KSPPythonSetType</a>(ksp,<font color="#666666">"example1.py:ConjGrad"</font>);
<a name="line36"> 36: </a> <a href="tutorials/../../docs/manualpages/KSP/KSPGetPC.html#KSPGetPC">KSPGetPC</a>(ksp,&pc);
<a name="line37"> 37: </a> <a href="tutorials/../../docs/manualpages/PC/PCSetType.html#PCSetType">PCSetType</a>(pc,PCPYTHON);
<a name="line38"> 38: </a> <a href="tutorials/../../docs/manualpages/PC/PCPythonSetType.html#PCPythonSetType">PCPythonSetType</a>(pc,<font color="#666666">"example1.py:Jacobi"</font>);
<a name="line40"> 40: </a> <a href="tutorials/../../docs/manualpages/KSP/KSPSetOperators.html#KSPSetOperators">KSPSetOperators</a>(ksp,A,A);
<a name="line41"> 41: </a> <a href="tutorials/../../docs/manualpages/KSP/KSPSetFromOptions.html#KSPSetFromOptions">KSPSetFromOptions</a>(ksp);
<a name="line42"> 42: </a> <a href="tutorials/../../docs/manualpages/KSP/KSPSolve.html#KSPSolve">KSPSolve</a>(ksp,b,x);
<a name="line44"> 44: </a> <a href="tutorials/../../docs/manualpages/Vec/VecDuplicate.html#VecDuplicate">VecDuplicate</a>(b,&r);
<a name="line45"> 45: </a> <a href="tutorials/../../docs/manualpages/Mat/MatMult.html#MatMult">MatMult</a>(A,x,r);
<a name="line46"> 46: </a> <a href="tutorials/../../docs/manualpages/Vec/VecAYPX.html#VecAYPX">VecAYPX</a>(r,-1,b);
<a name="line47"> 47: </a> <a href="tutorials/../../docs/manualpages/Vec/VecNorm.html#VecNorm">VecNorm</a>(r,<a href="tutorials/../../docs/manualpages/Vec/NORM_2.html#NORM_2">NORM_2</a>,&rnorm);
<a name="line48"> 48: </a> <a href="tutorials/../../docs/manualpages/Sys/PetscPrintf.html#PetscPrintf">PetscPrintf</a>(<a href="tutorials/../../docs/manualpages/Sys/PETSC_COMM_WORLD.html#PETSC_COMM_WORLD">PETSC_COMM_WORLD</a>,<font color="#666666">"error norm = %g\n"</font>,rnorm);
<a name="line50"> 50: </a> <font color="#4169E1">if</font> (draw) {
<a name="line51"> 51: </a> <a href="tutorials/../../docs/manualpages/Vec/VecView.html#VecView">VecView</a>(x,<a href="tutorials/../../docs/manualpages/Viewer/PETSC_VIEWER_DRAW_WORLD.html#PETSC_VIEWER_DRAW_WORLD">PETSC_VIEWER_DRAW_WORLD</a>);
<a name="line52"> 52: </a> <a href="tutorials/../../docs/manualpages/Sys/PetscSleep.html#PetscSleep">PetscSleep</a>(2);
<a name="line53"> 53: </a> }
<a name="line55"> 55: </a> <a href="tutorials/../../docs/manualpages/Vec/VecDestroy.html#VecDestroy">VecDestroy</a>(&x);
<a name="line56"> 56: </a> <a href="tutorials/../../docs/manualpages/Vec/VecDestroy.html#VecDestroy">VecDestroy</a>(&b);
<a name="line57"> 57: </a> <a href="tutorials/../../docs/manualpages/Vec/VecDestroy.html#VecDestroy">VecDestroy</a>(&r);
<a name="line58"> 58: </a> <a href="tutorials/../../docs/manualpages/Mat/MatDestroy.html#MatDestroy">MatDestroy</a>(&A);
<a name="line59"> 59: </a> <a href="tutorials/../../docs/manualpages/KSP/KSPDestroy.html#KSPDestroy">KSPDestroy</a>(&ksp);
<a name="line61"> 61: </a> <font color="#4169E1">return</font>(0);
<a name="line62"> 62: </a>}
<a name="line64"> 64: </a><font color="#B22222">/* ------------------------------------------------------- */</font>
<a name="line66"> 66: </a>static char help[] = <font color="#666666">"Python-implemented <a href="tutorials/../../docs/manualpages/Mat/Mat.html#Mat">Mat</a>/<a href="tutorials/../../docs/manualpages/KSP/KSP.html#KSP">KSP</a>/<a href="tutorials/../../docs/manualpages/PC/PC.html#PC">PC</a>.\n\n"</font>;
<a name="line68"> 68: </a><font color="#B22222">/*</font>
<a name="line69"> 69: </a><font color="#B22222">#define PYTHON_EXE "python2.5"</font>
<a name="line70"> 70: </a><font color="#B22222">#define PYTHON_LIB "/usr/lib/libpython2.5"</font>
<a name="line71"> 71: </a><font color="#B22222">*/</font>
<a name="line73"> 73: </a><font color="#A020F0">#if !defined(PYTHON_EXE)</font>
<a name="line74"> 74: </a><strong><font color="#228B22">#define PYTHON_EXE 0</font></strong>
<a name="line75"> 75: </a><font color="#A020F0">#endif</font>
<a name="line76"> 76: </a><font color="#A020F0">#if !defined(PYTHON_LIB)</font>
<a name="line77"> 77: </a><strong><font color="#228B22">#define PYTHON_LIB 0</font></strong>
<a name="line78"> 78: </a><font color="#A020F0">#endif</font>
<a name="line82"> 82: </a><strong><font color="#4169E1"><a name="main"></a>int main(int argc, char *argv[])</font></strong>
<a name="line83"> 83: </a>{
<a name="line86"> 86: </a> <a href="tutorials/../../docs/manualpages/Sys/PetscInitialize.html#PetscInitialize">PetscInitialize</a>(&argc,&argv,0,help);
<a name="line88"> 88: </a> PetscPythonInitialize(PYTHON_EXE,PYTHON_LIB);
<a name="line90"> 90: </a> RunTest();PetscPythonPrintError();
<a name="line92"> 92: </a> <a href="tutorials/../../docs/manualpages/Sys/PetscFinalize.html#PetscFinalize">PetscFinalize</a>();
<a name="line94"> 94: </a> <font color="#4169E1">return</font> 0;
<a name="line95"> 95: </a>}
<a name="line97"> 97: </a><font color="#B22222">/* ------------------------------------------------------- */</font>
</pre>
</body>
</html>
|
