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
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
|
<!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/32.html" />
<meta http-equiv="content-type" content="text/html;charset=utf-8">
<title>Documentation: Changes: 3.2</title>
</head>
<body bgcolor="#ffffff">
<div id="version" align=right><b>petsc-3.4.2 2013-07-02</b></div>
<h1>Documentation: Changes: 3.2</h1>
<style type="text/css">
h2 {
color: red;
}
h4 {
font-weight: bold;
text-decoration: underline;
}
</style>
<div id="main">
<div align="center">
<h2><a name="CHANGES">CHANGES in PETSc 3.2</a></h2>
</div>
<h4>General:</h4>
<ul>
<li>--download-c-blas-lapack is replaced with --download-f2cblaslapack</li>
<li>--with-precision=__float128 --download-f2cblaslapack is now support for newish GNU gcc compiler</li>
<li>Calling sequence of PetscBinarRead.m has been changed</li>
<li>
All XXXDestroy() functions now take a pointer to the object itself so that it can be nullified after destruction.
</li>
<li>
PetscError() and SETERRQX() now take a MPI_Comm as the first argument
to indicate where the error is known. If you don't know what
communicator use then pass in PETSC_COMM_SELF
</li>
<li>
Added PetscObjectAddOptionsHandler() that allows adding new routines
to any object that are called to process options when
XXXSetFromOptions() is called.
</li>
<li>
Changed PetscTruth to PetscBool, PETSC_TRUTH to PETSC_BOOL,
PetscOptionsTruth to PetscOptionsBool, etc.
</li>
<li>Introduced PetscCopyMode</li>
<li>
Changed -log_history to -history and PetscLogOpenHistory()
PetscLogCloseHistory() to PetscOpenHistory(), PetscCloseHistory()
because "log" is reserved for related to the performance logging
</li>
<li>Changed PetscOptionsPrint() to PetscOptionsView()</li>
<li>Changed PetscLogPrintSummary() to PetscLogView()</li>
<li>
Using gcc 4.6 you can now ./configure --with-precision=__float128
--download-qblaslapack to get computations in quad precision.
gfortran 4.6 is also supported. Warning some of the PETSc printing
of numbers in some places will print garbage. I haven't figured out
a way to handle the %G format. Can be run in parallel :-)
</li>
<li>
PetscFListFind() now takes a searchlibraries flag, if this is set and
the function is NOT found in the function list then the libraries and
executable are searched (when shared or dynamic libraries are used).
Previously it always search the libraries and executables if not
found in the list.
</li>
<li>PetscOpenMPxxx() utilities are now PetscHMPIxxx() for hierarchical MPI, the previous name was bad</li>
</ul>
<h4>Logging:</h4>
<ul>
<li>PetscCookie changed to PetscClassId in all forms.</li>
</ul>
<h4>config/configure.py:</h4>
<ul>
<li>
--with-64-bit-pointers option removed. It doesn't work properly
anway. The user should specify the correct 32bit or 64bit compilers
to configure. For eg: <code>configure --with-cc='gcc -m64'
--with-fc='gfortran -m64'</code>
</li>
<li>Python requirement is changed from 2.2 to 2.3</li>
<li>
Changed --with-shared to --with-shared-libraries, --known-mpi-shared
to --known-mpi-shared-libraries and --with-dynamic to --with-dynamic-loading
</li>
<li>
'ifneeded' support removed. '--download-package=ifneeded' should now
be '--download-package'
</li>
<li>The make macro PETSC_INCLUDE has been replaced by PETSC_CC_INCLUDES and PETSC_FC_INCLUDES for C and Fortran, respectively.</li>
</ul>
<h4>IS:</h4>
<ul>
<li>
ISBlock(),ISStride() are removed. Use
PetscTypeCompare((PetscObject)is,ISBLOCK,&flag),
PetscTypeCompare((PetscObject)is,ISSTRIDE,&flag) instead.
</li>
<li>
Added ISCreate(), ISSetType(), ISRegister() etc to match style of
other PETSc objects. Also added ISGeneralSetIndices(),
ISBlockSetIndices(), ISStrideSetStride()
</li>
<li>
ISCreateGeneral() and ISCreateBlock() now takes PetscCopyMode as an
argument and hence ISCreateGeneralNC() and ISCreateGeneralWithArray()
are not needed and removed. Use ISCreateGeneral() with the argument
PETSC_OWN_POINTER and PETSC_USE_POINTER instead.
</li>
<li>
ISLocalToGlobalMappingCreate() now takes PetscCopyMode as an
additional argument and ISLocalToGlobalMappingCreateNC() is no longer
needed and removed
</li>
<li>
ISStrideToGeneral() is replaced with ISToGeneral() that works for all
basic IS types.
</li>
<li>
ISCreateBlock() now takes indexing relative to block, no longer
relative to element. This is to match the paradigm of
Vec/MatSetValuesBlocked()
</li>
<li>
ISBlockGetBlockSize() is now ISGetBlockSize(). Block sizes can be set
for conforming ISGENERAL and ISSTRIDE using ISSetBlockSize().
</li>
<li>ISAllGatherIndices() is removed, use ISCreateGeneral() then ISAllGather().</li>
</ul>
<h4>PF:</h4>
<ul>
<li>The source array is marked constant in PFApply and the callback.</li>
</ul>
<h4>Vec:</h4>
<ul>
<li>changed VecSqrt() to VecSqrtAbs()</li>
<li>
VecLoad() and VecLoadIntoVector() have been merged into a single NEW
VecLoad() that takes a partially constructed vector as input and
loads according to any type or sizes that have previously been set
into the Vec.
<li>VecDestroyVecs(PetscInt n,Vec **) instead of VecDestroyVecs(Vec**,PetscInt n), also VecDestroyVecsF90(PetscInt n,{Vec, pointer :: x},PetscErrorCode ierr)</li>
<li>
VecLoad() does NOT access the options database to check for the
VecType, call VecSetFromOptions() first if you wish to check it.
</li>
<li>
VecDestroyVecs() now takes the pointer to the first array and zeros
it on return
</li>
<li>
VecDestroy() now zeros the pointer to the destroyed vector so it cannot be reused
</li>
<li>VecDestroyVecs(PetscInt n,Vec **) instead of VecDestroyVecs(Vec**,PetscInt n)</li>
</ul>
<h4>VecScatter:</h4>
<ul>
<li>
The ghost indices for VecCreateGhostBlock() and
VecCreateGhostBlockWithArray() are now by block instead of by
entry.
</li>
</ul>
<h4>Mat:</h4>
<ul>
<li>MAT_SOLVER_XXX is now MATSOLVERXXX</li>
<li>MAT_PARTIONING_XXX is now MATPARTITIONINGXXX</li>
<li>MATCOLORING_XXX is now MATCOLORINGXXX</li>
<li>MATORDERING_XXX is now MATORDERINGXXX</li>
<li>
It is an error to call MatSetOption() before the implementation has
been created (after MatSetSizes() and MatSetType()). Formerly,
options were silently ignored in this circumstance.
</li>
<li>Added MatSetMumpsIcntl()</li>
<li>
MatLoad() now takes a created Mat as input (and no longer a Mat
type), if you provide the size and type in the Mat then that is used
other defaults are used.
</li>
<li>MATCRL -> MATAIJCRL, MATSEQCRL -> MATSEQAIJCRL, MATMPICRL -> MATMPIAIJCRL</li>
<li>MatCreateSeqCRL() -> MatCreateSeqAIJCRL() MatCreateMPICRL() -> MatCreateMPIAIJCRL()</li>
<li>MATCSRPERM -> MATAIJPERM, MATSEQCSRPERM -> MATSEQAIJPERM, MATMPICSRPERM -> MATMPIAIJPERM</li>
<li>MatCreateSeqCSRPERM() -> MatCreateSeqAIJPERM() MatCreateMPICSRPERM() -> MatCreateMPIAIJPERM()</li>
<li>Added MatZeroRowsColumns() and MatZeroRowsColumnsIS().</li>
<li>
MatZeroRows() and MatZeroRowsIS() now take an x and b vector as
optional arguments; if these are provided then the b[idx[i]] is set
to diag*x[idx[i]]] for all rows listed in idx[].
</li>
<li>
MatSetLocalToGlobalMapping() and MatSetLocalToGlobalMappingBlock()
now take separate row and column maps. MatPreallocateSetLocal() has
an extra argument for the column mapping.
</li>
<li>
MatLoad() does NOT access the options database to check for the
MatType, call MatSetFromOptions() first if you wish to check it.
</li>
<li>Removed MatMFFDSetFromOptions() since one can call MatSetFromOptions()</li>
<li>
MatGetLocalMat() and MatGetLocalMatCondensed() are now
MatMPIAIJGetLocalMat() and MatMPIAIJGetLocalMatCondensed()
</li>
<li>MatOption MAT_USE_COMPRESSEDROW changed to MAT_CHECK_COMPRESSED_ROW<br></li>
</ul>
<h4>PC:</h4>
<ul>
<li>
PCFieldSplitSetIS() and PCFieldSplitSetFields() now take a split-name
parameter which is used to define the options database keys.
</li>
<li>CHOLMOD can be used for Cholesky factorization using AIJ or SBAIJ(1) matrix formats.</li>
<li>
PCREDUNDANT now uses the inner KSP with the inner PC. Default is
still preonly for KSP and LU for PC so default behavior is the same
but now you can use for example -redundant_ksp_gmres.
PCRedundantGetPC() is now changed to PCRedundantGetKSP()
</li>
<li>The header "petscmg.h" is renamed to "petscpcmg.h", "petscasa.h" is renamed to "petscpcasa.h".</li>
<li>The preconditioner PCBFBT has been removed. The same functionality is available in PCFIELDSPLIT.</li>
<li>PCOPENMP is now PCHMPI and its prefix is now -hmpi_ for heirarchical MPI</li>
<li>PCMGSetGalerkin() has a second argument of type PetscBool.</li>
</ul>
<h4>KSP:</h4>
<ul>
<li>KSPSetPreconditionerSide() changed to KSPSetPCSide() to match name of second argument</li>
<li>-ksp_right_pc -ksp_left_pc changed to -ksp_pc_side left,right,symmetric to match KSPSetPCSide()</li>
<li>Added KSPGMRESGetRestart() and KSPGMRESGetCGSRefinementType()</li>
<li>Added KSPGMRESGetOrthogonalization()</li>
<li>KSPAddOptionsChecker() is replaced with PetscObjectAddOptionsHandler()</li>
<li>
Added KSPSPECEST which estimates the spectrum on the first solve and
uses it to configure a reduction-free method for subsequent solves,
intended for use when reductions are expensive such as levels of
multigrid.
</li>
</ul>
<h4>SNES:</h4>
<ul>
<li>Added SNESLineSearchSetMonitor() and -snes_ls_monitor.</li>
<li>Add SNESKSPONLY, for solving linear problems with SNES.</li>
<li>Added the minlambda argument to SNESLineSearchGetParams() and SNESLineSearchSetParams().</li>
</ul>
<h4>TS:</h4>
<ul>
<li>Rename TSCRANK_NICHOLSON to TSCN and TSRUNGE_KUTTA to TSRK for consistency.</li>
<li>
TSSetIFunction() and TSSetRHSFunction() have an additional argument
for the vector to compute the function value in. If PETSC_NULL is
used, then one is created by the TS.
</li>
<li>
TSSetMatrices() has been removed, use TSSetIFunction() and
TSSetRHSFunction(), perhaps providing Jacobian matrices and
TSComputeRHSFunctionLinear() and/or TSComputeRHSJacobianConstant().
See src/ts/examples/tutorials/ex{2,4,5,6}.c for an example.
</li>
<li>Added TSARKIMEX: additive Runge-Kutta implicit-explicit methods for multi-rate systems.</li>
<li>
Changed TSStep() interface to only perform one step, added TSSolve()
to perform multiple steps with callbacks and monitors if desired.
</li>
<li>Added TSSetExactFinalTime(), replaces TSSundialsSetExactFinalTime().</li>
</ul>
<h4>DM/DA:</h4>
<ul>
<li>
Change array argument of DAGetArray and related functions from void**
to void* to avoid the need for an explicit cast (the argument still
has the meaning of a pointer to d-dimensionally indexed array of
user-defined node type, old code will still compile correctly but the
cast is no longer necessary).
</li>
<li>
DAGetCoordinates(), DAGetGhostedCoordinates(), and
DAGetCoordinateDA() now return borrowed references. The returned
object should not be destroyed by the user.
</li>
<li>Added DAVecGetArrayF90()</li>
<li>Changed DASetVertexDivision to DASetOwnershipRanges</li>
<li>The SDA object and all its method have been removed</li>
<li>The DA, ADDA, Slice and DMComposite objects are now all represented by a DM object.</li>
<li>Routines that began with ADDA and Slice now begin with DMADDA and DMSlice</li>
<li>
DA/ADDA/Slice/DMCompositeGetMatrix() -> DMGetMatrix(), similar for
all other routines that have a DM version
</li>
<li>
Removed DASetType(), DARegister() etc since they will all be handled
via DMSetType() etc
</li>
<li>
Added DMSetUp() that is called, for example, after all the parameters
are passed to the DM object to actually construct the data
structures; replace all calls to DASetType() with calls to
DMSetUp()
</li>
<li>DAXXX() routines are now either DMXXX() or DMDAXXX() see petscdmda.h for details</li>
<li>
The operation DALocalToGlobal() which performed no communication can
now be performed with DMLocalToGlobalBegin/End() with INSERT_VALUES.
</li>
<li>petscda.h, petscdadef.h and petscda.h90 are now petscdmda.h, petscdmdadef.h, and petscdmda.h90</li>
<li>
DAGetISLocalToGlobalMapping() and DAGetISLocalToGlobalMappingBlck()
are now DMGetLocalToGlobalMapping() and DMGetLocalToGlobalMappingBlock()
and are collective since they may be constructed lazily.
</li>
<li>DMGetBlockSize() is added.</li>
<li>
Changed DMDAPeriodicType to DMDABoundaryType, and changed this enum
to only specify one boundary (one for each dimension). This enum is
now one of DMDA_BOUNDARY_NONE, DMDA_BOUNDARY_GHOSTED, or
DMDA_BOUNDARY_PERIODIC. See DMDABoundaryType man page for more.
</li>
<li>
DMDASetPeriodicity() is now DMDASetBoundaryType(), and API is changed
to take three DMDABoundaryType arguments.
</li>
<li>
Changed API for DMDACreate*D(), DMDAGetInfo() to take the correct
number of DMDABoundaryType enums for the dimension.
</li>
<li>DMDASetGhostedCoordinates() is added to set coordinates of non-periodic ghost cells.</li>
</ul>
<h4>DMMG:</h4>
<h4>PetscViewer:</h4>
<ul>
<li>PETSC_VIEWER_XXX is now PETSCVIEWERXXX</li>
<li>Added PetscViewerASCIIOpenWithFILE() and PetscViewerASCIISetFILE()</li>
<li>
Added PetscViewerASCIISynchronizedAllow() which must be
called before using PetscViewerASCIISynchronizedPrintf()
</li>
</ul>
<h4>SYS:</h4>
<ul>
<li>
PetscOptionsAtoi(),PetscOptionsAtod(),PetscOptionsAtol() ->
PetscOptionsStringToInt(), PetscOptionsStringToReal(),
PetscOptionsStringToBool()
</li>
<li>
PetscFListFind() takes an additional argument to look up the symbol
in dynamically loaded libraries if not already in the list.
</li>
<li>Added PetscBagSetOptionsPrefix(), which allows a prefix for the option names of bag items.</li>
<li>
The preloading macros as well as developer-level logging and
profiling functions have been namespaced, e.g. PreLoadBegin() is now
PetscPreLoadBegin().
</li>
</ul>
<h4>AO:</h4>
<h4>Sieve:</h4>
<h4>Fortran:</h4>
<h4><a href="http://www.mcs.anl.gov/petsc/miscellaneous/external.html">ExternalPackages</a>:</h4>
</div>
</body>
</html>
|
