http://en.wikipedia.org/wiki/Machine_epsilon http://en.wikipedia.org/wiki/Machine_epsilon#Approximation_using_Fortran Rounding is a procedure for choosing the representation of a real number in a floating point number system. For a number system and a rounding procedure, machine epsilon is the maximum relative error of the chosen rounding procedure. Why doesn’t my application always give the same answer? http://software.intel.com/sites/default/files/article/164389/fp-consistency-122712_1.pdf

https://github.com/orange224/su_scripts/blob/master/eps2png.sh #/bin/bash # http://www.imagemagick.org/script/command-line-options.php#density # -depth: This option specifies the image resolution to store while encoding a raster image # or the canvas resolution while rendering (reading) vector formats such as Postscript, # PDF, WMF, and SVG into a raster image. # The default resolution is 72 dots per inch, which is equivalent to one […]

http://netlib.org/amos/ https://github.com/certik/fortran-utils/blob/master/src/amos.f90 http://dl.acm.org/citation.cfm?id=214331 Amos, D.E., “A Portable Package for Bessel Functions of a Complex Argument and Nonnegative Order,” Trans. Math. Software, 1986 fortran90코드의 계산값과 matlab결과값을 비교한 그림입니다. This figure shows the self-written Fortran 90 program’s results of Bessel function of third kind(Hankel funcion). The result has the difference with Matlab’s besslh function of ~ 10^-14 […]

생각보다 귀찮아서 제가 필요한 만큼만 딱 subroutines으로 만들어서 배포합니다. 코드는 이메일로 사용목적(ex. 학생이다 등등)을 간단하게 알려주시면 제가 판단해서 보내드립니다. ps) real(kind=8) 값이 정상적으로 입력되려면 compile할때 intel fortran complier의 경우에는 -r8 옵션을 주어야 합니다. 아래 두 예제에서 real(kind=8)로 선언된 dx값을 비교해보시면 무슨 이야기인지 감이 오실 듯. :D 예제 1. -r8 옵션이 없는 경우 [worker@local]$ ifort -o […]

http://people.sc.fsu.edu/~jburkardt/f_src/special_functions/special_functions.html

Collino, F. & Tsogka, C., 2001. Application of the perfectly matched absorbing layer model to the linear elastodynamic problem in anisotropic heterogeneous media. GEOPHYSICS, 66(1), pp.294–307. Available at: http://library.seg.org/doi/abs/10.1190/1.1444908 equation (20) and (21) for heterogeneous isotropic elastic medium. equation (23) and (xx) for homogeneous anisotropic elastic medium. Collino, F. & Monk, P., 1998. Optimizing the […]

acoustics 공부를 좀더 해봐야겠다.

Follow

Get every new post delivered to your Inbox.