Ivan Pribec ivan-pi
ivan-pi
/ boyd_bvp_example.f90
Created
March 18, 2019 00:59
Code example from book "Chebyshev and Fourier Spectral Methods" by John P. Boyd
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| program boyd_bvp | |
| dimension xi(20),aphi(20),g(20),h(20,20),ugraph(101) | |
| common/pbasis/phi(20),phix(20),phixx(20) | |
| common/thiele/th | |
| pi = 4.*atan(1.0) | |
| ! specify parameters |
ivan-pi
/ wort_cooling.py
Created
February 17, 2019 18:54
Some simple calculations of wort cooling with an internal coil
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| import numpy as np | |
| import matplotlib.pyplot as plt | |
| from scipy.integrate import solve_ivp | |
| # LMTD case | |
| def f(t,y,params): | |
| p2, p3, Tin = params |
ivan-pi
/ high_order_poisson_solver.f90
Created
February 11, 2019 16:09
Mehrstellenverfahren for the Poisson equation
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| ! See link below for details: | |
| ! https://scicomp.stackexchange.com/questions/6854/mehrstellenverfahren-for-poisson | |
| ! | |
| program poisson_higher_order | |
| implicit none | |
| integer, parameter :: wp = kind(1.0d0) | |
| real(wp), parameter :: pi = 4._wp*atan(0._wp) |
ivan-pi
/ test_banded.f90
Last active
November 6, 2023 10:30
Banded system of linear equations solved using "gbsv" driver from LAPACK (Intel MKL)
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| !LINK="${MKLROOT}/lib/intel64/libmkl_lapack95_lp64.a -Wl,--start-group ${MKLROOT}/lib/intel64/libmkl_intel_lp64.a ${MKLROOT}/lib/intel64/libmkl_sequential.a ${MKLROOT}/lib/intel64/libmkl_core.a -Wl,--end-group -lpthread -lm -ldl" | |
| !INC="-I${MKLROOT}/include/intel64/lp64/ -I${MKLROOT}/include/" | |
| !ifort -O3 -warn all -o test_banded $INC test_banded.f90 $LINK | |
| program test_banded | |
| use lapack95, only: gbtrf, gbtrs, gbsv | |
| use f95_precision, only: wp => dp | |
| implicit none | |
| integer, parameter :: kl = 2, ku = 1 |
ivan-pi
/ test_gttrf.f90
Last active
January 26, 2019 17:22
Tridiagonal linear system of equations example using gttrf from LAPACK (Intel MKL)
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| ! LINK="${MKLROOT}/lib/intel64/libmkl_lapack95_lp64.a -Wl,--start-group ${MKLROOT}/lib/intel64/libmkl_intel_lp64.a ${MKLROOT}/lib/intel64/libmkl_sequential.a ${MKLROOT}/lib/intel64/libmkl_core.a -Wl,--end-group -lpthread -lm -ldl" | |
| ! INC="-I${MKLROOT}/include/intel64/lp64/ -I${MKLROOT}/include/" | |
| ! ifort -warn all -o test_gttrf $INC test_gttrf.f90 $LINK | |
| ! ./test_gttrf > test_gttrf.out | |
| ! gnuplot | |
| ! gnuplot> Th = 1 | |
| ! gnuplot> plot "test_gttrf.out" u 1:2 w p pt 7, cosh((1-x)*Th)/cosh(Th) | |
| module rd_setup |
ivan-pi
/ mkl_least_squares_example.f90
Last active
December 28, 2018 23:14
Sample least squares problem using LAPACK routines from Intel MKL
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| ! # Check out the link below for specific compile instructions for your system: | |
| ! # https://software.intel.com/en-us/articles/intel-mkl-link-line-advisor/ | |
| ! | |
| ! To compile use: | |
| ! LINK="${MKLROOT}/lib/intel64/libmkl_lapack95_lp64.a -Wl,--start-group ${MKLROOT}/lib/intel64/libmkl_intel_lp64.a ${MKLROOT}/lib/intel64/libmkl_sequential.a ${MKLROOT}/lib/intel64/libmkl_core.a -Wl,--end-group -lpthread -lm -ldl" | |
| ! INC="-I${MKLROOT}/include/intel64/lp64/ -I${MKLROOT}/include/" | |
| ! ifort -warn all -o example3 $INC example3.f90 $LINK | |
| module m_least_squares |
ivan-pi
/ kalman.f90
Last active
October 24, 2018 08:19
— forked from mrocklin/kalman.f90
Kalman filter in BLAS/LAPACK Fortran; see [blogpost](http://matthewrocklin.com/blog/work/2012/11/24/Kalman-Filter) for more information
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| subroutine f(mu, Sigma, H, INFO, R, Sigmavar_2, data, muvar_2, k, n) | |
| implicit none | |
| integer, intent(in) :: k | |
| integer, intent(in) :: n | |
| real*8, intent(in) :: Sigma(n, n) ! Sigma | |
| real*8, intent(in) :: H(k, n) ! H | |
| real*8, intent(in) :: mu(n) ! mu | |
| real*8, intent(in) :: R(k, k) ! R, H*Sigma*H' + R | |
| real*8, intent(in) :: data(k) ! (H*Sigma*H' + R)^-1*((-1)*data + H*mu), data, (-1)* data + H*mu |
NewerOlder