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
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
|
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head> <link rel="canonical" href="http://www.mcs.anl.gov/petsc/petsc-current/docs/changes/2022.html" />
<meta http-equiv="content-type" content="text/html;charset=utf-8">
<title>Documentation: Changes: 2.0.22</title>
</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 docs/changes/2022.html "><small>Report Typos and Errors</small></a></div>
<h1>Documentation: Changes: 2.0.22</h1>
<style type="text/css">
h2 {
color: red;
}
h4 {
font-weight: bold;
text-decoration: underline;
}
</style>
<div id="main">
<div align="center">
<h2>NEW FEATURES and <a name="CHANGES">CHANGES</a> in PETSc 2.0.22</h2>
</div>
<h4>General:</h4>
<ul>
<li>
The Fortran90-specific routines such as <code>VecGetArrayF90()</code>
now work with the IBM Fortran compiler, the SGI Fortran 90 compiler,
and the Cray T3E Fortran compiler.
</li>
<li>
The third argument for the Fortran versions of the routines
<code>VecGetArray()/VecRestoreArray(),</code>
<code>ISGetIndices()/ISRestoreIndices(),
MatGetArray()/MatRestoreArray()</code> and
<code>DAGetGlobalIndices()</code> has been changed from
<code>integer</code> to <code>PetscOffset</code>. See
src/snes/examples/tutorials/ex5f.F for example usage.
</li>
<li>
Changed the location of the manual pages from docs/www/www.html and
docs/www to <code>docs/manualpages/manualpages.html</code> and
<code>docs/manualpages</code>.
</li>
<li>
The location of the mpiuni include files (used in the base.site file
when not compiling with an MPI) has moved from
<code>src/mpiuni</code> to <code>src/sys/src/mpiuni</code>
</li>
<li>
Added the option <code>-log_summary_exclude <vec,mat,sles,snes></code> to
limit the information printed in the summary table
</li>
<li>Added the option <code>-log_info_exclude <vec,mat,sles,snes></code></li>
<li>
Changed the calling sequences of <code>PetscGetTime(),
PetscGetCPUTime(), PetscGetFlops()</code> to return error codes
like all other PETSc functions.
</li>
<li>
Changed the internal structure of the PETSc objects. This will not
effect most users but makes PETSc more extensible for the
future.
</li>
<li>
Removed the PETSc object child support and replaced it with
<ul>
<li><code>PetscObjectCompose(PetscObject,char *name,PetscObject);</code></li>
<li><code>PetscObjectQuery(PetscObject,char *name,PetscObject *);</code></li>
<li><code>PetscObjectComposeFunction(PetscObject,char *name, void *);</code></li>
<li><code>PetscObjectQueryFunction(PetscObject,char *name,void **);</code></li>
</ul>
These allow the user to attach objects and functions to any PETSc object.
</li>
<li>
Added a feeble example that demonstrates how to manage partitioning
a grid in parallel in <code>src/mat/examples/tutorials/ex2.c</code>
</li>
<li>
The build command 'make all' in ${PETSC_DIR} will build the PETSc
libraries, including the fortran interface.
</li>
</ul>
<h4>AO (Application Orderings):</h4>
<h4>TS (Timestepping Solvers):</h4>
<ul>
<li>The calling sequence of <code>TSSetType()</code> has changed; see the manual page.</li>
<li><code>TSType</code> is now a string, rather than enum type</li>
</ul>
<h4>SNES (Nonlinear Solvers):</h4>
<ul>
<li>The calling sequence of <code>SNESSetType()</code> has changed; see the manual page.</li>
<li><code>SNESType</code> is now a string, rather than enum type</li>
</ul>
<h4>SLES (Linear Solvers):</h4>
<ul>
<li>See PC and KSP</li>
</ul>
<h4>KSP (Krylov Subspace Methods):</h4>
<ul>
<li>The calling sequence of <code>KSPSetType()</code> has changed; see the manual page.</li>
<li><code>KSPType</code> is now a string, rather than enum type.</li>
</ul>
<h4>PC (Preconditioners):</h4>
<ul>
<li>
Added two new PC classes, <code>PCSLES</code> and
<code>PCCOMPOSITE</code> that allow use of any
<code>SLES</code> object as a preconditioner and enable combining
several preconditioners.
</li>
<li>The calling sequence of <code>PCSetType()</code> has changed; see the manual page.</li>
<li><code>PCType</code> is now a string, rather than enum type.</li>
</ul>
<h4>MAT (Matrices):</h4>
<ul>
<li>
Added support for partitioning using the ParMETIS parallel
partitioning package, see the manual page for
<code>PartitioningCreate()</code>
</li>
<li>Added <code>MatGetColumnVector()</code></li>
<li>Added argument column size to <code>MatGetSubMatrix()</code></li>
<li>
Changed the memory allocation algorithm in
<code>MatLUFactorSymbolic()</code> and
<code>MatILUFactorSymbolic()</code> based on a suggestion by David
Hysom of ODU.
</li>
<li>
The third argument for the Fortran version of the routines
<code>MatGetArray()/MatRestoreArray()</code> has been changed from
<code>integer</code> to <code>PetscOffset</code>.
</li>
</ul>
<h4>DA (Distributed Arrays):</h4>
<ul>
<li>
Changed <code>DAGetDistributedVector()</code> and
<code>DAGetLocalVector()</code> to
<code>DACreateGlobalVector()</code> and
<code>DACreateLocalVector()</code>. You should destroy these
vectors with a usual call to <code>VecDestroy()</code>.
</li>
<li>
The third argument for the Fortran version of the routine
<code>DAGetGlobalIndices()</code> has been changed from
<code>integer</code> to <code>PetscOffset</code>.
</li>
</ul>
<h4>VEC (Vectors):</h4>
<ul>
<li>
Added <code>VecSetBlockSize(), VecSetValuesBlocked(),
VecSetValuesBlockedLocal(),
VecSetLocalToGlobalMappingBlocked()</code>
</li>
<li>
Added<code> VecCreatedShared(),</code> which creates
a parallel vector in shared memory on the SGI machines
</li>
<li>
Changed
<code>VEC_IGNORE_OFF_PROCESSOR_ENTRIES</code> to
<code>VEC_IGNORE_OFF_PROC_ENTRIES</code>
</li>
<li>
The third argument for the Fortran versions of the routines
<code>VecGetArray()/VecRestoreArray(),</code> has been changed from
<code>integer</code> to <code>PetscOffset</code>. See
src/snes/examples/tutorials/ex5f.F for example usage.
</li>
<li>
<code>VecGetArray(), VecRestoreArray()</code> now work from Fortran on
all machines, finally after 3 long years!
</li>
<li>
Changed the calling sequence of
<code>VecCreate()</code>, now takes both local and global size,
see the manual page.
</li>
</ul>
<h4>IS (Index Sets):</h4>
<ul>
<li>
The third argument for the Fortran versions of the routines
<code>ISGetIndices()/ISRestoreIndices()</code> has been changed from
<code>integer</code> to <code>PetscOffset</code>.
</li>
</ul>
<h4>Draw (Graphics):</h4>
<h4>Viewers:</h4>
<h4>System:</h4>
<h4>Error Handling:</h4>
<h4>Event Logging:</h4>
<h4>Fortran Interface:</h4>
<ul>
<li>
<code>PETSC_NULL</code> has been replaced with
<code>PETSC_NULL_INTEGER</code>,
<code>PETSC_NULL_SCALAR</code>,
<code>PETSC_NULL_DOUBLE</code> or
<code>PETSC_NULL_CHARACTER</code> depending on the context for
usage. If the function expects an integer, a scalar (double or
complex depending if you compile with BOPT=*_complex), a double
precision array, or a string.
</li>
<li>
The Fortran include files are now located in the
directory "include/finclude/*.", rather than
"include/FINCLUDE/*.h". The include files that used to be
in "include/finclude/*.h" (not recommended) are now in
"include/foldinclude/*.h"
</li>
<li>
If you have stored PETSc objects in Fortran integer
arrays, such as <code>user(3) = vec</code> you have
to declare those arrays to be of type
<code>PetscFortranAddr</code> instead of as integer. Note that
<code>PetscFortranAddr</code> is simply integer on 32 bit
machines and integer*8 on 64 bit machines.
</li>
<li>
The Fortran90-specific routines such as
<code>VecGetArrayF90()</code> now work with the IBM Fortran
compiler, the SGI Fortran 90 compiler, and the Cray T3E Fortran
compiler.
</li>
<li>
PETSc objects can now be passed between C and Fortran directly. You do
not have to use the commands
<code>PetscCObjectToFortranObject()</code> and
<code>PetscFortranObjectToCObject()</code> to translate the objects
between languages.
</li>
<li>
Changed Fortran object argument from type int to type PetscFortranAddr
in the routines <code>PetscCObjectToFortranObject()</code> and
<code>PetscFortranObjectToCObject().</code>
</li>
<li>
Added the routine <code>PetscInitializeFortran()</code> to assist
in mixed-language use of PETSc. See the manual page for details.
</li>
</ul>
</div>
</body>
</html>
|
