- c - 在位数组中找到第一个零
- linux - Unix 显示有关匹配两种模式之一的文件的信息
- 正则表达式替换多个文件
- linux - 隐藏来自 xtrace 的命令
我已经写下Fortran代码来计算距离然后排序,但是调用可执行命令有一些问题。
这是代码
program sort
implicit none
character CN*8,O*7
integer j,iconf,nconf
integer i,m
integer n,nmax,num
parameter (n=5)
double precision xbox,rq
parameter (nmax=3091,nconf=1)
double precision atom(nmax),id(nmax),ox(nmax),oy(nmax),oz(nmax)
double precision xij,yij,zij,rij,t
double precision r(n,n)
open(unit=1,status='unknown',file='a.gro')
do iconf= 1,nconf
read(1,*)
read(1,*)
do i=1,n
read(1,'(A8,A7,1i5,3f8.3)')CN,O,num,ox(i),oy(i),oz(i)
enddo
read(1,*)xbox
open(unit=3,file='dist.txt')
do i=1,n
do j=1,n
if(i .ne. j) then
xij=ox(i)-ox(j)
yij=oy(i)-oy(j)
zij=oz(i)-oz(j)
xij=xij - nint(xij/xbox)*xbox
yij=yij - nint(yij/xbox)*xbox
zij=zij - nint(zij/xbox)*xbox
r(i,j)=dsqrt(xij**2 + yij**2 + zij**2)
write(3,'(i3,2x,i3,4x,f17.15)') i,j, r(i,j)
call execute_command_line(sort -t, -k1 -g r(i,j))
write(*,*)
endif
enddo
enddo
enddo
END program
输入文件是a.gro
Generated by trjconv : 360 water t= 1000.00000
216
1water OW1 1 0.764 0.617 0.582
2water OW1 2 0.865 1.469 1.696
3water OW1 3 0.423 1.400 1.324
4water OW1 4 0.381 1.464 0.392
5water OW1 5 1.279 0.872 0.131
1.87759 1.87759 1.87759
outfile 文件 3,dist.txt
1 2 1.148553302245917
1 3 1.131341681367747
1 4 0.948787647474397
1 5 0.730514202462895
2 1 1.148553302245917
2 3 0.581815262776768
2 4 0.750524142249935
2 5 0.790896648178509
3 1 1.131341681367747
3 2 0.581815262776768
3 4 0.935138492417032
3 5 1.216627908647504
4 1 0.948787647474397
4 2 0.750524142249935
4 3 0.935138492417032
4 5 1.106792211754311
5 1 0.730514202462895
5 2 0.790896648178509
5 3 1.216627908647504
5 4 1.106792211754311
所以,我想对 r(i,j)
进行排序,保持 i 相同 j 不同。但是调用行在 fortran 代码中不起作用。
即将发生的错误
tetra.f(48): error #6413: This global name is invalid in this context. [SORT]
call execute_command_line(sort -t, -k1 -g r(i,j))
----------------------------------^
tetra.f(48): error #6404: This name does not have a type, and must have an explicit type. [K1]
call execute_command_line(sort -t, -k1 -g r(i,j))
--------------------------------------------^
tetra.f(48): error #6404: This name does not have a type, and must have an explicit type. [GR]
call execute_command_line(sort -t, -k1 -g r(i,j))
------------------------------------------------^
tetra.f(48): error #6362: The data types of the argument(s) are invalid. [EXECUTE_COMMAND_LINE]
call execute_command_line(sort -t, -k1 -g r(i,j))
---------------------------------------^
tetra.f(48): error #6362: The data types of the argument(s) are invalid. [EXECUTE_COMMAND_LINE]
call execute_command_line(sort -t, -k1 -g r(i,j))
-----------------------------------------------^
compilation aborted for tetra.f (code 1)
请告诉我如何在 Fortran 代码中执行 shell 命令。
最佳答案
我知道您想要做的是为 i 的每个可能值(我在下文中称为 irow)排序 r(i,j) 的四个值。如果是这种情况,并且您需要 Fortran 答案(而不是 Linux 答案),那么以下内容应该有效。排序是计算机科学中的一个核心问题,并且有很多算法。我选择冒泡排序是因为它相当容易理解,并且有很好的文档记录,例如 https://en.wikipedia.org/wiki/Bubble_sort .这绝对不是很快。请注意,排序后的结果值最小的在前,顺序为升序。
Module BSort
use, intrinsic :: iso_c_binding
contains
Subroutine f_Bubble_Sort(a,order)
implicit none
integer (c_int),intent(out) :: order(:)
integer (c_int) :: kx,ky,length,otemp
real (c_float),intent(in) :: a(:)
real (c_float),allocatable :: locala(:)
real (c_float) :: temp
logical(c_bool) :: swapped
length=size(a)
allocate(locala(length));order=(/(kx,kx=1,length)/);locala=a
do kx = length-1, 1, -1 !from top -> down
swapped = .false.
do ky = 1, kx
if (locala(ky) > locala(ky+1)) then
temp = locala(ky)
locala(ky) = locala(ky+1)
locala(ky+1) = temp
otemp=order(ky)
order(ky)=order(ky+1)
order(ky+1)=otemp
swapped = .true.
end if
end do
if (.not. swapped) exit
end do
End Subroutine f_Bubble_Sort
End Module BSort
Program Q52001740
use, intrinsic :: iso_c_binding
use BSort
implicit none
real(kind=c_float) :: r(5,4)
integer(c_int) :: out(4),irow,jcol,column(5,4),order(5,4)
open(unit=3,file='Q52001740.dist.txt')
! note no diagonal values in r
do irow=1,5
do jcol=1,4
read(3,'(5x,i3,4x,f17.15)') column(irow,jcol), r(irow,jcol)
end do
call f_Bubble_Sort(r(irow,:),out)
order(irow,:)=out
write(*,fmt='("irow= ",i2)')irow
write(*,fmt='(i2,4f8.5,4i3)')irow,r(irow,:),out
write(*,fmt='(a10,i2,4f8.5)')"sorted r:",irow,r(irow,out)
write(*,fmt='(a16,4i4)')"sorted columns:",column(irow,out)
end do
End Program Q52001740
输出是:
irow= 1
1 1.14855 1.13134 0.94879 0.73051 4 3 2 1
sorted r: 1 0.73051 0.94879 1.13134 1.14855
sorted columns: 5 4 3 2
irow= 2
2 1.14855 0.58182 0.75052 0.79090 2 3 4 1
sorted r: 2 0.58182 0.75052 0.79090 1.14855
sorted columns: 3 4 5 1
irow= 3
3 1.13134 0.58182 0.93514 1.21663 2 3 1 4
sorted r: 3 0.58182 0.93514 1.13134 1.21663
sorted columns: 2 4 1 5
irow= 4
4 0.94879 0.75052 0.93514 1.10679 2 3 1 4
sorted r: 4 0.75052 0.93514 0.94879 1.10679
sorted columns: 2 3 1 5
irow= 5
5 0.73051 0.79090 1.21663 1.10679 1 2 4 3
sorted r: 5 0.73051 0.79090 1.10679 1.21663
sorted columns: 1 2 4 3
关于linux - 在 Fortran 代码中执行 execute_command_line() 时出错,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/52001740/
关闭。这个问题不符合Stack Overflow guidelines .它目前不接受答案。 要求我们推荐或查找工具、库或最喜欢的场外资源的问题对于 Stack Overflow 来说是偏离主题的,
Linux 管道可以缓冲多少数据?这是可配置的吗? 如果管道的两端在同一个进程中,但线程不同,这会有什么不同吗? 请注意:这个“同一个进程,两个线程”的问题是理论上的边栏,真正的问题是关于缓冲的。 最
我找到了here [最后一页] 一种有趣的通过 Linux 启动 Linux 的方法。不幸的是,它只是被提及,我在网上找不到任何有用的链接。那么有人听说过一种避免引导加载程序而使用 Linux 的方法
很难说出这里要问什么。这个问题模棱两可、含糊不清、不完整、过于宽泛或夸夸其谈,无法以目前的形式得到合理的回答。如需帮助澄清此问题以便重新打开,visit the help center . 关闭 1
我试图了解 ld-linux.so 如何在 Linux 上解析对版本化符号的引用。我有以下文件: 测试.c: void f(); int main() { f(); } a.c 和 b.c:
与 RetroPie 的工作原理类似,我可以使用 Linux 应用程序作为我的桌面环境吗?我实际上并不需要像实际桌面和安装应用程序这样的东西。我只需要一种干净简单的方法来在 RaspberryPi 上
关闭。这个问题不符合Stack Overflow guidelines .它目前不接受答案。 这个问题似乎不是关于 a specific programming problem, a softwar
关闭。这个问题是off-topic .它目前不接受答案。 想改进这个问题吗? Update the question所以它是on-topic用于堆栈溢出。 关闭 10 年前。 Improve thi
有什么方法可以覆盖现有的源代码,我应该用 PyQt、PyGTK、Java 等从头开始构建吗? 最佳答案 如果您指的是软件本身而不是它所连接的存储库,那么自定义应用程序的方法就是 fork 项目。据我所
我的情况是:我在一个磁盘上安装了两个 linux。我将第一个安装在/dev/sda1 中,然后在/dev/sda2 中安装第二个然后我运行第一个系统,我写了一个脚本来在第一个系统运行时更新它。
我在 i2c-0 总线上使用地址为 0x3f 的系统监视器设备。该设备在设备树中配置有 pmbus 驱动程序。 问题是,加载 linux 内核时,这个“Sysmon”设备没有供电。因此,当我在总线 0
关闭。这个问题是off-topic .它目前不接受答案。 想改进这个问题吗? Update the question所以它是on-topic用于堆栈溢出。 关闭 11 年前。 Improve thi
我正试图在 linux 模块中分配一大块内存,而 kalloc 做不到。 我知道唯一的方法是使用 alloc_bootmem(unsigned long size) 但我只能从 linux 内核而不是
关闭。这个问题不符合Stack Overflow guidelines .它目前不接受答案。 这个问题似乎不是关于 a specific programming problem, a softwar
我有 .sh 文件来运行应用程序。在该文件中,我想动态设置服务器名称,而不是每次都配置。 我尝试了以下方法,它在 CentOS 中运行良好。 nohup /voip/java/jdk1.8.0_71/
我是在 Linux 上开发嵌入式 C++ 程序的新手。我有我的 Debian 操作系统,我在其中开发和编译了我的 C++ 项目(一个简单的控制台进程)。 我想将我的应用程序放到另一个 Debian 操
关闭。这个问题需要多问focused 。目前不接受答案。 想要改进此问题吗?更新问题,使其仅关注一个问题 editing this post . 已关闭 4 年前。 Improve this ques
我使用4.19.78版本的稳定内核,我想找到带有企鹅二进制数据的C数组。系统启动时显示。我需要在哪里搜索该内容? 我在 include/linux/linux_logo.h 文件中只找到了一些 Log
我知道可以使用 gdb 的服务器模式远程调试代码,我知道可以调试针对另一种架构交叉编译的代码,但是是否可以更进一步,从远程调试 Linux 应用程序OS X 使用 gdbserver? 最佳答案 当然
是否有任何可能的方法来运行在另一个 Linux 上编译的二进制文件?我知道当然最简单的是在另一台机器上重建它,但假设我们唯一能得到的是一个二进制文件,那么这可能与否? (我知道这可能并不容易,但我只是
我是一名优秀的程序员,十分优秀!