- android - 多次调用 OnPrimaryClipChangedListener
- android - 无法更新 RecyclerView 中的 TextView 字段
- android.database.CursorIndexOutOfBoundsException : Index 0 requested, 光标大小为 0
- android - 使用 AppCompat 时,我们是否需要明确指定其 UI 组件(Spinner、EditText)颜色
我使用 Fortran 来补充 Python,并且在一些方法中我使用 Python 方法作为 Fortran 子例程中的回调。一切似乎都正常,直到我将一个数组输入回调函数,如下所示。
RECURSIVE SUBROUTINE RECURSIVE_CURVE_SUBDIVISION(CPW, N, TOL, FUNC1)
IMPLICIT NONE
!F2PY INTENT(IN) CPW, N, TOL
!F2PY DEPEND(N) CPW
!F2PY (CALLBACK) FUNC1
INTEGER, INTENT(IN) :: N
DOUBLE PRECISION, INTENT(IN) :: CPW(0:N, 0:3), TOL
INTEGER :: I
DOUBLE PRECISION :: QP(0:N, 0:2), LP, LC, TEMP, &
AW(0:N, 0:3), BW(0:N, 0:3), V(0:2)
EXTERNAL :: FUNC1
DO I = 0, N
QP(I, :) = CPW(I, 0:2) / CPW(I, 3)
END DO
LP = 0.0D0
DO I = 0, N - 1
V = QP(I + 1, :) - QP(I, :)
CALL NORM(V, TEMP)
LP = LP + TEMP
END DO
V = QP(N, :) - QP(0, :)
CALL NORM(V, LC)
IF (ABS(LP - LC) .LE. TOL) THEN
CALL FUNC1(CPW, QP, LC, LP) !<-- here is the problem
! CALL FUNC1(LC, LP) !<-- this works
! CALL FUNC1(CPW=CPW, QP=QP, LC=LC, LP=LP)
! Added bonus if anyone can figure out how to use keyword arguements in
! the callback. For cleanliness, I'm trying to use func1(**kwargs) in Python.
ELSE
CALL SPLIT_BEZIER_CURVE(CPW, N, 0.50D0, AW, BW)
CALL RECURSIVE_CURVE_SUBDIVISION(AW, N, TOL / 2.0D0, FUNC1)
CALL RECURSIVE_CURVE_SUBDIVISION(BW, N, TOL / 2.0D0, FUNC1)
END IF
END SUBROUTINE RECURSIVE_CURVE_SUBDIVISION
这是尝试使用 f2py(使用 gfortran)编译时的一些输出:
warning C4244: '=' : conversion from 'npy_intp' to 'npy_int', possible loss of data
warning C4244: '=' : conversion from 'npy_intp' to 'npy_int', possible loss of data
warning C4244: '+=' : conversion from 'Py_ssize_t' to 'int', possible loss of data
warning C4244: '+=' : conversion from 'Py_ssize_t' to 'int', possible loss of data
warning C4244: '=' : conversion from 'Py_ssize_t' to 'int', possible loss of data
warning C4244: '=' : conversion from 'Py_ssize_t' to 'int', possible loss of data
warning C4244: '=' : conversion from 'Py_ssize_t' to 'int', possible loss of data
error C2065: 'n' : undeclared identifier
error C2065: 'n' : undeclared identifier
该模块本身可以很好地使用 gfortran 进行编译。我想我在 !F2PY 部分中没有足够的信息,但还没有弄清楚我缺少什么。
非常感谢任何提示。
更新 1:
所以我注意到我可以返回一个 1 x n 数组,但是 m x n 返回虚假结果。例如,我可以执行 CALL FUNC1(V) 并返回 1 x 3 数组 V 并将其打印到屏幕上(回调函数 FUNC1 现在只是打印到屏幕上进行测试)。当我用 CP 代替 V 时,它给出了上面显示的警告并且不会编译,所以它与数组的形状有关?
我不记得是在哪里看到的,但是如果我将子例程顶部附近的语句修改为:
!F2PY INTENT(IN) N, CPW, TOL
!F2PY DEPEND(N) CPW
!F2PY (CALLBACK) FUNC1
!F2PY CALL FUNC1(CP)
EXTERNAL :: FUNC1
它将编译并运行,但回调的输出(只是将数组打印到屏幕上)是伪造的。它是一个单一的 float ,每次迭代的幅度变化很大。某种段错误?
最佳答案
你必须在调用函数时传递 N
以便它知道数组有多大(你不必在 Python 中处理这个参数,但它在 C 级别需要它).所以把电话改成
CALL FUNC1(CPW, QP, LC, LP,N)
如果您不这样做,但查看生成的 .pyf 签名文件(f2py -m thingy -h thingy.pyf thingy.f90
)相关部分(它是自动生成的签名)是
python module recursive_curve_subdivision__user__routines
interface recursive_curve_subdivision_user_interface
subroutine func1(cpw,qp) ! in :thingy:thingy.f90:recursive_curve_subdivision:unknown_interface
double precision dimension(n + 1,4),intent(in),depend(n) :: cpw
double precision dimension(n + 1,3) :: qp
end subroutine func1
end interface recursive_curve_subdivision_user_interface
end python module recursive_curve_subdivision__user__routines
(请注意,为了简单起见,我使用仅包含两个参数的稍微精简的函数调用对其进行了测试,因此它与您的代码不完全匹配)。您会注意到它取决于所有数组大小的 N
,但您永远不会传递 N
。如果您确实添加了 N
作为参数,它知道大小并且 Python 的打印工作正常。
关于python - f2py 在 Fortran 的 Python 回调函数中使用数组,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/34005860/
C语言sscanf()函数:从字符串中读取指定格式的数据 头文件: ?
最近,我有一个关于工作预评估的问题,即使查询了每个功能的工作原理,我也不知道如何解决。这是一个伪代码。 下面是一个名为foo()的函数,该函数将被传递一个值并返回一个值。如果将以下值传递给foo函数,
CStr 函数 返回表达式,该表达式已被转换为 String 子类型的 Variant。 CStr(expression) expression 参数是任意有效的表达式。 说明 通常,可以
CSng 函数 返回表达式,该表达式已被转换为 Single 子类型的 Variant。 CSng(expression) expression 参数是任意有效的表达式。 说明 通常,可
CreateObject 函数 创建并返回对 Automation 对象的引用。 CreateObject(servername.typename [, location]) 参数 serv
Cos 函数 返回某个角的余弦值。 Cos(number) number 参数可以是任何将某个角表示为弧度的有效数值表达式。 说明 Cos 函数取某个角并返回直角三角形两边的比值。此比值是
CLng 函数 返回表达式,此表达式已被转换为 Long 子类型的 Variant。 CLng(expression) expression 参数是任意有效的表达式。 说明 通常,您可以使
CInt 函数 返回表达式,此表达式已被转换为 Integer 子类型的 Variant。 CInt(expression) expression 参数是任意有效的表达式。 说明 通常,可
Chr 函数 返回与指定的 ANSI 字符代码相对应的字符。 Chr(charcode) charcode 参数是可以标识字符的数字。 说明 从 0 到 31 的数字表示标准的不可打印的
CDbl 函数 返回表达式,此表达式已被转换为 Double 子类型的 Variant。 CDbl(expression) expression 参数是任意有效的表达式。 说明 通常,您可
CDate 函数 返回表达式,此表达式已被转换为 Date 子类型的 Variant。 CDate(date) date 参数是任意有效的日期表达式。 说明 IsDate 函数用于判断 d
CCur 函数 返回表达式,此表达式已被转换为 Currency 子类型的 Variant。 CCur(expression) expression 参数是任意有效的表达式。 说明 通常,
CByte 函数 返回表达式,此表达式已被转换为 Byte 子类型的 Variant。 CByte(expression) expression 参数是任意有效的表达式。 说明 通常,可以
CBool 函数 返回表达式,此表达式已转换为 Boolean 子类型的 Variant。 CBool(expression) expression 是任意有效的表达式。 说明 如果 ex
Atn 函数 返回数值的反正切值。 Atn(number) number 参数可以是任意有效的数值表达式。 说明 Atn 函数计算直角三角形两个边的比值 (number) 并返回对应角的弧
Asc 函数 返回与字符串的第一个字母对应的 ANSI 字符代码。 Asc(string) string 参数是任意有效的字符串表达式。如果 string 参数未包含字符,则将发生运行时错误。
Array 函数 返回包含数组的 Variant。 Array(arglist) arglist 参数是赋给包含在 Variant 中的数组元素的值的列表(用逗号分隔)。如果没有指定此参数,则
Abs 函数 返回数字的绝对值。 Abs(number) number 参数可以是任意有效的数值表达式。如果 number 包含 Null,则返回 Null;如果是未初始化变量,则返回 0。
FormatPercent 函数 返回表达式,此表达式已被格式化为尾随有 % 符号的百分比(乘以 100 )。 FormatPercent(expression[,NumDigitsAfterD
FormatNumber 函数 返回表达式,此表达式已被格式化为数值。 FormatNumber( expression [,NumDigitsAfterDecimal [,Inc
我是一名优秀的程序员,十分优秀!