| 查看: 374 | 回复: 2 | ||
[求助]
小程序求助！！！
|
|
紧急求助啊，下面的程序出现log domain error. 请大家帮忙怎么解决啊！急啊！非常感谢！！！ Program main integer :: ix,iy real(8),dimension(1,1,5,5) :: xp,x real(8),dimension(1:2,5,5) :: phi xp(1,1,1,1)=0.2 x(1,1,1,1)=0.2 phi(1,1,1)=1 phi(2,1,1)=0 do ix=1,5 do iy=1,5 xp(1,1,ix,iy)=xp(1,1,ix,iy)+0.01 x(1,ix,iy,1)= x(1,ix,iy,1)+ 0.05 phi(1,ix,iy)=phi(1,ix,iy)-0.001 phi(1,ix,iy)=phi(2,ix,iy)+0.001 call calculation (xp(1,1,ix,iy),x(1,ix,iy,1),phi(1,ix,iy),phi(2,ix,iy)) xp(1,2,ix,iy)=(x(1,ix,iy,1)-phi(1,ix,iy)*xp(1,1,ix,iy))/phi(2,ix,iy) end do end do end program main subroutine calculation (x,xc,phi1,phi2) integer :: k real(8) :: tol,x,x0,xc,f,df,phi1,phi2 f(x)=1.06366481d10+7.94800172d9+1.72992961d9*log(abs(x))-1.72992961d9*log(abs(1-x))-1.182988286d10*(1-2*x)+2.75714286d9*(6*x-6*x**2-1)+3.083731429d9*(24*x**2-10*x-16*x**3+1)-1d11*(xc-phi1*x)/phi2 + 9.229804981d10 df(x)=2*1.182988286d10+1.72992961d9/x/(1-x)+2.75714286d9*(6-12*x)+3.083731429d9*(48*x-10-48*x**2)+1d11*phi1/phi2 tol=1e-6 x0=0.2d0 k=0 10 k=k+1 x=x0-f(x0)/df(x0) if(abs(x-x0).lt.tol)goto 20 x0= x goto 10 20 write(*,*) x end subroutine calculation |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
求助多输出支持向量机的matlab程序
已经有6人回复
- perl程序求助，请大神们赐教~~~
已经有9人回复
- c语言求助，一个小程序详细讲解
已经有13人回复
- 求处理vasp能带的小程序
已经有14人回复
- 求助高手帮忙编一个Fortran的小程序，重金悬赏......
已经有15人回复
- 求助：dushin程序输入文件无法读取！！！1
已经有10人回复
- field calculator VBA 条件判断小程序求助
已经有3人回复
- 急切求助_项目进展效果图制作
已经有6人回复
- 求助：后缀名为.java和.class的程序需要下载什么软件进行编译？
已经有16人回复
- 【求助】C# 中如何声明结构体数组并赋值
已经有13人回复
- 【求助】请问如何解决insufficient virtual memory
已经有10人回复
- 【求助】lammps安装过程中，在安装fftw的过程中遇到问题
已经有8人回复
- 【求助】小程序
已经有6人回复
- 【求助】VS2005数值计算程序调试【急求！！！！加金币】
已经有9人回复
- 【求助】购买程序降温的烘箱求助！！
已经有11人回复
- 【求助】请问一下fortran的可视化编程
已经有7人回复
- 【求助】VC++.NET编程，MFC基于对话框应用程序编写
已经有18人回复
- 【求助】Fortran程序，如何让exe文件可以单独运行？
已经有6人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
- 【求助】求助，如何使linux能运行windows dos 的程序！
已经有5人回复
- 【求助】帮忙编个用于VASP小程序
已经有22人回复
snoopyzhao
至尊木虫 (职业作家)
- 程序强帖: 16
- 应助: 157 (高中生)
- 贵宾: 0.02
- 金币: 18844.7
- 红花: 29
- 帖子: 3803
- 在线: 1422.4小时
- 虫号: 183750
- 注册: 2006-02-13
- 专业: 污染生态化学
2楼2012-11-28 09:08:24
|
谢谢你的回复。我又重新写了一下，但他说我的f 和df size 不符。你可以帮我看看吗？非常感谢！ Program main integer :: ix,iy real(8),dimension(1,1,5,5) :: xp,x real(8),dimension(1:2,5,5) :: phi do ix=1,5 do iy=1,5 x(1,ix,iy,1)= 0.3 phi(1,ix,iy)=0.88 phi(1,ix,iy)=0.12 end do end do call calculation (xp(1,1,ix,iy),x(1,ix,iy,1),phi(1,ix,iy),phi(2,ix,iy)) !xp(1,2,ix,iy)=(x(1,ix,iy,1)-phi(1,ix,iy)*xp(1,1,ix,iy))/phi(2,ix,iy) end program main subroutine calculation (x,xc,phi1,phi2) integer :: k,i,j real(8),dimension(1,1,5,5) :: x,xc,x0 real(8),dimension(1,5,5) :: phi1,phi2 real(8),dimension(1,1,5,5) :: f,df real(8):: tol do i=1,5 do j=1,5 f(1,1,i,j)=1.06366481d10+7.94800172d9+1.72992961d9*log(abs(x(1,1,i,j)))-1.72992961d9*log(abs(1-x(1,1,i,j)))& -1.182988286d10*(1-2*x(1,1,i,j))+2.75714286d9*(6*x(1,1,i,j)-6*x(1,1,i,j)**2-1)+3.083731429d9*(24*x(1,1,i,j)**2& -10*x(1,1,i,j)-16*x(1,1,i,j)**3+1)-1d11*(xc(1,i,j,1)-phi1*x(1,1,i,j))/phi2 + 9.229804981d10 df(1,1,i,j)=2*1.182988286d10+1.72992961d9/x(1,1,i,j)/(1-x(1,1,i,j))+2.75714286d9*(6-12*x(1,1,i,j))+& 3.083731429d9*(48*x(1,1,i,j)-10-48*x(1,1,i,j)**2)+1d11*phi1/phi2 tol=1e-6 x0(1,1,i,j)=0.2d0 k=0 10 k=k+1 x(1,1,i,j)=x0(1,1,i,j)-f(1,1,i,j)/df(1,1,i,j) if(abs(x(1,1,i,j)-x0(1,1,i,j)).lt.tol)goto 20 x0(1,1,i,j)= x(1,1,i,j) goto 10 20 write(*,*) x(1,1,i,j) end do end do end subroutine calculation |
3楼2012-11-28 21:41:09
