我正在尝试用Go编写一个包,使用“通用”类型计算方程。具体来说,我想实现龙格库塔五阶逼近法。
这种逼近法仅使用未知函数y在t0时刻的值、起始时间t0、步长h和形式为dy/dt = g(t,y)的微分方程dgl来计算点t0 + h处y的值。
这种逼近法在使用标量类型和向量(甚至矩阵)时表现完全相同。更一般地说:它适用于任何可以与相同类型的值相加/相减并且可以通过标量(我使用float64)进行缩放的东西。
因此,我尝试将其表示为Go接口:
但是当我尝试“实现”这个接口时,由于参数类型的原因,我遇到了麻烦:
这导致我出现错误。
这种逼近法仅使用未知函数y在t0时刻的值、起始时间t0、步长h和形式为dy/dt = g(t,y)的微分方程dgl来计算点t0 + h处y的值。
这种逼近法在使用标量类型和向量(甚至矩阵)时表现完全相同。更一般地说:它适用于任何可以与相同类型的值相加/相减并且可以通过标量(我使用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中不可能。