| 查看: 1058 | 回复: 0 | |||
[求助]
求助：ludcmp&lubksb
|
|
subroutine ludcmp(coeff, numvar, np, indx, d) integer :: numvar, np, indx(numvar) doubleprecision :: d, coeff(np,np) integer, parameter :: nmax=500 doubleprecision, parameter :: tiny=1.0d-20 integer :: i, imax, j, k doubleprecision :: aamax, dum, sum, vv(nmax) d=1 do i = 1, numvar aamax = 0.d0 do j = 1, numvar if (abs(coeff(i,j)).gt.aamax) aamax = abs(coeff(i,j)) end do !if (aamax .eq. 0.d0) pause 'singular matrix in ludcmp' vv(i) = 1./aamax end do do j = 1, numvar do i= 1, j-1 sum = coeff(i,j) do k=1,i-1 sum = sum-coeff(i,k)*coeff(k,j) end do coeff(i,j) = sum end do aamax = 0.d0 do i = j,numvar sum = coeff(i,j) do k=1, j-1 sum = sum - coeff(i,k) * coeff(k,j) end do coeff(i,j) = sum dum = vv(i) * abs(sum) if (dum .ge. aamax) then imax = i aamax = dum end if end do if (j .ne. imax) then do k = 1, numvar dum = coeff(imax,k) coeff(imax,k) = coeff(j,k) coeff(j,k) = dum end do d = -d vv(imax) = vv(j) end if indx(j) = imax if (coeff(j,j) .eq. 0.d0) coeff(j,j) = tiny if (j .ne. numvar) then dum = 1./coeff(j,j) do i = j+1, numvar coeff(i,j) = coeff(i,j) * dum end do end if end do return end subroutine ludcmp 请教，这个的功能是怎样的，谢谢！ |
回复此楼» 猜你喜欢
综述论文作为代表作会不会影响评审专家的印象分？
已经有12人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有29人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有9人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
-
为什么网上很多人说本周 12号出结果
已经有7人回复
-
读博的好处
已经有3人回复
找到一些相关的精华帖子，希望有用哦~
【求助】除去2-溴丙酸乙酯
已经有3人回复
- 求助啊，SCI修稿后一月有余，询问编辑，给了一封回信，大家帮忙看一下是什么意思啊
已经有24人回复
- 求助：气相色谱测定醇醚平行性差如何解决
已经有8人回复
- 【求助】请教合成CdTe@SiO2
已经有6人回复
- 【求助/交流】NCBI给了我登录号但是用号查不到，这样算不算登录成功呢？
已经有8人回复
- 求助：Materials Chemistry and Physics的word 模板
已经有5人回复
- 【求助】分离黄酮和酚酸
已经有12人回复
- 【求助】羧酸的胺盐转化为酰胺
已经有5人回复
- 【求助】XRD几个小问题求助
已经有3人回复
- 【求助】请求帮助：测得一个晶体数据有这么多错误！怎么办？
已经有4人回复
- 【求助/交流】克隆Amp/Kan选择以及蓝白筛选
已经有10人回复
- 【求助】关于顺磁性物质和抗磁性物质
已经有12人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
