4

我正在使用 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

AddVecgonum 实现的方法:

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满足Vector来自 gonum 的接口,因此将这个内部字段传递给方法工作正常,如上面的示例所示。

4

2 回答 2

1

在 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具有Vectorstruct 的值类型参数,但是您将指向Vector字段的指针传递为

x, y := NewVector(3, []float64{1, 2, 3}), NewVector(3, []float64{4, 5, 6}) 

在上面的代码x,y中是从返回的指针类型NewVector

x.AddVec(x, y)
于 2018-08-21T14:34:56.863 回答
0

问题是在进行阴影检测时会比较指针。这是国际海事组织的一个错误。我刚刚提交了https://github.com/gonum/gonum/issues/945

于 2019-03-31T21:21:07.877 回答