| #
a492ac7d
|
| 22-Sep-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
bug fixes for MatRelax_AIJ with flop count when no zero initial guess and passing PETSC_NULL when should pass in 1. (This is already fixed in petsc-dev but should be fixed in 3.0.0
Hg-commit: 70028a
bug fixes for MatRelax_AIJ with flop count when no zero initial guess and passing PETSC_NULL when should pass in 1. (This is already fixed in petsc-dev but should be fixed in 3.0.0
Hg-commit: 70028ad51299b48effba84a6e0a427b77d04d800
show more ...
|
| #
887ee2ca
|
| 01-Sep-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
fixed Eisentat preconditioner for AIJ and SBAIJ to work for all omega, instead of just omega== 1
Hg-commit: d40d325f2fcbde49d908cb600d0aae362304a789
|
| #
9d124a3d
|
| 25-Aug-2009 |
Matthew Knepley <knepley@mcs.anl.gov> |
Merge
Hg-commit: 9346f4fc0b2ed71af84b78752b7cd19df6d53a22
|
| #
ea7ccfb4
|
| 24-Aug-2009 |
Hong Zhang <hzhang@mcs.anl.gov> |
merge
Hg-commit: b384e8625b6324308f70798b433de65795605cb9
|
| #
a2b30743
|
| 20-Aug-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
added support for nonzero initial guess for Eisentat preconditioner; what UNIC needs and I forgot to put in before.
Hg-commit: c320a1da02e44812ec1e106a6b5c452116794b2e
|
| #
0c4e8012
|
| 19-Aug-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
removed Contributed by: in documentation and source for PETSc developers
Hg-commit: 19546b65bd4769308070770659a516cb6b59906a
|
| #
1e480f41
|
| 16-Aug-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
set maxnz to zero in MatLoad_MPIAIJ() to prevent a buggy compiler from complaining that it might be used uninitialized
Hg-commit: 5a81ba3fd8a41c320b07697ffd0bdeec18fd2d62
|
| #
7b55108e
|
| 14-Aug-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
print ' ' for default Prefetch macro because '' causes the BuildSystem to not print anything for the macro so PETSC_Prefetch() becomes undefined symbol
Hg-commit: bc8ef2a4dfd7f7982cd51f6417d5fe7edc2
print ' ' for default Prefetch macro because '' causes the BuildSystem to not print anything for the macro so PETSC_Prefetch() becomes undefined symbol
Hg-commit: bc8ef2a4dfd7f7982cd51f6417d5fe7edc2f0d94
show more ...
|
| #
df1bf9be
|
| 13-Aug-2009 |
Matthew Knepley <knepley@mcs.anl.gov> |
Merge
Hg-commit: b86fe0a90f488cd8a57ad47660550b9b828135b1
|
| #
effcda25
|
| 13-Aug-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
minor optimizations to MatRelax_MPISBAIJ()
removed global length checks on some vector operations so that one of the two vectors can be sequential
Hg-commit: 2a706c91dfbb7e840e89f07b1e42608aec372742
|
| #
0230a03c
|
| 31-Jul-2009 |
Matthew Knepley <knepley@mcs.anl.gov> |
Merge
Hg-commit: 349518ceeca84623e101fa73fa2e7aa9b404d154
|
| #
f6d58c54
|
| 25-Jul-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
finished support for MatFDColoringCreate() etc for BAIJ matrices
Hg-commit: 61c9609b5ffa140392733ce982b1260ca27cbd92
|
| #
3acb8795
|
| 23-Jul-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
more work on adding FD computation of Jacobian for BAIJ matrices
Hg-commit: 89075e85186765069f32cddc814471a4b6a51957
|
| #
952b491a
|
| 20-Jul-2009 |
Matthew Knepley <knepley@mcs.anl.gov> |
Merge
Hg-commit: 042472a0480b67a718f2ad368cd38c5dbee71cec
|
| #
907d6cfa
|
| 20-Jul-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
fixed bug introdued with last bug fix related to C++, PetscTruth and bools
Hg-commit: 50e52fe1bfb25b17ff847d83f791aab7bb49635f
|
| #
2767d60e
|
| 20-Jul-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
commit after merge
Hg-commit: 3130227ca1db4e998011a05bda2522c634b70f5b
|
| #
3b024144
|
| 18-Jul-2009 |
Hong Zhang <hzhang@mcs.anl.gov> |
casting
Hg-commit: d2a7d6d80cb5e8817053b6b6a100e94e23e6a8a0
|
| #
ae3dafa0
|
| 17-Jul-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
pull bug fix for The situation is the following: I'm calling MatSetValues with the flag ADD_VALUES and with matrix option MAT_IGNORE_ZERO_ENTRIES. I'm inserting multiple values at once, multiple colu
pull bug fix for The situation is the following: I'm calling MatSetValues with the flag ADD_VALUES and with matrix option MAT_IGNORE_ZERO_ENTRIES. I'm inserting multiple values at once, multiple columns and rows, so I provide a rank-2 matrix of values. As I'm calling this from fortran I'm also using MAT_COLUMN_ORIENTED. Now for provided rows that are not owned by the process, it jumps to mpiaij.c:394 (line numbers as in petsc-dev). On line 399, it checks for zero entries, but only checks the very first entry of the (non-owned) row. If however other entries of that same row are nonzero, the entire row is still dropped. Note that this is independent of row_oriented/column_oriented as line 396 does exactly the same.
If I don't set the option MAT_IGNORE_ZERO_ENTRIES the problem disappears. In that case however we would either have to preallocate substantially more nonzeros, or complicate the matrix assembly in our code by taking out the zero entries ourselves and call MatSetValues for each entry seperately.
Hg-commit: f09f0ab443e582d40fb2b646b2461bcc61f6e436
show more ...
|
| #
b400d20c
|
| 17-Jul-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
From: s.kramer@imperial.ac.uk
Subject: non-local values being dropped in MatSetValues
Date: July 17, 2009 7:44:09 AM CDT
To: petsc-users@mcs.anl.gov
Reply-To: petsc-users@mcs.anl.gov
The si
From: s.kramer@imperial.ac.uk
Subject: non-local values being dropped in MatSetValues
Date: July 17, 2009 7:44:09 AM CDT
To: petsc-users@mcs.anl.gov
Reply-To: petsc-users@mcs.anl.gov
The situation is the following: I'm calling MatSetValues with the flag ADD_VALUES and with matrix option MAT_IGNORE_ZERO_ENTRIES. I'm inserting multiple values at once, multiple columns and rows, so I provide a rank-2 matrix of values. As I'm calling this from fortran I'm also using MAT_COLUMN_ORIENTED. Now for provided rows that are not owned by the process, it jumps to mpiaij.c:394 (line numbers as in petsc-dev). On line 399, it checks for zero entries, but only checks the very first entry of the (non-owned) row. If however other entries of that same row are nonzero, the entire row is still dropped. Note that this is independent of row_oriented/column_oriented as line 396 does exactly the same.
If I don't set the option MAT_IGNORE_ZERO_ENTRIES the problem disappears. In that case however we would either have to preallocate substantially more nonzeros, or complicate the matrix assembly in our code by taking out the zero entries ourselves and call MatSetValues for each entry seperately.
Hg-commit: d28deff5e1b8c4e99cf2d919d60e6073ac99448f
show more ...
|
| #
2d0b4cb3
|
| 17-Jul-2009 |
Matthew Knepley <knepley@mcs.anl.gov> |
Merge
Hg-commit: df1ddc3c1dde73a02fe57cd49cb672f5616a135d
|
| #
c9568adc
|
| 17-Jul-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
commit after merge
Hg-commit: a72eb0b78c9d43fb9fa2d1a141d61473d55624cb
|
| #
85911e72
|
| 16-Jul-2009 |
Jed Brown <jed@59A2.org> |
Fix error in MatRelax_MPIAIJ for its=1 with nonzero initial guess.
This is critical for -pc_type ml (which uses -mg_levels_ksp_type
richardson by default) to work in parallel. It has been broken si
Fix error in MatRelax_MPIAIJ for its=1 with nonzero initial guess.
This is critical for -pc_type ml (which uses -mg_levels_ksp_type
richardson by default) to work in parallel. It has been broken since
changeset: c8eb7c4070f0
description:
only create work vector in MatRelax_MPIAIJ() if it is needed.
Hg-commit: 0e67bd686af5bb4b1b472881c2fe48e2026d7eac
show more ...
|
| #
be7314b0
|
| 15-Jul-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
removed unneeded inline spops.h file
removed unused macros from some inline files
moved axpy.h to include/private/petscaxpy.h
moved first Fortran kernels from central location to sys/utils/ftn-kernel
removed unneeded inline spops.h file
removed unused macros from some inline files
moved axpy.h to include/private/petscaxpy.h
moved first Fortran kernels from central location to sys/utils/ftn-kernels where it is used
Hg-commit: c9e7a1a8dc579cf824e1f087165c3df9932aefb9
show more ...
|
| #
1cec2a20
|
| 05-Jul-2009 |
Matthew Knepley <knepley@mcs.anl.gov> |
Merge
Hg-commit: 13c4ee66e6bd736c1a3c42fd00663fff1c898f67
|
| #
bd0c2dcb
|
| 05-Jul-2009 |
Barry Smith <bsmith@mcs.anl.gov> |
finished and tested MatRelax_MPIAIJ for AIJ matrices with inodes
Hg-commit: 44291efb1fe0ae95fc6834338633a87bc59425f4
|