monolish
MONOlithic LInear equation Solvers for Highly-parallel architecture
|
Go to the documentation of this file.
10 #include <initializer_list>
13 #define MONOLISH_SOLVER_SUCCESS 0
14 #define MONOLISH_SOLVER_SIZE_ERROR -1
15 #define MONOLISH_SOLVER_MAXITER -2
16 #define MONOLISH_SOLVER_BREAKDOWN -3
17 #define MONOLISH_SOLVER_RESIDUAL_NAN -4
18 #define MONOLISH_SOLVER_NOT_IMPL -10
273 template <
typename T>
281 template <
typename T,
typename U>
296 template <
typename T>
311 template <
typename T>
325 template <
typename T>
339 template <
typename T>
353 template <
typename T>
360 template <
typename T,
typename... types>
362 const types &... args) {
376 template <
typename T,
typename U>
378 return x.size() == y.size();
391 template <
typename T>
405 template <
typename T>
419 template <
typename T>
433 template <
typename T>
440 template <
typename T,
typename U,
typename... types>
442 const types &... args) {
453 template <
typename T,
typename U>
455 return arg1.get_device_mem_stat() == arg2.get_device_mem_stat();
465 template <
typename T,
typename U,
typename... types>
467 const types &... args) {
486 template <
typename T>
488 const T diag_val,
const T val);
501 template <
typename T>
535 template <
typename T>
548 template <
typename T>
563 template <
typename T>
575 template <
typename T>
587 template <
typename T>
600 template <
typename T>
616 template <
typename T>
632 template <
typename T>
634 T a0, T a1, T a2, T b0,
642 template <
typename T>
auto send(T &x) { x.send(); }
647 template <
typename T,
typename... Types>
auto send(T &x, Types &... args) {
656 template <
typename T>
auto recv(T &x) { x.recv(); }
661 template <
typename T,
typename... Types>
auto recv(T &x, Types &... args) {
676 template <
typename T,
typename... Types>
matrix::COO< T > random_structure_matrix(const int M, const int N, const int nnzrow, const T val)
create random structure matrix (column number is decided by random)
bool is_same_size(const T &x, const U &y)
compare size of vector or 1Dview (same as is_same_structure())
bool build_with_avx()
get build option (true: with avx, false: without avx)
matrix::COO< T > laplacian_matrix_2D_5p(const int M, const int N)
create two dimensional Laplacian matrix using the five point central difference scheme
Linear Operator imitating Matrix.
void min(const matrix::CRS< double > &A, const matrix::CRS< double > &B, matrix::CRS< double > &C)
Create a new CRS matrix with smallest elements of two matrices (C[0:nnz] = min(A[0:nnz],...
bool build_with_gpu()
get build option (true: enable gpu, false: disable gpu)
matrix::COO< T > laplacian_matrix_1D(const int &M)
create 1D Laplacian matrix
bool build_with_cblas()
get build option (true: with cblas, false: without cblas (=with intel mkl))
bool is_same_structure(const T A, const U B)
compare matrix structure
void max(const matrix::CRS< double > &A, const matrix::CRS< double > &B, matrix::CRS< double > &C)
Create a new CRS matrix with greatest elements of two matrices (C[0:nnz] = max(A[0:nnz],...
double get_residual_l2(const matrix::Dense< double > &A, const vector< double > &x, const vector< double > &y)
get nrm |b-Ax|_2
T tridiagonal_toeplitz_matrix_eigenvalue(const int &M, int N, T a, T b)
Nth smallest eigenvalue of MxM tridiagonal Toeplitz matrix.
auto device_free(T &x)
free data of GPU
T frank_matrix_eigenvalue(const int &M, const int &N)
Nth eigenvalue from the bottom of MxM Frank matrix.
auto size() const
get vector size
bool build_with_lapack()
get build option (true: with lapack, false: without lapack (=with intel mkl))
bool is_same_device_mem_stat(const T &arg1, const U &arg2)
compare same device memory status
matrix::COO< T > frank_matrix(const int &M)
create Frank matrix
matrix::COO< T > tridiagonal_toeplitz_matrix(const int &M, T a, T b)
create tridiagonal Toeplitz matrix
T laplacian_matrix_1D_eigenvalue(const int &M, int N)
Nth smallest eigenvalue of 1D Laplacian matrix.
T toeplitz_plus_hankel_matrix_eigenvalue(const int &M, int N, T a0, T a1, T a2, T b0, T b1, T b2)
Nth smallest eigenvalue of GEVP Ax=lBx of Toeplitz-plus-Hankel matrixes A, B.
bool set_default_device(size_t device_num)
set default device number
matrix::COO< T > band_matrix(const int M, const int N, const int W, const T diag_val, const T val)
create band matrix
Coodinate (COO) format Matrix (need to sort)
void set_log_filename(const std::string filename)
Specifying the log finename.
auto recv(T &x)
recv. and free data from GPU
bool solver_check(const int err)
check error
void random_vector(vector< T > &vec, const T min, const T max)
create random vector
bool build_with_avx2()
get build option (true: with avx2, false: without avx2)
void set_log_level(const size_t Level)
Logger utils ///////////////////////////////.
int get_default_device()
get default device number
matrix::COO< T > toeplitz_plus_hankel_matrix(const int &M, T a0, T a1, T a2)
create Toeplitz-plus-Hankel matrix
bool build_with_mpi()
get build option (true: enable MPI, false: disable MPI)
int get_num_devices()
get the number of devices
auto send(T &x)
send data to GPU
matrix::COO< T > eye(const int M)
create band matrix
Compressed Row Storage (CRS) format Matrix.
bool build_with_avx512()
get build option (true: with avx512, false: without avx512)
bool build_with_mkl()
get build option (true: with intel mkl, false: without intel mkl)