質問編集履歴

追記の作成

2019/05/19 01:15

投稿

syoukera

スコア10

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -165,3 +165,207 @@
 ```
+# 追記20190519
+Python上で整数型の値を管理し，Fortranに渡して装置番号として使用するということを考えました．テストのためにPythonからOpen文を含むコードを書いて実行を．結果は以下．
+|実行環境|結果|
+|:---|:---|
+|Pythonで装置番号を定義，Fortranに渡す|Segmentation fault|
+|Pythonから引数なしでFortranを呼ぶ，Fortranで装置番号を定義|実行可能|
+Fortran上で整数型を定義して渡すうまくいきますが，Pythonで管理しようとするとOPENの過程で怒られました．
+参考に使用したコードを置いておきます．
+```python
+import ctypes
+import numpy as np
+# device number as argument
+def call_openfiles(readDN, writeDN):
+    f = np.ctypeslib.load_library("libfort.so", ".")
+    f.openfiles_.argtypes = [
+        ctypes.c_int32,
+        ctypes.c_int32
+        ]
+    f.openfiles_.restype  = ctypes.c_void_p
+    f.openfiles_(readDN, writeDN)
+# non argument
+def call_fortran():
+    f = np.ctypeslib.load_library("libfort.so", ".")
+    f.main_.restype  = ctypes.c_void_p
+    f.main_()
+# 装置番号を定義して引数として渡すとき
+readDN  = 10
+writeDN = 11
+print("**call fortran from python with device number as argument")
+call_openfiles(readDN, writeDN)
+# 引数なしで渡すとき
+# print("**call fortran from python with non argument")
+# call_fortran()
+```
+```fortran90
+subroutine main()
+    implicit none
+    integer :: readDN = 10
+    integer :: writeDN = 11
+    call openFiles(readDN, writeDN)
+    call rwNumbers(readDN, writeDN)
+    call closeFiles(readDN, writeDN)
+end subroutine main
+subroutine rwNumbers(readDN, writeDN)
+    implicit none
+    real,dimension(30)::numbers
+    real temp
+    integer i, j
+    integer readDN
+    integer writeDN
+    ! read numbers
+    do i = 1, 30
+        read(readDN, *) numbers(i)
+    end do
+    close(10)
+    ! baggle sort
+    do i = 2, 30
+        do j = 2, 32-i
+            if (numbers(j) < numbers(j-1)) then
+                temp = numbers(j)
+                numbers(j) = numbers(j-1)
+                numbers(j-1) = temp
+            end if
+        end do
+    end do
+    ! print numbers
+    print *, "Sorted numbers is written in sortedNumbers.txt"
+    do i=1, 30
+        write(writeDN,*) numbers(i)
+    end do
+end subroutine
+subroutine openFiles(readDN, writeDN)
+    implicit none
+    integer(4),intent(in) :: readDN
+    integer(4),intent(in) :: writeDN
+    print *, "Open files from fortran"
+    open(readDN , file='numbers.txt', status='old')
+    open(writeDN , file='sortedNumbers.txt', status='replace')
+end subroutine
+subroutine closeFiles(readDN, writeDN)
+    implicit none
+    integer readDN, writeDN
+    close(readDN)
+    close(writeDN)
+end subroutine
+```