双直線補間について
http://vldb.gsi.go.jp/sokuchi/patchjgd/helpweb/words.html
→バイリニア補間の項を参照
プログラムとmakefile
コンパイル例
$ make
ifort -c -CB -traceback -fpe0 main.f90
ifort -c -CB -traceback -fpe0 bilnr.f90
ifort -o bl main.o bilnr.o
real 0m0.111s
user 0m0.071s
sys 0m0.040s
実行例
[Fri Sep 16 10:49:01 JST 2011]
[am@aofd30 processor=x86_64]
[~/teaching/bilinear]
$ ./bl
Program name: ./bl
For given xa(i),ya(i),za(i),x & y, z is estimated.
(xa(4),ya(4),za(4)) (xa(3),ya(3),za(3))
+------------------+
| |
| x |
| (x,y,z) |
| |
| |
| |
| |
| |
+------------------+
(xa(1),ya(1),za(1)) (xa(2),ya(2),za(2))
INPUT
i, xa(i),ya(i),za(i) = 1 0.0 0.0 1.0
i, xa(i),ya(i),za(i) = 2 1.0 0.0 2.0
i, xa(i),ya(i),za(i) = 3 1.0 1.0 2.0
i, xa(i),ya(i),za(i) = 4 0.0 1.0 1.0
OUTPUT
x, y, z = 0.9 0.5 1.9
./bl terminated nomally.