Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Transfer linear algebra used in DDC splines to kokkos kernel #393

Open
jbigot opened this issue Apr 8, 2024 · 5 comments
Open

Transfer linear algebra used in DDC splines to kokkos kernel #393

jbigot opened this issue Apr 8, 2024 · 5 comments
Assignees

Comments

@jbigot
Copy link
Member

jbigot commented Apr 8, 2024

No description provided.

@cedricchevalier19
Copy link
Member

Meeting 5/13.

@blegouix
Copy link
Collaborator

blegouix commented May 6, 2024

Should I join ? Give me the time of the meeting if so

@yasahi-hpc
Copy link
Member

@tpadioleau
Copy link
Member

tpadioleau commented May 27, 2024

General dense

LAPACKE_dgetrs
LAPACKE_dgetrf

General banded

LAPACKE_dgbtrs
LAPACKE_dgbtrf

PDS banded

LAPACKE_dpbtrs
LAPACKE_dpbtrf

PDS tridiagonal

LAPACKE_dpttrs
LAPACKE_dpttrf

@yasahi-hpc
Copy link
Member

yasahi-hpc commented May 28, 2024

I have made a small investigation.

Function dgetrs dgbtrs dpbtrs dpttrs
Brief general N-by-N matrix solver. A is factorized by dgetrf general band matrix solver. A is factorized by dgbtrf symmetric positive definite band matrix solver. A is factorized by dpbtrf tridiagonal system solver. A is factorized by dpttrf
Dependency dlaswp, dtrsm dgemv, dger, dswap, dtbsv dtbsv dptts2
In Kokkos kernels - - - -

These functions depend on following subroutines.

Subroutines dlaswp dtrsm dgemv dger dswap dtbsv dptts2 dscal
Brief performs series of row interchanges on the matrix A solves one of the matrix equations op( A )*X = alpha*B, or X*op( A ) = alpha*B performs one of the matrix-vector operations y := alpha*A*x + beta*y, or y := alpha*A**T*x + beta*y performs the rank 1 operation A := alpha*x*y**T + A interchanges two vectors solves one of the systems of equations A*x = b, or A**T*x = b solves a tridiagonal system of the form A * X = B using the L*D*L' factorization of A computed by dpttrf scales a vector by a constant
Serial - x x - - x - -
Team - x x - - - - -
Team Vector - x - - - - - -
Note blas implementation available In review Depends on dscal blas implementation available

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

5 participants