translate debug

first make:

 

administrator@ubuntu:~/dan_code/translate$ make
rm -f ./aaa
test -d ./Obj/ || mkdir ./Obj/
gfortran -g -c -o Obj/reorder.o reorder.f90
gfortran -g -c -o Obj/testtededg.o testtededg.f90
gfortran -g -o ./aaa ./Obj/main.o ./Obj/readdata.o ./Obj/nodetoedge.o ./Obj/convert.o ./Obj/reorder.o ./Obj/testtededg.o 
./Obj/reorder.o: In function `calrec':
/home/administrator/dan_code/translate/reorder.f90:212: undefined reference to `eof_'
./Obj/testtededg.o: In function `testtededg':
/home/administrator/dan_code/translate/testtededg.f90:43: undefined reference to `eof_'
/home/administrator/dan_code/translate/testtededg.f90:49: undefined reference to `eof_'
/home/administrator/dan_code/translate/testtededg.f90:55: undefined reference to `eof_'
/home/administrator/dan_code/translate/testtededg.f90:61: undefined reference to `eof_'
collect2: ld returned 1 exit status
make: *** [aaa] 错误 1

eof文件尾函数不适用的缘故，参照文末fortran读文件时判断文件结束的两种方法

1. 修改reorder.f90

do while (.true.) n=n+1 read(fino,*,iostat=stat) (a(i),i=1,opno) if(stat <= 0) goto 100 end do 100 no=n

出错结果no=0

修改为：

do while (.true.) n=n+1 read(fino,*,end=100) (a(i),i=1,opno) end do 100 no=n

结果为：

*******reorder the edges******* nted 173 At line 95 of file reorder.f90 (unit = 60, file = 'cubic.tedfulledge') Fortran runtime error: End of file

nted为173大于172，读下一个文件的时候超过文件行数

修改为no=n-1，结果为：

*******reorder the edges******* nted 172 nsrcedge= 0 nbdedge= 156 total free edge= 156 *******test ted edge******* Nted 172 Nedge 288 173 173 289 289 段错误

 

段错误定位于testtededg.f90中Line 69循环当中。

 

用原Translate_code，做类似本文修改，没有出现错误。执行结果为：

 

shanon@ubuntu:~/dan_code/Translate_code$ ./translate !!!!!!!!!!!!!!!!!!!!!now second step!!!!!!!!!!!!! *******begin to generate .facet,.node******** 65 172 the lower bd 0.0000000 0.0000000 0.0000000 the upper bd 1.0000000 0.75000000 0.50000000 *******begin to generate .fulledge,.tedfulledge******* Number of nodes = 65 Number of teds = 172 Number of edges = 288 edgemax,edgmin 0.48369488 0.14733800 213 48 *******begin to generate .patch,.tedpatch,bdpatch..******* srcbd= 1000000.0 1000000.0 1000000.0 -1000000.0 -1000000.0 -1000000.0 Number of nodes = 65 Number of teds = 172 1000000.0 -1000000.0 1000000.0 -1000000.0 1000000.0 -1000000.0 Number of active facets = 292 *******reorder the edges******* nted= 172 nsrcedge= 0 nbdedge= 156 total free edge=nbd+nsrc 156 *******test ted edge******* Nted, Nedge 172 288 testing tededg, tedfulledg 172 172 288 288 STOP ok,go on to tdfem_pml shanon@ubuntu:~/dan_code/Translate_code$

2. 修改testtededg.f90

do while(.not.eof(300)) ip1=ip1+1 read(300,*) (tededg(i,ip1), i=1,6) end do ip2=0 do while(.not.eof(400)) ip2=ip2+1 read(400,*) (tedfull(i,ip2), i=1,6) end do ip3=0 do while(.not.eof(500)) ip3=ip3+1 read(500,*) (testedg(i,ip3), i=1,2) ! .edge file end do ip4=0 do while(.not.eof(600)) ip4=ip4+1 read(600,*) (full(i,ip4), i=1,2) ! .fulledge file end do close(300);close(400);close(500);close(600)

改为

 

do while(.true.) ip1=ip1+1 read(300,*,end=10) (tededg(i,ip1), i=1,6) end do 10 ip2=0 do while(.true.) ip2=ip2+1 read(400,*,end=20) (tedfull(i,ip2), i=1,6) end do 20 ip3=0 do while(.true.) ip3=ip3+1 read(500,*,end=30) (testedg(i,ip3), i=1,2) ! .edge file end do 30 ip4=0 do while(.true.) ip4=ip4+1 read(600,*,end=40) (full(i,ip4), i=1,2) ! .fulledge file end do 40 close(300);close(400);close(500);close(600)

 

 

 

 

Reference:

fortran读文件时判断文件结束的两种方法

1，

...

open(11,file='example')

do while(.true.)

  read(11,*,iostat=status)record

  if(status.le.0)goto 100

  print*,record

enddo

100 close(11)

...

 

2,

...

open(11,file='example')

do while(.true.)

  read(11,*,end=100)record

  print*,record

enddo

100 close(11)