Sparse Matrix
Exercise
I implemented
- BiCG and BiCGstab
- CRS method(CSR) for understaning the mechanism.
Future plan
- Multi Color ordering
- Cuthill_Mckee method
Fortran magic
You can use Fortran code in Jupyter notebook if you install fortran-magic through python-pip.
pip3 install -U fortran-magic
In jupyter notebook,
%load_ext fortranmagic # activating fortran-magic
In[1] %%fortran
subroutine solve(x, y, z)
real, intent(in) :: x(:), y(:)
real, intent(out) :: z(size(x))
z(:) = sin(x(:) + y(:))
end subroutine solve
In[2] solve([1, 2, 3], [4, 5, 6])
array([-0.9589243, 0.6569866, 0.4121185], dtype=float32)