libcKrylov

A generic implementation of linear solvers in Krylov space

Author: Pengliang Yang, Harbin Institute of Technology

E-mail: ypl.2100@gmail.com

Linear solvers in Krylov space

• Conjugate gradient (CG)

• Preconditioned conjugate gradient (PCG)

• BiCGStab (Bi-conjugate gradient with stabalization)

• BiCGStab with right preconditioning

• GMRES (generalized minimum residual)

• GMRES with right preconditioning

• CGNR

• CGNE

How to install and run

1. First, compile source code for Krylov solvers library

   cd src;

   make

   After this step, the implementation of these linear solvers will be compiled to form a library in /lib, named libcKrylov.a. To remove the old library and compile for a new one, 'make clean' and 'make'.

2. Then, test with the application examples!

Run the Toeplitz matrix inversion using CGNR and CGNE:

cd demo_toeplitz;

make;

./main

Run the demo for frequency domain wave equation using BiCGStab and GMRES:

cd demo_wave;
	
make;

./main

Run the demo for Poisson equation using CG and PCG:

cd demo_poisson;

make;

./main

Note that you are free to modify the parameters in main.c to test different algorithm options and the performance.