14 #ifndef DECOMPOSE_LU_CMPLX_INCLUDED
15 #define DECOMPOSE_LU_CMPLX_INCLUDED
33 void solve(
double *vec);
48 std::valarray<double>
m_lu;
50 inline size_t re(
int i,
int j)
52 return N2 * i + 2 * j;
55 inline size_t im(
int i,
int j)
57 return N2 * i + 2 * j + 1;
60 inline size_t re(
int i)
65 inline size_t im(
int i)
Decompose_LU_Cmplx(size_t N)
void mult_inverse(double *mat)
std::valarray< double > m_lu
void get_inverse(double *mat_inv)
void set_matrix(const double *mat)