我正在使用gonum来执行一些线性代数计算。扩展原始mat.VecDense
结构后,我在自身上应用方法时出现"坏区域:相同"的恐慌。当我使用 gonum 提供的原始设置时,不会发生此错误。
这是我的实现:
type Vector struct {
mat.VecDense
}
func NewVector(n int, data []float64) *Vector {
return &Vector{*mat.NewVecDense(n, data)}
}
我正在使用以下代码片段对其进行测试:
func main() {
u, v := mat.NewVecDense(3, []float64{1, 2, 3}), mat.NewVecDense(3, []float64{4, 5, 6})
fmt.Printf("[U - NewVecDense]tADDRESS: %v, VALUE: %vn", &u, u)
fmt.Printf("[V - NewVecDense]tADDRESS: %v, VALUE: %vn", &v, v)
u.AddVec(u, v)
fmt.Println("-------------------------")
x, y := NewVector(3, []float64{1, 2, 3}), NewVector(3, []float64{4, 5, 6})
fmt.Printf("[X - NewVector]tADDRESS: %v, VALUE: %vn", &x, x)
fmt.Printf("[Y - NewVector]tADDRESS: %v, VALUE: %vn", &y, y)
x.AddVec(x, y)
fmt.Println(x)
}
虽然第一次添加执行正常,但第二次添加失败:
[U - NewVecDense] ADDRESS: 0xc42000c028, VALUE: &{{[1 2 3] 1} 3}
[V - NewVecDense] ADDRESS: 0xc42000c030, VALUE: &{{[4 5 6] 1} 3}
-------------------------
[X - NewVector] ADDRESS: 0xc42000c040, VALUE: &{{{[1 2 3] 1} 3}}
[Y - NewVector] ADDRESS: 0xc42000c048, VALUE: &{{{[4 5 6] 1} 3}}
panic: mat: bad region: identical
AddVec
是 gonum 实现的方法:
func (v *VecDense) AddVec(a, b Vector)
为什么会这样,手是什么正确的实施方法?
编辑:
多亏了@Himanshu我设法解决了这个问题。
我为我正在使用的每个方法创建了直通方法,传递结构的正确级别:
type Vector struct {
*mat.VecDense
}
func NewVector(n int, data []float64) Vector {
return Vector{mat.NewVecDense(n, data)}
}
func (v *Vector) AddVec(a, b Vector) {
v.VecDense.AddVec(a.VecDense, b.VecDense)
}
func (v *Vector) SubVec(a, b Vector) {
v.VecDense.SubVec(a.VecDense, b.VecDense)
}
func (v *Vector) ScaleVec(alpha float64, a Vector) {
v.VecDense.ScaleVec(alpha, a.VecDense)
}
func (v *Vector) AddScaledVec(a Vector, alpha float64, b Vector) {
v.VecDense.AddScaledVec(a.VecDense, alpha, b.VecDense)
}
此外 - 我不确定这是否是正确的方法 - 我还将NewVector
的返回类型从指针更改为值,因为它无论如何都包含指向mat.VecDense
的指针。请注意,*mat.VecDense
满足 gonum 的Vector
接口,因此将此内部字段传递给方法工作正常,如上例所示。
在 Golang 中,它被描述为提升的方法为
升级字段的行为类似于结构的普通字段,只是它们 不能用作结构的复合文本中的字段名称。
给定结构类型 S 和定义的类型 T,提升的方法包含在结构的方法集中,如下所示:
如果 S 包含- 嵌入字段 T,则 S 和 *S 的方法集都包含具有接收方 T 的提升方法。*S 的方法集还包括带有接收器 *T 的提升方法。
- 如果 S 包含嵌入字段 *T,则 S 和 *S 的方法集都包括具有接收器 T 或 *T 的提升方法。
问题是您正在将指针类型参数传递给AddVec
函数。但是您在第二种情况下使用指针类型字段。
func (v *VecDense) AddVec(a, b Vector)
需要注意的另一件事是,AddVec
具有Vector
结构的值类型参数,但您将指向Vector
字段的指针传递为
x, y := NewVector(3, []float64{1, 2, 3}), NewVector(3, []float64{4, 5, 6})
在上面的代码中x,y
是从NewVector
返回的指针类型
x.AddVec(x, y)
问题是在进行阴影检测时会比较指针。这是一个错误 IMO。我刚刚提交了 https://github.com/gonum/gonum/issues/945