我正在使用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: %v\n", &u, u)
fmt.Printf("[V - NewVecDense]\tADDRESS: %v, VALUE: %v\n", &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: %v\n", &x, x)
fmt.Printf("[Y - NewVector]\tADDRESS: %v, VALUE: %v\n", &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
是一种方法implemented by 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,提升的方法包含在结构的方法集中,如下所示:
问题是你将指针类型参数传递给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