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 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352
|
Other PETSc Features
--------------------
PETSc on a process subset
~~~~~~~~~~~~~~~~~~~~~~~~~
Users who wish to employ PETSc on only a subset of MPI processes
within a larger parallel job, or who wish to use a “manager” process to
coordinate the work of “worker” PETSc processes, should specify an
alternative communicator for ``PETSC_COMM_WORLD`` by directly setting
its value, for example to use an existing MPI communicator ``comm``,
.. code-block::
PETSC_COMM_WORLD = comm; /* To use a previously-defined MPI_Comm */
*before* calling ``PetscInitialize()``, but, obviously, after calling
``MPI_Init()``.
.. _sec_options:
Runtime Options
~~~~~~~~~~~~~~~
Allowing the user to modify parameters and options easily at runtime is
very desirable for many applications. PETSc provides a simple mechanism
to enable such customization. To print a list of available options for a
given program, simply specify the option ``-help`` at
runtime, e.g.,
.. code-block:: console
$ mpiexec -n 1 ./ex1 -help
Note that all runtime options correspond to particular PETSc routines
that can be explicitly called from within a program to set compile-time
defaults. For many applications it is natural to use a combination of
compile-time and runtime choices. For example, when solving a linear
system, one could explicitly specify use of the Krylov subspace
solver BiCGStab by calling
.. code-block::
KSPSetType(ksp, KSPBCGS);
One could then override this choice at runtime with the option
::
-ksp_type tfqmr
to select the Transpose-Free QMR algorithm. (See
:any:`ch_ksp` for details.)
The remainder of this section discusses details of runtime options.
.. _the-options-database-1:
The Options Database
^^^^^^^^^^^^^^^^^^^^
Each PETSc process maintains a database of option names and values
(stored as text strings). This database is generated with the command
``PetscInitialize()``, which is listed below in its C/C++ and Fortran
variants, respectively:
.. code-block::
PetscInitialize(int *argc, char ***args, const char *file, const char *help); // C
.. code-block:: fortran
call PetscInitialize(integer ierr) ! Fortran
The arguments ``argc`` and ``args`` (in the C/C++ version only) are the
addresses of the usual command line arguments, while the ``file`` is a name
of an optional file that can contain additional options. By default this file is
called ``.petscrc`` in the user’s home directory. The user can also
specify options via the environmental variable ``PETSC_OPTIONS``. The
options are processed in the following order:
#. file
#. environmental variable
#. command line
Thus, the command line options supersede the environmental variable
options, which in turn supersede the options file.
The file format for specifying options is
.. code-block:: none
-optionname possible_value
-anotheroptionname possible_value
...
All of the option names must begin with a dash (-) and have no
intervening spaces. Note that the option values cannot have intervening
spaces either, and tab characters cannot be used between the option
names and values. For
uniformity throughout PETSc, we employ the format
``-[prefix_]package_option`` (for instance, ``-ksp_type``,
``-mat_view ::info``, or ``-mg_levels_ksp_type``).
Users can specify an alias for any option name (to avoid typing the
sometimes lengthy default name) by adding an alias to the ``.petscrc``
file in the format
.. code-block:: none
alias -newname -oldname
For example,
.. code-block:: none
alias -kspt -ksp_type
alias -sd -start_in_debugger
Comments can be placed in the ``.petscrc`` file by using ``#`` in the
first column of a line.
Options Prefixes
^^^^^^^^^^^^^^^^
Options prefixes allow specific objects to be controlled from the
options database. For instance, ``PCMG`` gives prefixes to its nested
``KSP`` objects; one may control the coarse grid solver by adding the
``mg_coarse`` prefix, for example ``-mg_coarse_ksp_type preonly``. One
may also use ``KSPSetOptionsPrefix()``,\ ``DMSetOptionsPrefix()`` ,
``SNESSetOptionsPrefix()``, ``TSSetOptionsPrefix()``, and similar
functions to assign custom prefixes, useful for applications with
multiple or nested solvers.
Adding options from a file
^^^^^^^^^^^^^^^^^^^^^^^^^^
PETSc can load additional options from a file using ``PetscOptionsInsertFile()``,
which can also be used from the command line, e.g. ``-options_file my_options.opts``.
One can also use YAML files with ``PetscOptionsInsertFileYAML()``.
For example, the following file:
.. literalinclude:: /../src/sys/tests/ex47-options.yaml
:language: yaml
corresponds to the following PETSc options:
.. literalinclude:: /../src/sys/tests/output/ex47_3_options.out
:language: none
:start-after: #
:end-before: #End
With ``-options_file``, PETSc will parse the file as YAML if it ends in a standard
YAML or JSON [#json]_ extension or if one uses a ``:yaml`` postfix,
e.g. ``-options_file my_options.yaml`` or ``-options_file my_options.txt:yaml``
PETSc will also check the first line of the options file itself and
parse the file as YAML if it matches certain criteria, for example.
.. literalinclude:: /../src/sys/tests/ex47-yaml_tag
:language: yaml
and
.. literalinclude:: /../src/sys/tests/ex47-yaml_doc
:language: yaml
both correspond to options
.. literalinclude:: /../src/sys/tests/output/ex47_2_auto.out
:language: none
:start-after: #
:end-before: #End
User-Defined PetscOptions
^^^^^^^^^^^^^^^^^^^^^^^^^
Any subroutine in a PETSc program can add entries to the database with
the command
.. code-block::
PetscOptionsSetValue(PetscOptions options, char *name, char *value);
though this is rarely done. To locate options in the database, one
should use the commands
.. code-block::
PetscOptionsHasName(PetscOptions options, char *pre, char *name, PetscBool *flg);
PetscOptionsGetInt(PetscOptions options, char *pre, char *name, PetscInt *value, PetscBool *flg);
PetscOptionsGetReal(PetscOptions options, char *pre, char *name, PetscReal *value, PetscBool *flg);
PetscOptionsGetString(PetscOptions options, char *pre, char *name, char *value, size_t maxlen, PetscBool *flg);
PetscOptionsGetStringArray(PetscOptions options, char *pre, char *name, char **values, PetscInt *nmax, PetscBool *flg);
PetscOptionsGetIntArray(PetscOptions options, char *pre, char *name, PetscInt *value, PetscInt *nmax, PetscBool *flg);
PetscOptionsGetRealArray(PetscOptions options, char *pre, char *name, PetscReal *value, PetscInt *nmax, PetscBool *flg);
All of these routines set ``flg=PETSC_TRUE`` if the corresponding option
was found, ``flg=PETSC_FALSE`` if it was not found. The optional
argument ``pre`` indicates that the true name of the option is the given
name (with the dash “-” removed) prepended by the prefix ``pre``.
Usually ``pre`` should be set to ``NULL`` (or ``PETSC_NULL_CHARACTER``
for Fortran); its purpose is to allow someone to rename all the options
in a package without knowing the names of the individual options. For
example, when using block Jacobi preconditioning, the ``KSP`` and ``PC``
methods used on the individual blocks can be controlled via the options
``-sub_ksp_type`` and ``-sub_pc_type``.
Keeping Track of Options
^^^^^^^^^^^^^^^^^^^^^^^^
One useful means of keeping track of user-specified runtime options is
use of ``-options_view``, which prints to ``stdout`` during
``PetscFinalize()`` a table of all runtime options that the user has
specified. A related option is ``-options_left``, which prints the
options table and indicates any options that have *not* been requested
upon a call to ``PetscFinalize()``. This feature is useful to check
whether an option has been activated for a particular PETSc object (such
as a solver or matrix format), or whether an option name may have been
accidentally misspelled.
The option ``-options_monitor`` ``<viewer>`` turns on the default monitoring of options.
``PetscOptionsMonitorSet()`` can be used to provide custom monitors.
The option ``-options_monitor_cancel`` prevents any monitoring by monitors set with ``PetscOptionsMonitorSet()`` (but not that set with ``-options_monitor``).
.. _sec_viewers:
Viewers: Looking at PETSc Objects
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
PETSc employs a consistent scheme for examining, printing, and saving
objects through commands of the form
.. code-block::
XXXView(XXX obj, PetscViewer viewer);
Here ``obj`` is a PETSc object of type ``XXX``, where ``XXX`` is
``Mat``, ``Vec``, ``SNES``, etc. There are several predefined viewers.
- Passing in a zero (``0``) for the viewer causes the object to be
printed to the screen; this is useful when viewing an object in a
debugger but should be avoided in source code.
- ``PETSC_VIEWER_STDOUT_SELF`` and ``PETSC_VIEWER_STDOUT_WORLD`` causes
the object to be printed to the screen.
- ``PETSC_VIEWER_DRAW_SELF`` ``PETSC_VIEWER_DRAW_WORLD`` causes the
object to be drawn in a default X window.
- Passing in a viewer obtained by ``PetscViewerDrawOpen()`` causes the
object to be displayed graphically. See
:any:`sec_graphics` for more on PETSc’s graphics support.
- To save an object to a file in ASCII format, the user creates the
viewer object with the command
``PetscViewerASCIIOpen(MPI_Comm comm, char* file, PetscViewer *viewer)``.
This object is analogous to ``PETSC_VIEWER_STDOUT_SELF`` (for a
communicator of ``MPI_COMM_SELF``) and ``PETSC_VIEWER_STDOUT_WORLD``
(for a parallel communicator).
- To save an object to a file in binary format, the user creates the
viewer object with the command
``PetscViewerBinaryOpen(MPI_Comm comm, char* file, PetscViewerBinaryType type, PetscViewer *viewer)``.
Details of binary I/O are discussed below.
- Vector and matrix objects can be passed to a running MATLAB process
with a viewer created by
``PetscViewerSocketOpen(MPI_Comm comm, char *machine, int port, PetscViewer *viewer)``.
See :any:`sec_matlabsocket`.
The user can control the format of ASCII printed objects with viewers
created by ``PetscViewerASCIIOpen()`` by calling
.. code-block::
PetscViewerPushFormat(PetscViewer viewer, PetscViewerFormat format);
Formats include ``PETSC_VIEWER_DEFAULT``, ``PETSC_VIEWER_ASCII_MATLAB``,
and ``PETSC_VIEWER_ASCII_IMPL``. The implementation-specific format,
``PETSC_VIEWER_ASCII_IMPL``, displays the object in the most natural way
for a particular implementation.
The routines
.. code-block::
PetscViewerPushFormat(PetscViewer viewer, PetscViewerFormat format);
PetscViewerPopFormat(PetscViewer viewer);
allow one to temporarily change the format of a viewer.
As discussed above, one can output PETSc objects in binary format by
first opening a binary viewer with ``PetscViewerBinaryOpen()`` and then
using ``MatView()``, ``VecView()``, etc. The corresponding routines for
input of a binary object have the form ``XXXLoad()``. In particular,
matrix and vector binary input is handled by the following routines:
.. code-block::
MatLoad(Mat newmat, PetscViewer viewer);
VecLoad(Vec newvec, PetscViewer viewer);
These routines generate parallel matrices and vectors if the viewer’s
communicator has more than one process. The particular matrix and vector
formats are determined from the options database; see the manual pages
for details.
One can provide additional information about matrix data for matrices
stored on disk by providing an optional file ``matrixfilename.info``,
where ``matrixfilename`` is the name of the file containing the matrix.
The format of the optional file is the same as the ``.petscrc`` file and
can (currently) contain the following:
.. code-block:: none
-matload_block_size <bs>
The block size indicates the size of blocks to use if the matrix is read
into a block oriented data structure (for example, ``MATMPIBAIJ``). The
diagonal information ``s1,s2,s3,...`` indicates which (block) diagonals
in the matrix have nonzero values. The info file is automatically created
when ``VecView()`` or ``MatView()`` is used with a binary viewer; hence if you
save a matrix with a given block size with ``MatView()``, then a ``MatLoad()``
on that file will automatically use the saved block size.
.. _sec_viewfromoptions:
Viewing From Options
^^^^^^^^^^^^^^^^^^^^
Command-line options provide a particularly convenient way to view PETSc
objects. All options of the form ``-xxx_view`` accept
colon(``:``)-separated compound arguments which specify a viewer type,
format, and/or destination (e.g. file name or socket) if appropriate.
For example, to quickly export a binary file containing a matrix, one
may use ``-mat_view binary:matrix.out``, or to output to a
MATLAB-compatible ASCII file, one may use
``-mat_view ascii:matrix.m:ascii_matlab``. See the
``PetscOptionsCreateViewer()`` man page for full details, as well as the
``XXXViewFromOptions()`` man pages (for instance,
``PetscDrawSetFromOptions()``) for many other convenient command-line
options.
Using Viewers to Check Load Imbalance
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
The ``PetscViewer`` format ``PETSC_VIEWER_LOAD_BALANCE`` will cause certain
objects to display simple measures of their imbalance. For example
.. code-block:: none
-n 4 ./ex32 -ksp_view_mat ::load_balance
will display
.. code-block:: none
Nonzeros: Min 162 avg 168 max 174
indicating that one process has 162 nonzero entries in the matrix, the
average number of nonzeros per process is 168 and the maximum number of
nonzeros is 174. Similar for vectors one can see the load balancing
with, for example,
.. code-block:: none
-n 4 ./ex32 -ksp_view_rhs ::load_balance
The measurements of load balancing can also be done within the program
with calls to the appropriate object viewer with the viewer format
``PETSC_VIEWER_LOAD_BALANCE``.
.. _sec_saws:
Using SAWs with PETSc
~~~~~~~~~~~~~~~~~~~~~
The Scientific Application Web server, SAWs [#saws]_, allows one to monitor
running PETSc applications from a browser. To use SAWs you must ``configure`` PETSc with
the option ``--download-saws``. Options to use SAWs include
- ``-saws_options`` - allows setting values in the PETSc options
database via the browser (works only on one process).
- ``-stack_view saws`` - allows monitoring the current stack frame that
PETSc is in; refresh to see the new location.
- ``-snes_monitor_saws, -ksp_monitor_saws`` - monitor the solvers’
iterations from the web browser.
For each of these you need to point your browser to
``http://hostname:8080``, for example ``http://localhost:8080``. Options
that control behavior of SAWs include
- ``-saws_log filename`` - log all SAWs actions in a file.
- ``-saws_https certfile`` - use HTTPS instead of HTTP with a
certificate.
- ``-saws_port_auto_select`` - have SAWs pick a port number instead of
using 8080.
- ``-saws_port port`` - use ``port`` instead of 8080.
- ``-saws_root rootdirectory`` - local directory to which the SAWs
browser will have read access.
- ``-saws_local`` - use the local file system to obtain the SAWS
javascript files (they much be in ``rootdirectory/js``).
Also see the manual pages for ``PetscSAWsBlock()``,
``PetscObjectSAWsTakeAccess()``, ``PetscObjectSAWsGrantAccess()``,
``PetscObjectSAWsSetBlock()``, ``PetscStackSAWsGrantAccess()``
``PetscStackSAWsTakeAccess()``, ``KSPMonitorSAWs()``, and
``SNESMonitorSAWs()``.
.. _sec-debugging:
Debugging
~~~~~~~~~
PETSc programs may be debugged using one of the two options below.
- ``-start_in_debugger`` ``[noxterm,dbx,xxgdb,xdb,xldb,lldb]``
``[-display name]`` - start all processes in debugger
- ``-on_error_attach_debugger`` ``[noxterm,dbx,xxgdb,xdb,xldb,lldb]``
``[-display name]`` - start debugger only on encountering an error
Note that, in general, debugging MPI programs cannot be done in the
usual manner of starting the programming in the debugger (because then
it cannot set up the MPI communication and remote processes).
By default on Linux systems the GNU debugger ``gdb`` is used, on macOS systems ``lldb`` is used
By default, the debugger will be started in a new
xterm (Apple Terminal on macOS), to enable running separate debuggers on each process, unless the
option ``noxterm`` is used. In order to handle the MPI startup phase,
the debugger command ``cont`` should be used to continue execution of
the program within the debugger. Rerunning the program through the
debugger requires terminating the first job and restarting the
processor(s); the usual ``run`` option in the debugger will not
correctly handle the MPI startup and should not be used. Not all
debuggers work on all machines, the user may have to experiment to find
one that works correctly.
You can select a subset of the processes to be debugged (the rest just
run without the debugger) with the option
.. code-block:: none
-debugger_ranks rank1,rank2,...
where you simply list the ranks you want the debugger to run with.
.. _sec_errors:
Error Handling
~~~~~~~~~~~~~~
Errors are handled through the routine ``PetscError()``. This routine
checks a stack of error handlers and calls the one on the top. If the
stack is empty, it selects ``PetscTraceBackErrorHandler()``, which tries
to print a traceback. A new error handler can be put on the stack with
.. code-block::
PetscPushErrorHandler(PetscErrorCode (*HandlerFunction)(int line, char *dir, char *file, char *message, int number, void*), void *HandlerContext)
The arguments to ``HandlerFunction()`` are the line number where the
error occurred, the file in which the error was detected, the
corresponding directory, the error message, the error integer, and the
``HandlerContext.`` The routine
.. code-block::
PetscPopErrorHandler()
removes the last error handler and discards it.
PETSc provides two additional error handlers besides
``PetscTraceBackErrorHandler()``:
.. code-block::
PetscAbortErrorHandler()
PetscAttachDebuggerErrorHandler()
The function ``PetscAbortErrorHandler()`` calls abort on encountering an
error, while ``PetscAttachDebuggerErrorHandler()`` attaches a debugger to the
running process if an error is detected. At runtime, these error
handlers can be set with the options ``-on_error_abort`` or
``-on_error_attach_debugger`` ``[noxterm, dbx, xxgdb, xldb]``
``[-display DISPLAY]``.
All PETSc calls can be traced (useful for determining where a program is
hanging without running in the debugger) with the option
.. code-block:: none
-log_trace [filename]
where ``filename`` is optional. By default the traces are printed to the
screen. This can also be set with the command
``PetscLogTraceBegin(FILE*)``.
It is also possible to trap signals by using the command
.. code-block::
PetscPushSignalHandler(PetscErrorCode (*Handler)(int, void *), void *ctx);
The default handler ``PetscSignalHandlerDefault()`` calls
``PetscError()`` and then terminates. In general, a signal in PETSc
indicates a catastrophic failure. Any error handler that the user
provides should try to clean up only before exiting. By default all
PETSc programs turn on the default PETSc signal handler in ``PetscInitialize()``,
this can be prevented with the option ``-no_signal_handler`` that can be provided on the command line,
in the ~./petscrc file, or with the call
.. code-block::
PetscCall(PetscOptionsSetValue(NULL, "-no_signal_handler", "true"));
Once the first PETSc signal handler has been pushed it is impossible to go back to
to a signal handler that was set directly by the user with the UNIX signal handler API or by
the loader.
Some Fortran compilers/loaders cause, by default, a traceback of the Fortran call stack when a
segmentation violation occurs to be printed. This is handled by them setting a special signal handler
when the program is started up. This feature is useful for debugging without needing to start up a debugger.
If ``PetscPushSignalHandler()`` has been called this traceback will not occur, hence if the Fortran traceback
is desired one should put
.. code-block::
PetscCallA(PetscOptionsSetValue(PETSC_NULL_OPTIONS,"-no_signal_handler","true",ierr))
**before** the call to ``PetscInitialize()``. This prevents PETSc from defaulting to using a signal handler.
There is a separate signal handler for floating-point exceptions. The
option ``-fp_trap`` turns on the floating-point trap at runtime, and the
routine
.. code-block::
PetscFPTrapPush(PetscFPTrap flag);
can be used within a program. A ``flag`` of ``PETSC_FP_TRAP_ON`` indicates that
floating-point exceptions should be trapped, while a value of
``PETSC_FP_TRAP_OFF`` (the default) indicates that they should be
ignored.
.. code-block::
PetscFPTrapPop(void);
should be used to revert to the previous handling of floating point exceptions before the call to ``PetscFPTrapPush()``.
A small set of macros is used to make the error handling lightweight.
These macros are used throughout the PETSc libraries and can be employed
by the application programmer as well. When an error is first detected,
one should set it by calling
.. code-block::
SETERRQ(MPI_Comm comm, PetscErrorCode flag, char *message);
The user should check the return codes for all PETSc routines (and
possibly user-defined routines as well) with
.. code-block::
PetscCall(PetscRoutine(...));
Likewise, all memory allocations should be checked with
.. code-block::
PetscCall(PetscMalloc1(n, &ptr));
If this procedure is followed throughout all of the user’s libraries and
codes, any error will by default generate a clean traceback of the
location of the error.
Note that the macro ``PETSC_FUNCTION_NAME`` is used to keep track of
routine names during error tracebacks. Users need not worry about this
macro in their application codes; however, users can take advantage of
this feature if desired by setting this macro before each user-defined
routine that may call ``SETERRQ()``, ``PetscCall()``. A simple example of
usage is given below.
.. code-block::
PetscErrorCode MyRoutine1()
{
/* Declarations Here */
PetscFunctionBeginUser;
/* code here */
PetscFunctionReturn(PETSC_SUCCESS);
}
.. _sec_complex:
Numbers
~~~~~~~
PETSc supports the use of complex numbers in application programs
written in C, C++, and Fortran. To do so, we employ either the C99
``complex`` type or the C++ versions of the PETSc libraries in which the
basic “scalar” datatype, given in PETSc codes by ``PetscScalar``, is
defined as ``complex`` (or ``complex<double>`` for machines using
templated complex class libraries). To work with complex numbers, the
user should run ``configure`` with the additional option
``--with-scalar-type=complex``. The
:doc:`installation instructions </install/index>`
provide detailed instructions for installing PETSc. You can use
``--with-clanguage=c`` (the default) to use the C99 complex numbers or
``--with-clanguage=c++`` to use the C++ complex type [#cxx_note]_.
Recall that each configuration of the PETSc libraries is stored in a different
directory, given by ``$PETSC_DIR/$PETSC_ARCH``
according to the architecture. Thus, the libraries for complex numbers
are maintained separately from those for real numbers. When using any of
the complex numbers versions of PETSc, *all* vector and matrix elements
are treated as complex, even if their imaginary components are zero. Of
course, one can elect to use only the real parts of the complex numbers
when using the complex versions of the PETSc libraries; however, when
working *only* with real numbers in a code, one should use a version of
PETSc for real numbers for best efficiency.
The program
`KSP Tutorial ex11 <PETSC_DOC_OUT_ROOT_PLACEHOLDER/src/ksp/ksp/tutorials/ex11.c.html>`__
solves a linear system with a complex coefficient matrix. Its Fortran
counterpart is
`KSP Tutorial ex11f <PETSC_DOC_OUT_ROOT_PLACEHOLDER/src/ksp/ksp/tutorials/ex11f.F90.html>`__.
Parallel Communication
~~~~~~~~~~~~~~~~~~~~~~
When used in a message-passing environment, all communication within
PETSc is done through MPI, the message-passing interface standard
:cite:`mpi-final`. Any file that includes ``petscsys.h`` (or
any other PETSc include file) can freely use any MPI routine.
.. _sec_graphics:
Graphics
~~~~~~~~
The PETSc graphics library is not intended to compete with high-quality
graphics packages. Instead, it is intended to be easy to use
interactively with PETSc programs. We urge users to generate their
publication-quality graphics using a professional graphics package. If a
user wants to hook certain packages into PETSc, he or she should send a
message to
`petsc-maint@mcs.anl.gov <mailto:petsc-maint@mcs.anl.gov>`__; we
will see whether it is reasonable to try to provide direct interfaces.
Windows as PetscViewers
^^^^^^^^^^^^^^^^^^^^^^^
For drawing predefined PETSc objects such as matrices and vectors, one
may first create a viewer using the command
.. code-block::
PetscViewerDrawOpen(MPI_Comm comm, char *display, char *title, int x, int y, int w, int h, PetscViewer *viewer);
This viewer may be passed to any of the ``XXXView()`` routines.
Alternately, one may use command-line options to quickly specify viewer
formats, including ``PetscDraw``-based ones; see
:any:`sec_viewfromoptions`.
To draw directly into the viewer, one must obtain the ``PetscDraw``
object with the command
.. code-block::
PetscViewerDrawGetDraw(PetscViewer viewer, PetscDraw *draw);
Then one can call any of the ``PetscDrawXXX()`` commands on the ``draw``
object. If one obtains the ``draw`` object in this manner, one does not
call the ``PetscDrawOpenX()`` command discussed below.
Predefined viewers, ``PETSC_VIEWER_DRAW_WORLD`` and
``PETSC_VIEWER_DRAW_SELF``, may be used at any time. Their initial use
will cause the appropriate window to be created.
Implementations using OpenGL, TikZ, and other formats may be selected
with ``PetscDrawSetType()``. PETSc can also produce movies; see
``PetscDrawSetSaveMovie()``, and note that command-line options can also
be convenient; see the ``PetscDrawSetFromOptions()`` man page.
By default, PETSc drawing tools employ a private colormap, which
remedies the problem of poor color choices for contour plots due to an
external program’s mangling of the colormap. Unfortunately, this may
cause flashing of colors as the mouse is moved between the PETSc windows
and other windows. Alternatively, a shared colormap can be used via the
option ``-draw_x_shared_colormap``.
Simple PetscDrawing
^^^^^^^^^^^^^^^^^^^
With the default format, one can open a window that is not associated
with a viewer directly under the X11 Window System with the
command
.. code-block::
PetscDrawCreate(MPI_Comm comm, char *display, char *title, int x, int y, int w, int h, PetscDraw *win);
PetscDrawSetFromOptions(win);
All drawing routines are performed relative to the window’s coordinate
system and viewport. By default, the drawing coordinates are from
``(0,0)`` to ``(1,1)``, where ``(0,0)`` indicates the lower left corner
of the window. The application program can change the window coordinates
with the command
.. code-block::
PetscDrawSetCoordinates(PetscDraw win, PetscReal xl, PetscReal yl, PetscReal xr, PetscReal yr);
By default, graphics will be drawn in the entire window. To restrict the
drawing to a portion of the window, one may use the command
.. code-block::
PetscDrawSetViewPort(PetscDraw win, PetscReal xl, PetscReal yl, PetscReal xr, PetscReal yr);
These arguments, which indicate the fraction of the window in which the
drawing should be done, must satisfy
:math:`0 \leq {\tt xl} \leq {\tt xr} \leq 1` and
:math:`0 \leq {\tt yl} \leq {\tt yr} \leq 1.`
To draw a line, one uses the command
.. code-block::
PetscDrawLine(PetscDraw win, PetscReal xl, PetscReal yl, PetscReal xr, PetscReal yr, int cl);
The argument ``cl`` indicates the color (which is an integer between 0
and 255) of the line. A list of predefined colors may be found in
``include/petscdraw.h`` and includes ``PETSC_DRAW_BLACK``,
``PETSC_DRAW_RED``, ``PETSC_DRAW_BLUE`` etc.
To ensure that all graphics actually have been displayed, one should use
the command
.. code-block::
PetscDrawFlush(PetscDraw win);
When displaying by using double buffering, which is set with the command
.. code-block::
PetscDrawSetDoubleBuffer(PetscDraw win);
*all* processes must call
.. code-block::
PetscDrawFlush(PetscDraw win);
in order to swap the buffers. From the options database one may use
``-draw_pause`` ``n``, which causes the PETSc application to pause ``n``
seconds at each ``PetscDrawPause()``. A time of ``-1`` indicates that
the application should pause until receiving mouse input from the user.
Text can be drawn with commands
.. code-block::
PetscDrawString(PetscDraw win, PetscReal x, PetscReal y, int color, char *text);
PetscDrawStringVertical(PetscDraw win, PetscReal x, PetscReal y, int color, const char *text);
PetscDrawStringCentered(PetscDraw win, PetscReal x, PetscReal y, int color, const char *text);
PetscDrawStringBoxed(PetscDraw draw, PetscReal sxl, PetscReal syl, int sc, int bc, const char text[], PetscReal *w, PetscReal *h);
The user can set the text font size or determine it with the commands
.. code-block::
PetscDrawStringSetSize(PetscDraw win, PetscReal width, PetscReal height);
PetscDrawStringGetSize(PetscDraw win, PetscReal *width, PetscReal *height);
Line Graphs
^^^^^^^^^^^
PETSc includes a set of routines for manipulating simple two-dimensional
graphs. These routines, which begin with ``PetscDrawAxisDraw()``, are
usually not used directly by the application programmer. Instead, the
programmer employs the line graph routines to draw simple line graphs.
As shown in the :ref:`listing below <listing_draw_test_ex3>`, line
graphs are created with the command
.. code-block::
PetscDrawLGCreate(PetscDraw win, PetscInt ncurves, PetscDrawLG *ctx);
The argument ``ncurves`` indicates how many curves are to be drawn.
Points can be added to each of the curves with the command
.. code-block::
PetscDrawLGAddPoint(PetscDrawLG ctx, PetscReal *x, PetscReal *y);
The arguments ``x`` and ``y`` are arrays containing the next point value
for each curve. Several points for each curve may be added with
.. code-block::
PetscDrawLGAddPoints(PetscDrawLG ctx, PetscInt n, PetscReal **x, PetscReal **y);
The line graph is drawn (or redrawn) with the command
.. code-block::
PetscDrawLGDraw(PetscDrawLG ctx);
A line graph that is no longer needed can be destroyed with the command
.. code-block::
PetscDrawLGDestroy(PetscDrawLG *ctx);
To plot new curves, one can reset a linegraph with the command
.. code-block::
PetscDrawLGReset(PetscDrawLG ctx);
The line graph automatically determines the range of values to display
on the two axes. The user can change these defaults with the command
.. code-block::
PetscDrawLGSetLimits(PetscDrawLG ctx, PetscReal xmin, PetscReal xmax, PetscReal ymin, PetscReal ymax);
It is also possible to change the display of the axes and to label them.
This procedure is done by first obtaining the axes context with the
command
.. code-block::
PetscDrawLGGetAxis(PetscDrawLG ctx, PetscDrawAxis *axis);
One can set the axes’ colors and labels, respectively, by using the
commands
.. code-block::
PetscDrawAxisSetColors(PetscDrawAxis axis, int axis_lines, int ticks, int text);
PetscDrawAxisSetLabels(PetscDrawAxis axis, char *top, char *x, char *y);
It is possible to turn off all graphics with the option ``-nox``. This
will prevent any windows from being opened or any drawing actions to be
done. This is useful for running large jobs when the graphics overhead
is too large, or for timing.
The full example, `Draw Test ex3 <PETSC_DOC_OUT_ROOT_PLACEHOLDER/src/sys/classes/draw/tests/ex3.c.html>`__,
follows.
.. _listing_draw_test_ex3:
.. admonition:: Listing: ``src/classes/draw/tests/ex3.c``
:name: snes-ex1
.. literalinclude:: /../src/sys/classes/draw/tests/ex3.c
:end-before: /*TEST
Graphical Convergence Monitor
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
For both the linear and nonlinear solvers default routines allow one to
graphically monitor convergence of the iterative method. These are
accessed via the command line with ``-ksp_monitor draw::draw_lg`` and
``-snes_monitor draw::draw_lg``. See also
:any:`sec_kspmonitor` and :any:`sec_snesmonitor`.
Disabling Graphics at Compile Time
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
To disable all X-window-based graphics, run ``configure`` with the
additional option ``--with-x=0``
.. _sec-developer-environments:
Developer Environments
----------------------
Emacs Users
~~~~~~~~~~~
Many PETSc developers use Emacs, which can be used as a "simple" text editor or a comprehensive development environment.
For a more integrated development environment, we recommend using `lsp-mode <https://emacs-lsp.github.io/lsp-mode/>`_ (or `eglot <https://github.com/joaotavora/eglot>`_) with `clangd <https://clangd.llvm.org/>`_.
The most convenient way to teach clangd what compilation flags to use is to install `Bear <https://github.com/rizsotto/Bear>`_ ("build ear") and run::
bear make -B
which will do a complete rebuild (``-B``) of PETSc and capture the compilation commands in a file named ``compile_commands.json``, which will be automatically picked up by clangd.
You can use the same procedure when building examples or your own project.
It can also be used with any other editor that supports clangd, including VS Code and Vim.
When lsp-mode is accompanied by `flycheck <https://www.flycheck.org/en/latest/>`_, Emacs will provide real-time feedback and syntax checking, along with refactoring tools provided by clangd.
The easiest way to install packages in recent Emacs is to use the "Options" menu to select "Manage Emacs Packages".
Tags
^^^^
It is sometimes useful to cross-reference tags across projects.
Regardless of whether you use lsp-mode, it can be useful to use `GNU Global <https://www.gnu.org/software/global/>`_ (install ``gtags``) to provide reverse lookups (e.g. find all call sites
for a given function) across all projects you might work on/browse.
Tags for PETSc can be generated by running ``make allgtags`` from ``$PETSC_DIR``, or one can generate tags for all projects by running a command such as
.. code-block:: none
find $PETSC_DIR/{include,src,tutorials,$PETSC_ARCH/include} any/other/paths \
-regex '.*\.\(cc\|hh\|cpp\|cxx\|C\|hpp\|c\|h\|cu\)$' \
| grep -v ftn-auto | gtags -f -
from your home directory or wherever you keep source code. If you are
making large changes, it is useful to either set this up to run as a
cron job or to make a convenient alias so that refreshing is easy. Then
add the following to ``~/.emacs`` to enable gtags and specify key bindings.
.. code-block:: emacs
(when (require 'gtags)
(global-set-key (kbd "C-c f") 'gtags-find-file)
(global-set-key (kbd "C-c .") 'gtags-find-tag)
(global-set-key (kbd "C-c r") 'gtags-find-rtag)
(global-set-key (kbd "C-c ,") 'gtags-pop-stack))
(add-hook 'c-mode-common-hook
'(lambda () (gtags-mode t))) ; Or add to existing hook
A more basic alternative to the GNU Global (``gtags``) approach that does not require adding packages is to use
the builtin ``etags`` feature. First, run ``make alletags`` from the
PETSc home directory to generate the file ``$PETSC_DIR/TAGS``, and
then from within Emacs, run
.. code-block:: none
M-x visit-tags-table
where ``M`` denotes the Emacs Meta key, and enter the name of the
``TAGS`` file. Then the command ``M-.`` will cause Emacs to find the
file and line number where a desired PETSc function is defined. Any
string in any of the PETSc files can be found with the command ``M-x tags-search``.
To find repeated occurrences, one can simply use ``M-,`` to find the next occurrence.
VS Code Users
~~~~~~~~~~~~~
`VS Code <https://code.visualstudio.com/>`_ (unlike :ref:`sec-visual-studio`, described below) is an open-source editor with a rich extension ecosystem.
It has `excellent integration <https://marketplace.visualstudio.com/items?itemName=llvm-vs-code-extensions.vscode-clangd>`_ with clangd and will automatically pick up ``compile_commands.json``
as produced by a command such as ``bear make -B`` (see :ref:`sec-developer-environments`).
If you have no prior attachment to a specific code editor, we recommend trying VS Code.
Vi and Vim Users
~~~~~~~~~~~~~~~~
This section lists helpful Vim commands for PETSc. Ones that configure Vim can be placed
in a ``.vimrc`` file in the top of the PETSc directory and will be loaded automatically.
Vim has configurable keymaps: all of the "command mode" commands given that start with
a colon (such as ``:help``) can be assigned to short sequences in "normal mode," which
is how most Vim users use their most frequently used commands.
See the :ref:`sec-developer-environments` discussion above for configuration of clangd, which
provides integrated development environment.
Tags
^^^^
The ``tags`` feature can be used to search PETSc files quickly and efficiently.
To use this feature, one should first check if the file, ``$PETSC_DIR/CTAGS``
exists. If this file is not present, it should be generated by running ``make
alletags`` from the PETSc home directory. Once the file exists, from Vi/Vim the
user should issue the command
.. code-block:: none
:set tags=CTAGS
from the ``$PETSC_DIR`` directory and enter the name of the ``CTAGS`` file. The
command ``:tag functionname`` will cause Vi/Vim to open the file and line
number where a desired PETSc function is defined in the current window.
``<Ctrl-o>`` will return the screen to your previous location.
The command ``:stag functionname`` will split the current window and then open
the file and line number for that function in one half. Some prefer this because
it is easier to compare the file you are editing to the function definition this way.
Cscope and gtags
^^^^^^^^^^^^^^^^
Vim can also use the ``cscope`` utility to navigate source code. One useful thing
it can do that the basic ``tags`` feature can't is search for references to a symbol,
rather than its definition, which can be useful for refactoring. The command
.. code-block:: none
:cs find s functionname
opens a list of all of the places the function is called in PETSc, and opens
the file and line that you choose. The variant ``:scs find s functionname``
does the same but splits the window like ``stag``.
The PETSc makefile does not have a command for building a cscope database, but
GNU Global is cross-compatible with cscope: call ``make allgtags`` to make the
gtags database, and run the commands
.. code-block:: none
:set csprg=gtags-cscope
:cs add GTAGS
Quickfix
^^^^^^^^
Rather than exiting editing a file to build the library and check for errors or
warnings, calling ``:make`` runs the make command without leaving Vim and
collects the errors and warnings in a "quickfix" window. Move the cursor to
one of the errors or warnings in the quickfix window and press ``<Enter>`` and
the main window will jump to the file and line with the error. The following
commands filter lines of out PETSc's make output that can clutter the quickfix window:
.. code-block:: none
:set efm^=%-GStarting\ make\ run\ on\ %.%#
:set efm^=%-GMachine\ characteristics:\ %.%#
:set efm^=%-G#define\ PETSC%.%#
Autocompletion and snippets
^^^^^^^^^^^^^^^^^^^^^^^^^^^
Autocompletion of long function names can be helpful when working with PETSc.
If you have a tags file, you can press ``<Ctrl-N>`` when you have partially
typed a word to bring up a list of potential completions that you can choose
from with ``<Tab>``.
More powerful autocompletion, such as completing the fieldname of a struct, is
available from external plugins that can be added to Vim, such as `SuperTab
<https://github.com/ervandew/supertab>`__, `VimCompletesMe
<https://github.com/ackyshake/VimCompletesMe>`__, or `YouCompleteMe
<https://github.com/ycm-core/YouCompleteMe>`__.
Along the same lines, plugins can be added that fill in the boilerplate
associated with PETSc programming with code snippets. One such tool is
`UltiSnips <https://github.com/sirver/UltiSnips>`__.
LSP for Vim
^^^^^^^^^^^
Several plugins provide the equivalent of emacs' lsp-mode: YouCompleteMe,
mentioned above, is one; another popular one is `ale
<https://github.com/dense-analysis/ale>`__. These can check for syntax errors,
check for compilation errors in the background, and provide sophisticated tools
for refactoring. Like lsp-mode, they also rely on a compilation database, so
``bear -- make -B`` should be used as well to generate the file
``compile_commands.json``.
See `online tutorials <http://www.yolinux.com/TUTORIALS/LinuxTutorialAdvanced_vi.html>`__
for additional Vi/Vim options.
Eclipse Users
~~~~~~~~~~~~~
If you are interested in developing code that uses PETSc from Eclipse or
developing PETSc in Eclipse and have knowledge of how to do indexing and
build libraries in Eclipse, please contact us at
`petsc-dev@mcs.anl.gov <mailto:petsc-dev@mcs.anl.gov>`_.
One way to index and build PETSc in Eclipse is as follows.
#. Open
“File\ :math:`\rightarrow`\ Import\ :math:`\rightarrow`\ Git\ :math:`\rightarrow`\ Projects
from Git”. In the next two panels, you can either add your existing
local repository or download PETSc from Bitbucket by providing the
URL. Most Eclipse distributions come with Git support. If not,
install the EGit plugin. When importing the project, select the
wizard “Import as general project”.
#. Right-click on the project (or the “File” menu on top) and select
“New :math:`\rightarrow` Convert to a C/C++ Project (Adds C/C++
Nature)”. In the setting window, choose “C Project” and specify the
project type as “Shared Library”.
#. Right-click on the C project and open the “Properties” panel. Under
“C/C++ Build :math:`\rightarrow` Builder Settings”, set the Build
directory to ``$PETSC_DIR`` and make sure “Generate Makefiles
automatically” is unselected. Under the section “C/C++
General\ :math:`\rightarrow`\ Paths and Symbols”, add the PETSc paths
to “Includes”.
.. code-block:: none
$PETSC_DIR/include
$PETSC_DIR/$PETSC_ARCH/include
Under the section “C/C++ General\ :math:`\rightarrow`\ index”, choose
“Use active build configuration”.
#. Configure PETSc normally outside Eclipse to generate a makefile and
then build the project in Eclipse. The source code will be parsed by
Eclipse.
If you launch Eclipse from the Dock on Mac OS X, ``.bashrc`` will not be
loaded (a known OS X behavior, for security reasons). This will be a
problem if you set the environment variables ``$PETSC_DIR`` and
``$PETSC_ARCH`` in ``.bashrc``. A solution which involves replacing the
executable can be found at
```/questions/829749/launch-mac-eclipse-with-environment-variables-set`` </questions/829749/launch-mac-eclipse-with-environment-variables-set>`__.
Alternatively, you can add ``$PETSC_DIR`` and ``$PETSC_ARCH`` manually
under “Properties :math:`\rightarrow` C/C++ Build :math:`\rightarrow`
Environment”.
To allow an Eclipse code to compile with the PETSc include files and
link with the PETSc libraries, a PETSc user has suggested the following.
#. Right-click on your C project and select “Properties
:math:`\rightarrow` C/C++ Build :math:`\rightarrow` Settings”
#. A new window on the righthand side appears with various settings
options. Select “Includes” and add the required PETSc paths,
.. code-block:: none
$PETSC_DIR/include
$PETSC_DIR/$PETSC_ARCH/include
#. Select “Libraries” under the header Linker and set the library search
path:
.. code-block:: none
$PETSC_DIR/$PETSC_ARCH/lib
and the libraries, for example
.. code-block:: none
m, petsc, stdc++, mpichxx, mpich, lapack, blas, gfortran, dl, rt,gcc_s, pthread, X11
Another PETSc user has provided the following steps to build an Eclipse
index for PETSc that can be used with their own code, without compiling
PETSc source into their project.
#. In the user project source directory, create a symlink to the PETSC
``src/`` directory.
#. Refresh the project explorer in Eclipse, so the new symlink is
followed.
#. Right-click on the project in the project explorer, and choose “Index
:math:`\rightarrow` Rebuild”. The index should now be built.
#. Right-click on the PETSc symlink in the project explorer, and choose
“Exclude from build...” to make sure Eclipse does not try to compile
PETSc with the project.
For further examples of using Eclipse with a PETSc-based application,
see the documentation for LaMEM [#lamem]_.
Qt Creator Users
~~~~~~~~~~~~~~~~
This information was provided by Mohammad Mirzadeh. The Qt Creator IDE
is part of the Qt SDK, developed for cross-platform GUI programming
using C++. It is available under GPL v3, LGPL v2 and a commercial
license and may be obtained, either as part of the Qt SDK or as
stand-alone software. It supports
automatic makefile generation using cross-platform ``qmake`` and
CMake build systems as well as allowing one to import projects based
on existing, possibly hand-written, makefiles. Qt Creator has a visual
debugger using GDB and LLDB (on Linux and OS X) or Microsoft’s CDB (on
Microsoft Windows) as backends. It also has an interface to Valgrind’s “memcheck”
and “callgrind” tools to detect memory leaks and profile code. It has
built-in support for a variety of version control systems including git,
mercurial, and subversion. Finally, Qt Creator comes fully equipped with
auto-completion, function look-up, and code refactoring tools. This
enables one to easily browse source files, find relevant functions, and
refactor them across an entire project.
Creating a Project
^^^^^^^^^^^^^^^^^^
When using Qt Creator with ``qmake``, one needs a ``.pro`` file. This
configuration file tells Qt Creator about all build/compile options and
locations of source files. One may start with a blank ``.pro`` file and
fill in configuration options as needed. For example:
.. code-block:: none
# The name of the application executable
TARGET = ex1
# There are two ways to add PETSc functionality
# 1-Manual: Set all include path and libs required by PETSc
PETSC_INCLUDE = path/to/petsc_includes # e.g. obtained via running `make getincludedirs'
PETSC_LIBS = path/to/petsc_libs # e.g. obtained via running `make getlinklibs'
INCLUDEPATH += $$PETSC_INCLUDES
LIBS += $$PETSC_LIBS
# 2-Automatic: Use the PKGCONFIG functionality
# NOTE: petsc.pc must be in the pkgconfig path. You might need to adjust PKG_CONFIG_PATH
CONFIG += link_pkgconfig
PKGCONFIG += PETSc
# Set appropriate compiler and its flags
QMAKE_CC = path/to/mpicc
QMAKE_CXX = path/to/mpicxx # if this is a cpp project
QMAKE_LINK = path/to/mpicxx # if this is a cpp project
QMAKE_CFLAGS += -O3 # add extra flags here
QMAKE_CXXFLAGS += -O3
QMAKE_LFLAGS += -O3
# Add all files that must be compiled
SOURCES += ex1.c source1.c source2.cpp
HEADERS += source1.h source2.h
# OTHER_FILES are ignored during compilation but will be shown in file panel in Qt Creator
OTHER_FILES += \
path/to/resource_file \
path/to/another_file
In this example, keywords include:
- ``TARGET``: The name of the application executable.
- ``INCLUDEPATH``: Used at compile time to point to required include
files. Essentially, it is used as an ``-I \$\$INCLUDEPATH`` flag for
the compiler. This should include all application-specific header
files and those related to PETSc (which may be found via
``make getincludedirs``).
- ``LIBS``: Defines all required external libraries to link with the
application. To get PETSc’s linking libraries, use
``make getlinklibs``.
- ``CONFIG``: Configuration options to be used by ``qmake``. Here, the
option ``link_pkgconfig`` instructs ``qmake`` to internally use
``pkgconfig`` to resolve ``INCLUDEPATH`` and ``LIBS`` variables.
- ``PKGCONFIG``: Name of the configuration file (the ``.pc`` file –
here ``petsc.pc``) to be passed to ``pkgconfig``. Note that for this
functionality to work, ``petsc.pc`` must be in path which might
require adjusting the ``PKG_CONFIG_PATH`` environment variable. For
more information see
`the Qt Creator documentation <https://doc.qt.io/qtcreator/creator-build-settings.html>`__.
- ``QMAKE_CC`` and ``QMAKE_CXX``: Define which C/C++ compilers use.
- ``QMAKE_LINK``: Defines the proper linker to be used. Relevant if
compiling C++ projects.
- ``QMAKE_CFLAGS``, ``QMAKE_CXXFLAGS`` and ``QMAKE_LFLAGS``: Set the
corresponding compile and linking flags.
- ``SOURCES``: Source files to be compiled.
- ``HEADERS``: Header files required by the application.
- ``OTHER_FILES``: Other files to include (source, header, or any other
extension). Note that none of the source files placed here are
compiled.
More options can be included in a ``.pro`` file; see
https://doc.qt.io/qt-5/qmake-project-files.html. Once the ``.pro`` file
is generated, the user can simply open it via Qt Creator. Upon opening,
one has the option to create two different build options, debug and
release, and switch between the two. For more information on using the
Qt Creator interface and other more advanced aspects of the IDE, refer
to https://www.qt.io/qt-features-libraries-apis-tools-and-ide/
.. _sec-visual-studio:
Visual Studio Users
~~~~~~~~~~~~~~~~~~~
To use PETSc from Microsoft Visual Studio, one would have to compile a PETSc
example with its corresponding makefile and then transcribe all compiler
and linker options used in this build into a Visual Studio project file,
in the appropriate format in Visual Studio project settings.
Xcode IDE Users
~~~~~~~~~~~~~~~
See :any:`doc_macos_install` for the standard Unix command line tools approach to development on macOS. The information
below is only if you plan to write code within the Xcode IDE.
Apple Xcode IDE for macOS Applications
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Follow the instructions in ``$PETSC_DIR/systems/Apple/OSX/bin/makeall``
to build the PETSc framework and documentation suitable for use in
Xcode.
You can then use the PETSc framework in
``$PETSC_DIR/arch-osx/PETSc.framework`` in the usual manner for Apple
frameworks. See the examples in
``$PETSC_DIR/systems/Apple/OSX/examples``. When working in Xcode, things
like function name completion should work for all PETSc functions as
well as MPI functions. You must also link against the Apple
``Accelerate.framework``.
Apple Xcode IDE for iPhone/iPad iOS Applications
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Follow the instructions in
``$PETSC_DIR/systems/Apple/iOS/bin/iosbuilder.py`` to build the PETSc
library for use on the iPhone/iPad.
You can then use the PETSc static library in
``$PETSC_DIR/arch-osx/libPETSc.a`` in the usual manner for Apple
libraries inside your iOS XCode projects; see the examples in
``$PETSC_DIR/systems/Apple/iOS/examples``. You must also link against
the Apple ``Accelerate.framework``.
A thorough discussion of the
procedure is given in `Comparison of Migration Techniques for High-Performance Code to Android and iOS
<https://www.researchgate.net/publication/308973080_Comparison_of_Migration_Techniques_for_High-Performance_Code_to_Android_and_iOS>`__.
For Android, you must have your standalone bin folder in the path, so that the compilers
are visible.
The installation process has not been tested for iOS or Android since 2017.
.. rubric:: Footnotes
.. [#saws]
`Saws wiki on Bitbucket <https://bitbucket.org/saws/saws/wiki/Home>`__
.. [#cxx_note]
Note that this option is not required to use PETSc with C++
.. [#lamem]
See the ``doc/`` directory at https://bitbucket.org/bkaus/lamem
.. [#json]
JSON is a subset of YAML
|