本文整理汇总了Golang中github.com/hrautila/matrix.FloatWithValue函数的典型用法代码示例。如果您正苦于以下问题:Golang FloatWithValue函数的具体用法?Golang FloatWithValue怎么用?Golang FloatWithValue使用的例子?那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了FloatWithValue函数的20个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于我们的系统推荐出更棒的Golang代码示例。
示例1: MakeData
func MakeData(M, N, P int, randomData, diagonal bool) (A, B, C *matrix.FloatMatrix) {
if diagonal && M != N {
diagonal = false
fmt.Printf("cannot make diagonal if B.rows != B.cols\n")
}
if randomData {
A = matrix.FloatNormal(M, P)
if diagonal {
d := matrix.FloatNormal(P, 1)
B := matrix.FloatDiagonal(P, 0.0)
B.SetIndexesFromArray(d.FloatArray(), matrix.DiagonalIndexes(B)...)
} else {
B = matrix.FloatNormal(P, N)
}
} else {
A = matrix.FloatWithValue(M, P, 1.0)
if diagonal {
B = matrix.FloatDiagonal(P, 1.0)
} else {
B = matrix.FloatWithValue(P, N, 1.0)
}
}
C = matrix.FloatZeros(M, N)
return
}
开发者ID:hrautila,项目名称:mperf,代码行数:27,代码来源:helpers.go
示例2: TestVectors
func TestVectors(t *testing.T) {
a := matrix.FloatWithValue(1, 5, 1.0)
b := matrix.FloatWithValue(1, 5, 2.0)
c := matrix.FloatWithValue(5, 1, 2.0)
ar := a.FloatArray()
br := b.FloatArray()
cr := c.FloatArray()
v := DDot(ar, br, 1.0, a.LeadingIndex(), b.LeadingIndex(), a.NumElements())
t.Logf("a*b = %.1f\n", v)
v = DDot(cr, br, 1.0, 1, b.LeadingIndex(), c.NumElements())
t.Logf("c*b = %.1f\n", v)
v = DNorm2(br, 1, b.NumElements())
t.Logf("norm2(b) = %.1f\n", v)
b0 := matrix.FloatNew(1, 5, []float64{-1.0, 0.0, 2.0, -3.0, 5.0})
br = b0.FloatArray()
ix := DIAMax(br, 1, b0.NumElements())
t.Logf("iamax = %d: %.1f %v\n", ix, b0.GetIndex(ix), b0)
b0 = matrix.FloatNew(1, 5, []float64{-8.0, 0.0, 2.0, -3.0, 5.0})
br = b0.FloatArray()
ix = DIAMax(br, 1, b0.NumElements())
t.Logf("iamax = %d: %.1f %v\n", ix, b0.GetIndex(ix), b0)
b0 = matrix.FloatNew(1, 5, []float64{-8.0, 0.0, 2.0, -9.0, 5.0})
br = b0.FloatArray()
ix = DIAMax(br, 1, b0.NumElements())
t.Logf("iamax = %d: %.1f %v\n", ix, b0.GetIndex(ix), b0)
}
开发者ID:hrautila,项目名称:matops,代码行数:29,代码来源:calgo_test.go
示例3: trmmTest
func trmmTest(t *testing.T, A *matrix.FloatMatrix, flags Flags, nb int) bool {
var B0 *matrix.FloatMatrix
N := A.Cols()
S := 0
E := A.Cols()
side := linalg.OptLeft
if flags&RIGHT != 0 {
B0 = matrix.FloatWithValue(2, A.Rows(), 2.0)
side = linalg.OptRight
E = B0.Rows()
} else {
B0 = matrix.FloatWithValue(A.Rows(), 2, 2.0)
E = B0.Cols()
}
B1 := B0.Copy()
trans := linalg.OptNoTrans
if flags&TRANSA != 0 {
trans = linalg.OptTransA
}
diag := linalg.OptNonUnit
if flags&UNIT != 0 {
diag = linalg.OptUnit
}
uplo := linalg.OptUpper
if flags&LOWER != 0 {
uplo = linalg.OptLower
}
blas.TrmmFloat(A, B0, 1.0, uplo, diag, trans, side)
if A.Rows() < 8 {
//t.Logf("..A\n%v\n", A)
t.Logf(" BLAS B0:\n%v\n", B0)
}
Ar := A.FloatArray()
Br := B1.FloatArray()
if nb != 0 {
DTrmmBlk(Br, Ar, 1.0, flags, B1.LeadingIndex(), A.LeadingIndex(),
N, S, E, nb)
} else {
DTrmmUnblk(Br, Ar, 1.0, flags, B1.LeadingIndex(), A.LeadingIndex(),
N, S, E, 0)
}
result := B0.AllClose(B1)
t.Logf(" B0 == B1: %v\n", result)
if A.Rows() < 8 {
t.Logf(" DTrmm B1:\n%v\n", B1)
}
return result
}
开发者ID:hrautila,项目名称:matops,代码行数:51,代码来源:ttrm_test.go
示例4: trmvTest
func trmvTest(t *testing.T, A *matrix.FloatMatrix, flags Flags, nb int) bool {
N := A.Cols()
//S := 0
//E := A.Cols()
X0 := matrix.FloatWithValue(A.Rows(), 1, 2.0)
X1 := X0.Copy()
trans := linalg.OptNoTrans
if flags&TRANS != 0 {
trans = linalg.OptTrans
}
diag := linalg.OptNonUnit
if flags&UNIT != 0 {
diag = linalg.OptUnit
}
uplo := linalg.OptUpper
if flags&LOWER != 0 {
uplo = linalg.OptLower
}
blas.TrmvFloat(A, X0, uplo, diag, trans)
Ar := A.FloatArray()
Xr := X1.FloatArray()
if nb == 0 {
DTrimvUnblkMV(Xr, Ar, flags, 1, A.LeadingIndex(), N)
}
result := X0.AllClose(X1)
t.Logf(" X0 == X1: %v\n", result)
if !result && A.Rows() < 8 {
t.Logf(" BLAS TRMV X0:\n%v\n", X0)
t.Logf(" DTrmv X1:\n%v\n", X1)
}
return result
}
开发者ID:hrautila,项目名称:matops,代码行数:35,代码来源:ttrm_test.go
示例5: _TestSolveBlockedSmall
func _TestSolveBlockedSmall(t *testing.T) {
L := matrix.FloatMatrixFromTable(lower5, matrix.RowOrder)
N := L.Rows()
nb := 4
U := L.Transpose()
X0 := matrix.FloatWithValue(L.Rows(), 1, 1.0)
X1 := X0.Copy()
xsum := 0.0
for i := 0; i < N; i++ {
xsum += float64(i)
X0.Add(xsum, i)
X1.Add(xsum, -(i + 1))
}
t.Logf("-- SOLVE LOWER, NON-UNIT ---\n")
solveMVTest(t, L, X0.Copy(), LOWER, N, nb)
t.Logf("-- SOLVE UPPER, NON-UNIT ---\n")
solveMVTest(t, U, X1.Copy(), UPPER, N, nb)
t.Logf("-- SOLVE LOWER, UNIT ---\n")
solveMVTest(t, L, X0.Copy(), LOWER|UNIT, N, nb)
t.Logf("-- SOLVE UPPER, UNIT ---\n")
solveMVTest(t, U, X1.Copy(), UPPER|UNIT, N, nb)
}
开发者ID:hrautila,项目名称:matops,代码行数:26,代码来源:tsolve_test.go
示例6: _TestLDLnoPiv
func _TestLDLnoPiv(t *testing.T) {
N := 42
nb := 8
A0 := matrix.FloatUniform(N, N)
A := matrix.FloatZeros(N, N)
Mult(A, A0, A0, 1.0, 1.0, TRANSB)
B := matrix.FloatNormal(A.Rows(), 2)
w := matrix.FloatWithValue(A.Rows(), 2, 1.0)
// B0 = A*B
B0 := B.Copy()
nb = 2
L, _ := DecomposeLDLnoPiv(A.Copy(), w, LOWER, nb)
Mult(B0, A, B, 1.0, 0.0, NOTRANS)
SolveLDLnoPiv(B0, L, LOWER)
t.Logf("L*D*L.T: ||B - A*X||_1: %e\n", NormP(B0.Minus(B), NORM_ONE))
U, _ := DecomposeLDLnoPiv(A.Copy(), w, UPPER, nb)
Mult(B0, A, B, 1.0, 0.0, NOTRANS)
SolveLDLnoPiv(B0, U, UPPER)
t.Logf("U*D*U.T: ||B - A*X||_1: %e\n", NormP(B0.Minus(B), NORM_ONE))
}
开发者ID:hrautila,项目名称:matops,代码行数:26,代码来源:ldl_test.go
示例7: _TestMultMVTransA
func _TestMultMVTransA(t *testing.T) {
bM := 1000 * M
bN := 1000 * N
A := matrix.FloatNormal(bN, bM)
X := matrix.FloatWithValue(bN, 1, 1.0)
Y1 := matrix.FloatZeros(bM, 1)
Y0 := matrix.FloatZeros(bM, 1)
Ar := A.FloatArray()
Xr := X.FloatArray()
Y1r := Y1.FloatArray()
blas.GemvFloat(A, X, Y0, 1.0, 1.0, linalg.OptTrans)
DMultMV(Y1r, Ar, Xr, 1.0, 1.0, TRANSA, 1, A.LeadingIndex(), 1, 0, bN, 0, bM, 4, 4)
ok := Y0.AllClose(Y1)
t.Logf("Y0 == Y1: %v\n", ok)
if !ok {
var y1, y0 matrix.FloatMatrix
Y1.SubMatrix(&y1, 0, 0, 5, 1)
t.Logf("Y1[0:5]:\n%v\n", y1)
Y0.SubMatrix(&y0, 0, 0, 5, 1)
t.Logf("Y0[0:5]:\n%v\n", y0)
}
}
开发者ID:hrautila,项目名称:matops,代码行数:25,代码来源:tmult_test.go
示例8: _TestMultSymmSmall
func _TestMultSymmSmall(t *testing.T) {
//bM := 5
bN := 7
bP := 7
Adata := [][]float64{
[]float64{1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0},
[]float64{0.0, 2.0, 2.0, 2.0, 2.0, 2.0, 2.0},
[]float64{0.0, 0.0, 3.0, 3.0, 3.0, 3.0, 3.0},
[]float64{0.0, 0.0, 0.0, 4.0, 4.0, 4.0, 4.0},
[]float64{0.0, 0.0, 0.0, 0.0, 5.0, 5.0, 5.0},
[]float64{0.0, 0.0, 0.0, 0.0, 0.0, 6.0, 6.0},
[]float64{0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 7.0}}
A := matrix.FloatMatrixFromTable(Adata, matrix.RowOrder)
B := matrix.FloatWithValue(bN, bP, 2.0)
C0 := matrix.FloatZeros(bN, bP)
C1 := matrix.FloatZeros(bN, bP)
Ar := A.FloatArray()
Br := B.FloatArray()
C1r := C1.FloatArray()
blas.SymmFloat(A, B, C0, 1.0, 1.0, linalg.OptUpper, linalg.OptRight)
DMultSymm(C1r, Ar, Br, 1.0, 1.0, UPPER|RIGHT, bN, A.LeadingIndex(), bN, bN, 0, bP, 0, bN, 2, 2, 2)
ok := C0.AllClose(C1)
t.Logf("C0 == C1: %v\n", ok)
if !ok {
t.Logf("A=\n%v\n", A)
t.Logf("blas: C=A*B\n%v\n", C0)
t.Logf("C1: C1 = A*X\n%v\n", C1)
}
}
开发者ID:hrautila,项目名称:matops,代码行数:33,代码来源:tmult_test.go
示例9: _TestRankSmall
func _TestRankSmall(t *testing.T) {
bM := 5
bN := 5
//bP := 5
Adata := [][]float64{
[]float64{1.0, 1.0, 1.0, 1.0, 1.0},
[]float64{2.0, 2.0, 2.0, 2.0, 2.0},
[]float64{3.0, 3.0, 3.0, 3.0, 3.0},
[]float64{4.0, 4.0, 4.0, 4.0, 4.0},
[]float64{5.0, 5.0, 5.0, 5.0, 5.0}}
A := matrix.FloatMatrixFromTable(Adata, matrix.RowOrder)
A0 := matrix.FloatMatrixFromTable(Adata, matrix.RowOrder)
X := matrix.FloatVector([]float64{1.0, 2.0, 3.0, 4.0, 5.0})
Y := matrix.FloatWithValue(bN, 1, 2.0)
Ar := A.FloatArray()
Xr := X.FloatArray()
Yr := Y.FloatArray()
blas.GerFloat(X, Y, A0, 1.0)
DRankMV(Ar, Xr, Yr, 1.0, A.LeadingIndex(), 1, 1, 0, bN, 0, bM, 4, 4)
ok := A0.AllClose(A)
t.Logf("A0 == A1: %v\n", ok)
if !ok {
t.Logf("blas ger:\n%v\n", A0)
t.Logf("A1: \n%v\n", A)
}
}
开发者ID:hrautila,项目名称:matops,代码行数:30,代码来源:trank_test.go
示例10: TestLDLlower
func TestLDLlower(t *testing.T) {
/*
Ldata := [][]float64{
[]float64{7.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0},
[]float64{7.0, 6.0, 0.0, 0.0, 0.0, 0.0, 0.0},
[]float64{7.0, 6.0, 5.0, 0.0, 0.0, 0.0, 0.0},
[]float64{7.0, 6.0, 5.0, 4.0, 0.0, 0.0, 0.0},
[]float64{7.0, 6.0, 5.0, 4.0, 6.0, 0.0, 0.0},
[]float64{7.0, 6.0, 5.0, 4.0, 3.0, 2.0, 0.0},
[]float64{7.0, 6.0, 5.0, 4.0, 3.0, 2.0, 1.0}}
A := matrix.FloatMatrixFromTable(Ldata, matrix.RowOrder)
N := A.Rows()
*/
N := 7
nb := 0
A0 := matrix.FloatUniform(N, N)
A := matrix.FloatZeros(N, N)
Mult(A, A0, A0, 1.0, 1.0, TRANSB)
B := matrix.FloatNormal(A.Rows(), 2)
B0 := B.Copy()
B1 := B.Copy()
Mult(B0, A, B, 1.0, 0.0, NOTRANS)
_, _, _ = B0, B1, A0
ipiv := make([]int, N, N)
L, _ := DecomposeLDL(A.Copy(), nil, ipiv, LOWER, 0)
//t.Logf("unblk: ipiv = %v\n", ipiv)
//t.Logf("unblk: L\n%v\n", L)
ApplyRowPivots(B, ipiv, FORWARD)
MultTrm(B, L, 1.0, LOWER|UNIT|TRANSA)
MultDiag(B, L, LEFT)
MultTrm(B, L, 1.0, LOWER|UNIT)
ApplyRowPivots(B0, ipiv, FORWARD)
t.Logf(" unblk: L*D*L.T %d pivots: ||A*B - L*D*L.T*B||_1: %e\n",
NumPivots(ipiv), NormP(B.Minus(B0), NORM_ONE))
t.Logf("pivots: %v\n", ipiv)
nb = 4
w := matrix.FloatWithValue(A.Rows(), nb, 1.0)
L, _ = DecomposeLDL(A.Copy(), w, ipiv, LOWER, nb)
//t.Logf("blk: ipiv = %v\n", ipiv)
//t.Logf("blk: L\n%v\n", L)
// B2 = A*B1 == A*B
B2 := B1.Copy()
Mult(B2, A, B1, 1.0, 0.0, NOTRANS)
ApplyRowPivots(B1, ipiv, FORWARD)
MultTrm(B1, L, 1.0, LOWER|UNIT|TRANSA)
MultDiag(B1, L, LEFT)
MultTrm(B1, L, 1.0, LOWER|UNIT)
ApplyRowPivots(B2, ipiv, FORWARD)
t.Logf(" blk: L*D*L.T %d pivots: ||A*B - L*D*L.T*B||_1: %e\n",
NumPivots(ipiv), NormP(B2.Minus(B1), NORM_ONE))
t.Logf("pivots: %v\n", ipiv)
}
开发者ID:hrautila,项目名称:matops,代码行数:59,代码来源:ldl_test.go
示例11: _TestRank
func _TestRank(t *testing.T) {
bM := M * 100
bN := N * 100
//bP := 5
A := matrix.FloatWithValue(bM, bN, 1.0)
A0 := matrix.FloatWithValue(bM, bN, 1.0)
X := matrix.FloatNormal(bM, 1)
Y := matrix.FloatNormal(bN, 1)
Ar := A.FloatArray()
Xr := X.FloatArray()
Yr := Y.FloatArray()
blas.GerFloat(X, Y, A0, 1.0)
DRankMV(Ar, Xr, Yr, 1.0, A.LeadingIndex(), 1, 1, 0, bN, 0, bM, 4, 4)
t.Logf("GER A0 == A1: %v\n", A0.AllClose(A))
}
开发者ID:hrautila,项目名称:matops,代码行数:19,代码来源:trank_test.go
示例12: _TestBKpivot1U
func _TestBKpivot1U(t *testing.T) {
Ldata := [][]float64{
[]float64{7.0, 6.0, 5.0, 4.0, 3.0, 2.0, 1.0},
[]float64{1.0, 6.0, 5.0, 4.0, 3.0, 2.0, 1.0},
[]float64{1.0, 2.0, 5.0, 4.0, 3.0, 5.0, 1.0},
[]float64{1.0, 2.0, 3.0, 4.0, 3.0, 2.0, 1.0},
[]float64{1.0, 5.0, 3.0, 4.0, 3.0, 2.0, 1.0},
[]float64{1.0, 2.0, 3.0, 4.0, 5.0, 2.0, 1.0},
[]float64{1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 1.0}}
Bdata := [][]float64{
[]float64{10.0, 20.0},
[]float64{10.0, 20.0},
[]float64{10.0, 20.0},
[]float64{10.0, 20.0},
[]float64{10.0, 20.0},
[]float64{10.0, 20.0},
[]float64{10.0, 20.0}}
A := matrix.FloatMatrixFromTable(Ldata, matrix.RowOrder)
X := matrix.FloatMatrixFromTable(Bdata, matrix.RowOrder)
N := A.Rows()
B := matrix.FloatZeros(N, 2)
MultSym(B, A, X, 1.0, 0.0, UPPER|LEFT)
//t.Logf("initial B:\n%v\n", B)
//N := 8
//A := matrix.FloatUniformSymmetric(N)
W := matrix.FloatWithValue(A.Rows(), 5, 0.0)
nb := 0
ipiv := make([]int, N, N)
L, _ := DecomposeBK(A.Copy(), W, ipiv, UPPER, nb)
t.Logf("ipiv: %v\n", ipiv)
t.Logf("L:\n%v\n", L)
nb = 4
//ipiv0 := make([]int, N, N)
//L0, _ := DecomposeBK(A.Copy(), W, ipiv0, LOWER, nb)
//t.Logf("ipiv: %v\n", ipiv0)
//t.Logf("L:\n%v\n", L0)
//B0 := B.Copy()
//SolveBK(B0, L0, ipiv0, LOWER)
//t.Logf("B0:\n%v\n", B0)
ipiv2 := make([]int32, N, N)
lapack.SytrfFloat(A, ipiv2, linalg.OptUpper)
t.Logf("ipiv2: %v\n", ipiv2)
t.Logf("lapack A:\n%v\n", A)
//lapack.Sytrs(A, B, ipiv2, linalg.OptLower)
//t.Logf("lapack B:\n%v\n", B)
//t.Logf("B == B0: %v\n", B.AllClose(B0))
}
开发者ID:hrautila,项目名称:matops,代码行数:53,代码来源:ldlbku_test.go
示例13: _TestMultSyr2Small
func _TestMultSyr2Small(t *testing.T) {
bN := 7
//A := matrix.FloatNormal(bN, bN)
//B := matrix.FloatNormal(bN, bP)
//A := matrix.FloatWithValue(bM, bP, 1.0)
X := matrix.FloatWithValue(bN, 1, 1.0)
Y := matrix.FloatWithValue(bN, 1, 1.0)
C0 := matrix.FloatZeros(bN, bN)
C1 := matrix.FloatZeros(bN, bN)
for i := 0; i < bN; i++ {
X.Add(1.0+float64(i), i)
Y.Add(2.0+float64(i), i)
}
t.Logf("X=\n%v\nY=\n%v\n", X, Y)
Xr := X.FloatArray()
Yr := Y.FloatArray()
C1r := C1.FloatArray()
blas.Syr2Float(X, Y, C0, 1.0, linalg.OptUpper)
DSymmRank2MV(C1r, Xr, Yr, 1.0, UPPER, C1.LeadingIndex(), 1, 1, 0, bN, 4)
ok := C0.AllClose(C1)
t.Logf("UPPER C0 == C1: %v\n", ok)
if !ok {
t.Logf("C1: C1 = A*X\n%v\n", C1)
t.Logf("blas: C0\n%v\n", C0)
}
blas.Syr2Float(X, Y, C0, 1.0, linalg.OptLower)
DSymmRank2MV(C1r, Xr, Yr, 1.0, LOWER, C1.LeadingIndex(), 1, 1, 0, bN, 4)
ok = C0.AllClose(C1)
t.Logf("LOWER C0 == C1: %v\n", ok)
if !ok {
t.Logf("blas: C0\n%v\n", C0)
t.Logf("C1: C1 = A*X\n%v\n", C1)
}
}
开发者ID:hrautila,项目名称:matops,代码行数:39,代码来源:trank_test.go
示例14: _TestBK2U
func _TestBK2U(t *testing.T) {
Bdata := [][]float64{
[]float64{10.0, 20.0},
[]float64{10.0, 20.0},
[]float64{10.0, 20.0},
[]float64{10.0, 20.0},
[]float64{10.0, 20.0},
[]float64{10.0, 20.0},
[]float64{10.0, 20.0}}
N := 7
A0 := matrix.FloatNormal(N, N)
A := matrix.FloatZeros(N, N)
// A is symmetric, posivite definite
Mult(A, A0, A0, 1.0, 1.0, TRANSB)
X := matrix.FloatMatrixFromTable(Bdata, matrix.RowOrder)
B := matrix.FloatZeros(N, 2)
MultSym(B, A, X, 1.0, 0.0, LOWER|LEFT)
t.Logf("initial B:\n%v\n", B)
nb := 0
W := matrix.FloatWithValue(A.Rows(), 5, 1.0)
A.SetAt(4, 1, A.GetAt(4, 1)+1.0)
A.SetAt(1, 4, A.GetAt(4, 1))
ipiv := make([]int, N, N)
L, _ := DecomposeBK(A.Copy(), W, ipiv, LOWER, nb)
t.Logf("ipiv: %v\n", ipiv)
t.Logf("L:\n%v\n", L)
ipiv0 := make([]int, N, N)
nb = 4
L0, _ := DecomposeBK(A.Copy(), W, ipiv0, LOWER, nb)
t.Logf("ipiv: %v\n", ipiv0)
t.Logf("L:\n%v\n", L0)
B0 := B.Copy()
SolveBK(B0, L0, ipiv0, LOWER)
t.Logf("B0:\n%v\n", B0)
ipiv2 := make([]int32, N, N)
lapack.Sytrf(A, ipiv2, linalg.OptLower)
t.Logf("ipiv2: %v\n", ipiv2)
t.Logf("lapack A:\n%v\n", A)
lapack.Sytrs(A, B, ipiv2, linalg.OptLower)
t.Logf("lapack B:\n%v\n", B)
t.Logf("B == B0: %v\n", B.AllClose(B0))
}
开发者ID:hrautila,项目名称:matops,代码行数:49,代码来源:ldlbku_test.go
示例15: _TestSolveRandom
func _TestSolveRandom(t *testing.T) {
bN := 22
nb := 4
L := matrix.FloatNormalSymmetric(bN, matrix.Lower)
U := L.Transpose()
X0 := matrix.FloatWithValue(L.Rows(), 1, 1.0)
t.Logf("-- BLOCKED SOLVE LOWER NON-UNIT ---\n")
solveMVTest(t, L, X0.Copy(), LOWER, bN, nb)
t.Logf("-- BLOCKED SOLVE LOWER UNIT ---\n")
solveMVTest(t, L, X0.Copy(), LOWER|UNIT, bN, nb)
t.Logf("-- BLOCKED SOLVE UPPER NON-UNIT ---\n")
solveMVTest(t, U, X0.Copy(), UPPER, bN, nb)
t.Logf("-- BLOCKED SOLVE UPPER UNIT ---\n")
solveMVTest(t, U, X0.Copy(), UPPER|UNIT, bN, nb)
}
开发者ID:hrautila,项目名称:matops,代码行数:16,代码来源:tsolve_test.go
示例16: main
func main() {
flag.Parse()
x := floorplan(matrix.FloatWithValue(5, 1, 100.0))
if x != nil {
W := x.GetIndex(0)
H := x.GetIndex(1)
xs := matrix.FloatVector(x.FloatArray()[2:7])
ys := matrix.FloatVector(x.FloatArray()[7:12])
ws := matrix.FloatVector(x.FloatArray()[12:17])
hs := matrix.FloatVector(x.FloatArray()[17:])
fmt.Printf("W = %.5f, H = %.5f\n", W, H)
fmt.Printf("x = \n%v\n", xs.ToString("%.5f"))
fmt.Printf("y = \n%v\n", ys.ToString("%.5f"))
fmt.Printf("w = \n%v\n", ws.ToString("%.5f"))
fmt.Printf("h = \n%v\n", hs.ToString("%.5f"))
check(x)
}
}
开发者ID:hrautila,项目名称:go.opt,代码行数:19,代码来源:testcpl.go
示例17: _TestBKSolve
func _TestBKSolve(t *testing.T) {
Ldata := [][]float64{
[]float64{1.0, 2.0, 3.0, 4.0},
[]float64{2.0, 2.0, 3.0, 4.0},
[]float64{3.0, 3.0, 3.0, 4.0},
[]float64{4.0, 4.0, 4.0, 4.0}}
Xdata := [][]float64{
[]float64{1.0, 2.0},
[]float64{1.0, 2.0},
[]float64{1.0, 2.0},
[]float64{1.0, 2.0}}
A := matrix.FloatMatrixFromTable(Ldata, matrix.RowOrder)
X := matrix.FloatMatrixFromTable(Xdata, matrix.RowOrder)
N := A.Rows()
B := matrix.FloatZeros(N, 2)
Mult(B, A, X, 1.0, 0.0, NOTRANS)
S := matrix.FloatZeros(N, 2)
MultSym(S, A, X, 1.0, 0.0, LOWER|LEFT)
t.Logf("B:\n%v\n", B)
t.Logf("S:\n%v\n", S)
//N := 8
//A := matrix.FloatUniformSymmetric(N)
nb := 0
W := matrix.FloatWithValue(A.Rows(), 5, 0.0)
ipiv := make([]int, N, N)
L, _ := DecomposeBK(A.Copy(), W, ipiv, LOWER, nb)
t.Logf("ipiv: %v\n", ipiv)
t.Logf("L:\n%v\n", L)
B0 := B.Copy()
SolveBK(B0, L, ipiv, LOWER)
t.Logf("B0:\n%v\n", B0)
ipiv2 := make([]int32, N, N)
lapack.Sytrf(A, ipiv2, linalg.OptLower)
t.Logf("ipiv2: %v\n", ipiv2)
t.Logf("lapack A:\n%v\n", A)
lapack.Sytrs(A, B, ipiv2, linalg.OptLower)
t.Logf("lapack B:\n%v\n", B)
}
开发者ID:hrautila,项目名称:matops,代码行数:42,代码来源:ldlbkl_test.go
示例18: _TestMultTransASmall
func _TestMultTransASmall(t *testing.T) {
bM := 7
bN := 7
bP := 7
D := matrix.FloatNormal(bM, bP)
E := matrix.FloatNormal(bP, bN)
C0 := matrix.FloatWithValue(bM, bN, 0.0)
C1 := C0.Copy()
Dt := D.Transpose()
Dr := Dt.FloatArray()
Er := E.FloatArray()
C1r := C1.FloatArray()
blas.GemmFloat(Dt, E, C0, 1.0, 1.0, linalg.OptTransA)
t.Logf("blas: C=D*E\n%v\n", C0)
DMult(C1r, Dr, Er, 1.0, 1.0, TRANSA, bM, bM, bP, bP, 0, bN, 0, bM, 4, 4, 4)
t.Logf("C0 == C1: %v\n", C0.AllClose(C1))
t.Logf("C1: C1=D*E\n%v\n", C1)
}
开发者ID:hrautila,项目名称:matops,代码行数:20,代码来源:tmult_test.go
示例19: _TestMultSmall
func _TestMultSmall(t *testing.T) {
bM := 6
bN := 6
bP := 6
D := matrix.FloatNormal(bM, bP)
E := matrix.FloatNormal(bP, bN)
C0 := matrix.FloatWithValue(bM, bN, 1.0)
C1 := C0.Copy()
Dr := D.FloatArray()
Er := E.FloatArray()
C1r := C1.FloatArray()
blas.GemmFloat(D, E, C0, 1.0, 1.0)
t.Logf("blas: C=D*E\n%v\n", C0)
DMult(C1r, Dr, Er, 1.0, 1.0, NOTRANS, bM, bM, bP, bP, 0, bN, 0, bM, 4, 4, 4)
t.Logf("C0 == C1: %v\n", C0.AllClose(C1))
t.Logf("C1: C1=D*E\n%v\n", C1)
}
开发者ID:hrautila,项目名称:matops,代码行数:20,代码来源:tmult_test.go
示例20: main
func main() {
Sdata := [][]float64{
[]float64{4e-2, 6e-3, -4e-3, 0.0},
[]float64{6e-3, 1e-2, 0.0, 0.0},
[]float64{-4e-3, 0.0, 2.5e-3, 0.0},
[]float64{0.0, 0.0, 0.0, 0.0}}
pbar := matrix.FloatVector([]float64{.12, .10, .07, .03})
S := matrix.FloatMatrixFromTable(Sdata)
n := pbar.Rows()
G := matrix.FloatDiagonal(n, -1.0)
h := matrix.FloatZeros(n, 1)
A := matrix.FloatWithValue(1, n, 1.0)
b := matrix.FloatNew(1, 1, []float64{1.0})
var solopts cvx.SolverOptions
solopts.MaxIter = 30
solopts.ShowProgress = true
mu := 1.0
Smu := matrix.Scale(S, mu)
pbarNeg := matrix.Scale(pbar, -1.0)
fmt.Printf("Smu=\n%v\n", Smu.String())
fmt.Printf("-pbar=\n%v\n", pbarNeg.String())
sol, err := cvx.Qp(Smu, pbarNeg, G, h, A, b, &solopts, nil)
fmt.Printf("status: %v\n", err)
if sol != nil && sol.Status == cvx.Optimal {
x := sol.Result.At("x")[0]
ret := blas.DotFloat(x, pbar)
risk := math.Sqrt(blas.DotFloat(x, S.Times(x)))
fmt.Printf("ret=%.3f, risk=%.3f\n", ret, risk)
fmt.Printf("x=\n%v\n", x)
}
}
开发者ID:hrautila,项目名称:go.opt,代码行数:37,代码来源:testqp.go
注:本文中的github.com/hrautila/matrix.FloatWithValue函数示例整理自Github/MSDocs等源码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。 |
请发表评论