| 1） 简单插值 do i=1,nxh-1 do j=1,nyh-1 f(i,j) = r(2*i,2*j) end do end do 测试结果 1 2.7965267341830442 0.99423048018935078 2 4.6530862847415096E-002 1.6542807028063357E-002 3 1.0567083868992405E-003 3.7568447820823846E-004 4 2.9529012639606583E-005 1.0498252728046983E-005 outer number of iteration = 4 normalized workload = 21.933593750000000 L2-norm = 1.8688712224537981E-006 Max-norm = 3.9565421108278542E-006 2） 5点格式插值 do i=1,nxh-1 do j=1,nyh-1 f(i,j) = 1.0d0/8.0d0*( 4.0d0*r(2*i,2*j) & + 1.0d0*(r(2*i+1,2*j)+r(2*i-1,2*j)+r(2*i,2*j+1)+r(2*i,2*j-1)) ) end do end do 1 2.7965267341830442 0.99423048018935078 2 4.6530862847415096E-002 1.6542807028063357E-002 3 1.0567083868992405E-003 3.7568447820823846E-004 4 2.9529012639606583E-005 1.0498252728046983E-005 outer number of iteration = 4 normalized workload = 21.933593750000000 L2-norm = 1.8688712224537981E-006 Max-norm = 3.9565421108278542E-006 3） 9点格式插值 1 2.7965267341830442 0.99423048018935078 2 0.10420304641789940 3.7046613433333320E-002 3 5.1698948876830048E-003 1.8380182151954728E-003 4 3.3861151941499344E-004 1.2038429292683420E-004 5 2.4816688111167853E-005 8.8229114479335366E-006 outer number of iteration = 5 normalized workload = 27.468750000000000 L2-norm = 7.8665572969769039E-007 Max-norm = 1.6521698827376419E-006 4）结论 简单的好像也可以哦 |
说点什么...