| 查看: 1040 | 回复: 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 请教，这个的功能是怎样的，谢谢！ |
回复此楼» 猜你喜欢
析晶
已经有5人回复
-
国基评审
已经有9人回复
-
国自然面上和省基金B类撒花
已经有22人回复
-
2026-博士申请
已经有4人回复
-
26级硕士毕业生求博导收留
已经有4人回复
-
湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息
已经有6人回复
-
考研调剂
已经有3人回复
-
急招9月入学博士，要有4级、最晚7月硕士毕业。精密电机驱控课题；学位材料
已经有5人回复
-
又一批高校组建人工智能学院 师资行吗 不是骗人吗
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人