Branches

jed/mat-taij

View source
jed/mat-taij
  • Contributors
    1. Loading...
Author Commit Message Labels Comments Date
250+ commits behind master.
Jed Brown
IRK test: add coarse level
Jed Brown
WIP KSP ex61.c: checkpoint 3D Poisson with implicit Runge-Kutta
Jed Brown
WIP: TAIJ with nonzero T
Jed Brown
Mat: expose PetscLINPACKgefa and gedi --with-visibility
Jed Brown
Merge branch 'master' into jed/mat-taij Get build system updates and PetscMalloc1. Fix %G and malloc in MatRestoreRow_*TAIJ. * master: (2564 commits) Doc: Added some references do not print clang deprecation messages for SSL (on appble) PetscTextBelt() requires socket viewer test: remove unused test scripts nightlybuilds: add a test for c++/complex with 64-bit-indices [on osx with clang] malloc: PETSC_USE_MALLOC_COALESCED,…
Jed Brown
Mat TAIJ: fix memory management for SOR and optimize I (x) S + J (x) I
Jed Brown
MatMPITAIJ: add MatGetDiagonalBlock_MPITAIJ, to support block Jacobi/SOR
Jed Brown
KSP ex60: factor assembly into separate function
Jed Brown
KSP ex60: update output for "report time steps and Krylov iterations"
Jed Brown
Mat kernel: fix comment
Debojyoti Ghosh
MatTAIJ: Optimized the backward sweep of MatSOR_SeqTAIJ by saving the application of lower triangular during forward sweep.
Debojyoti Ghosh
MatTAIJ: Implemented MatSOR_SeqTAIJ
Jed Brown
KSP ex60: report time steps and Krylov iterations
Debojyoti Ghosh
MATTAIJ: Fixed two compilation warnings resulting from not calling CHKERRQ(ierr).
Jed Brown
KSP ex60: report final time, update output files
Jed Brown
KSP ex60: add centered-difference discretization of advection
Jed Brown
KSP ex60: make new arbitrary-order Gauss method the default
Jed Brown
KSP ex60: add support for Gauss methods of arbitrary order Includes test with 4-stage Gauss method
Jed Brown
KSP ex60: use PetscFunctionList to factor out 2-stage Gauss method
Jed Brown
KSP ex60: change variable name I to Identity When C99 complex.h is included, "I" is defined to be the imaginary number. Consequently, "I" cannot be used for any other purpose in PETSc.
Jed Brown
KSP ex60: fix __FUNCT__, PETSC_NULL -> NULL, and PetscOptionsBegin
Jed Brown
PCPBJacobi: add PCApplyTranspose (not optimized)
Jed Brown
PCPBJacobi: add PCApply_PBJacobi_N (not unrolled/optimized)
Jed Brown
PCPBJacobi: use VecGetArrayRead for unmodified vector
Jed Brown
PCPBJacobi: fix flop counting
Jed Brown
DTGaussQuadrature: fix quadrature weights
Debojyoti Ghosh
MATTAIJ: Setting block sizes for rectangular S,T also
Debojyoti Ghosh
MatConvert_Basic: Setting block sizes so as to handle rectangular blocks as well.
Debojyoti Ghosh
KSP ex60: Updated the test output for ex60
Debojyoti Ghosh
MATTAIJ: For square S,T, set the blocksize for the TAIJ matrix.
Comments 3
Debojyoti Ghosh
MATTAIJ: Updated the InvertBlockDiagonal for TAIJ to account for T, as well as special cases like S,T being NULL, T being identity.
Debojyoti Ghosh
MAT ex171: Added a test for MATTAIJ where T is specified as an identity matrix; and setting S such that most of its values are not 0!
Debojyoti Ghosh
MATTAIJ: Treating T = identity matrix as a special case since it reduces operations
Debojyoti Ghosh
MATTAIJ and MAT ex 171: Minor corrections in the comments
Debojyoti Ghosh
MATTAIJ: Added if conditions to handle cases where either or both of S and T may be NULL (so that arrays full of zeros are not allocated and operated on)
Debojyoti Ghosh
MAT ex171: Using the matrix "small" for ex171 testing (MATTAIJ functions); testing with nprocs 1,2,3,4 and different square and rectangular sizes for S,T
Debojyoti Ghosh
KSP ex60: Bugfix in the index of local uarr in exact solution function; and removing an unused variable
Debojyoti Ghosh
MATTAIJ: Corrected the comments to be consistent with the mathematical definition of the Kronecker product
Debojyoti Ghosh
MATTAIJ: Implemented MatGetRow (Seq and MPI) for MATTAIJ. Commented out the implementations of MatConvert; MATTAIJ will use MatConvert_Basic (which uses MatGetRow).
Comments 3
Debojyoti Ghosh
MATTAIJ: Passing Vec instead of Vec* in MatMult and MatMultAdd
Debojyoti Ghosh
MATTAIJ: Bugfix in MatMult for non-square I in I \otimes S
Debojyoti Ghosh
MAT_MPIAIJ: In MatMPIAIJGetSeqAIJ, check if pointers are NULL before assigning them. This will allow passing NULL pointers for information that is not needed by the calling function
Debojyoti Ghosh
KSP ex60: Reduced the diffusion coefficient of the test run (runex60) so that the initial sine wave does not decay to round-off errors and updated the test output file
Debojyoti Ghosh
KSP ex60: Updated creation of MATTAIJ to match the new interface, Using the extended TAIJ implementation to create the right-hand-side for solving the stages, and for step completion.
Debojyoti Ghosh
MAT ex171: updated ex171 (MATTAIJ tests) to form TAIJ with S and T; and try non-square S,T; Added tests for MATTAIJ where S or T can be NULL
Debojyoti Ghosh
MATTAIJ: Extended the TAIJ type to handle operations of the form [S \otimes I + T \otimes A] where S and T can be rectangular matrices; If NULL passed as S or T to MatCreateTAIJ, then set them to zero.
Jed Brown
MatTAIJ: fix MatInvertBlockDiagonal_SeqTAIJ, ksp ex60 works with pbjacobi
Jed Brown
KSP ex60: fix for Complex and report errors with %G
Debojyoti Ghosh
KSP ex60: Added a KSP example that demonstrates MATTAIJ by solving a constant-coefficient heat equation (1D) with an implicit Runge-Kutta method.
Debojyoti Ghosh
MAT ex171: Added an example 171 to test the MATTAIJ implementation.
Debojyoti Ghosh
MATTAIJ: Handling MatMult and MatMultAdd in TAIJ through a single function
Debojyoti Ghosh
MATTAIJ: Created and implemented with MatMult, MatMultAdd and InvertBlockDiagonal operations