kshedden / dimred Goto Github PK

	// q3 is a symmetric matrix but not of type Symmetric, would
	// be better to use EigenSym here but hard to convert.
	var ei mat64.Eigen
	ok := ei.Factorize(q3, true)

y := cm.Data.GetPos(cm.ypos).([]float64)
	var mn [][]float64
	var ny []int

        // Calculate the sums and sample sizes
        for j, xp := range cm.xpos {
		x := cm.Data.GetPos(xp).([]float64)
                fmt.Printf("current x: %v\n", x)
                for i := 0; i < len(y); i++ {
                        yi := int(y[i])
                        mn = growf(mn, yi, p)
			mn[yi][j] += x[i]
                        if j == 0 {
                                ny = growi(ny, yi)
	                        ny[yi]++
			}
                }
        }