【发布时间】:2013-01-04 03:18:35
【问题描述】:
我正在尝试在 Go 中编写一个使用“通用”类型计算方程的包。具体来说,我想实现runge kutta 5 逼近。
这个近似值仅使用t0处的y的值、开始时间t0、步宽h和一个微分方程dgl,其形式为dy/dt = g(t,y),其中g 是某个函数。
当使用标量类型时,这种近似的行为与使用向量(甚至矩阵)时的行为完全相同。更一般地说:它适用于可以添加/减去相同类型的值并且可以通过标量缩放的所有内容(我使用float64)
所以我尝试将其表达为 Go 接口:
type Numeric interface {
Add(rhs Numeric) Numeric
Sub(rhs Numeric) Numeric
Mul(rhs float64) Numeric
}
但是当我尝试“实现”这个接口时,由于参数类型,我遇到了麻烦:
type Vec6F struct {
x, y, z float64
vx, vy, vz float64
}
func (lhs *Vec6F) Add(rhs *Vec6F) rk5.Numeric {
result := new(Vec6F)
result.x = lhs.x + rhs.x
result.y = lhs.y + rhs.y
result.z = lhs.z + rhs.z
result.vx = lhs.vx + rhs.vx
result.vy = lhs.vy + rhs.vy
result.vz = lhs.vz + rhs.vz
return result
}
这给了我错误
cannot use result (type *Vec6F) as type rk5.Numeric in return argument:
*Vec6F does not implement rk5.Numeric (wrong type for Add method
have Add(*Vec6F) rk5.Numeric
want Add(rk5.Numeric) rk5.Numeric
一方面,这对我来说绝对符合逻辑(因为 rhs 可能是另一个实现 Numeric 的对象)
但另一方面:我如何在 Go 中表达类似的东西?在 C++ 中,我可以改用运算符重载,但这在 go 中是不可能的。
【问题讨论】:
-
“off-topic”:求解的微分方程是描述(简单)火星轨道的微分方程。对于 quick'n'dirty 完整代码,请参阅play.golang.org/p/GwehylBbLK