ic0(3)
NAME
ic0 - incomplete Choleski factorization
DESCRIPTION
- The class implements the icomplete Choleski factorization IC0(A)
when $A$ is a square symmetric matrix stored in CSR format. - csr<double> a = ...;
ic0<double> fact(a);
IMPLEMENTATION
- template <class T>
class basic_ic0 : public csr<T> {
public: - typedef typename csr<T>::size_type size_type;
typedef T element_type; - // constructors:
- basic_ic0 ();
explicit basic_ic0 (const csr<T>& a); - // solver:
- void inplace_solve (vec<T>& x) const;
void solve (const vec<T>& b, vec<T>& x) const;
vec<T> solve (const vec<T>& b) const; - // accessors:
- size_type nrow () const { return csr<T>::nrow(); }
size_type ncol () const { return csr<T>::ncol(); }
size_type nnz () const { return csr<T>::nnz(); } - };
template <class T>
basic_ic0<T> ic0 (const csr<T>& a);