Hulk did something

[bytom/vapor.git] / vendor / gonum.org / v1 / gonum / lapack / gonum / dtrtrs.go

diff --git a/vendor/gonum.org/v1/gonum/lapack/gonum/dtrtrs.go b/vendor/gonum.org/v1/gonum/lapack/gonum/dtrtrs.go
new file mode 100644 (file)
index 0000000..e1782d2
--- /dev/null
+++ b/vendor/gonum.org/v1/gonum/lapack/gonum/dtrtrs.go
@@ -0,0 +1,30 @@
+// Copyright ©2015 The Gonum Authors. All rights reserved.
+// Use of this source code is governed by a BSD-style
+// license that can be found in the LICENSE file.
+
+package gonum
+
+import (
+       "gonum.org/v1/gonum/blas"
+       "gonum.org/v1/gonum/blas/blas64"
+)
+
+// Dtrtrs solves a triangular system of the form A * X = B or A^T * X = B. Dtrtrs
+// returns whether the solve completed successfully. If A is singular, no solve is performed.
+func (impl Implementation) Dtrtrs(uplo blas.Uplo, trans blas.Transpose, diag blas.Diag, n, nrhs int, a []float64, lda int, b []float64, ldb int) (ok bool) {
+       nounit := diag == blas.NonUnit
+       if n == 0 {
+               return false
+       }
+       // Check for singularity.
+       if nounit {
+               for i := 0; i < n; i++ {
+                       if a[i*lda+i] == 0 {
+                               return false
+                       }
+               }
+       }
+       bi := blas64.Implementation()
+       bi.Dtrsm(blas.Left, uplo, trans, diag, n, nrhs, 1, a, lda, b, ldb)
+       return true
+}