func sortEigenVec(D, U, V, C *cmat.FloatMatrix, updown int) { var sD, m0, m1 cmat.FloatMatrix N := D.Len() for k := 0; k < N-1; k++ { sD.SubVector(D, k, N-k) pk := vecMinMax(&sD, -updown) if pk != 0 { t0 := D.GetAt(k) D.SetAt(k, D.GetAt(pk+k)) D.SetAt(k+pk, t0) if U != nil { m0.Column(U, k) m1.Column(U, k+pk) blasd.Swap(&m1, &m0) } if V != nil { m0.Row(V, k) m1.Row(V, k+pk) blasd.Swap(&m1, &m0) } if C != nil { m0.Column(C, k) m1.Column(C, k+pk) blasd.Swap(&m1, &m0) } } } }