Home
last modified time | relevance | path

Searched refs:adjoint (Results 1 – 25 of 255) sorted by relevance

1234567891011

/external/eigen/test/
Dproduct_syrk.cpp39 … ((s1 * rhs2 * rhs2.adjoint()).eval().template triangularView<Lower>().toDenseMatrix())); in syrk()
41 …VERIFY_IS_APPROX(((m2.template triangularView<Lower>() += s1 * rhs2 * rhs22.adjoint()).nestedExpr… in syrk()
42 … ((s1 * rhs2 * rhs22.adjoint()).eval().template triangularView<Lower>().toDenseMatrix())); in syrk()
47 … (s1 * rhs2 * rhs2.adjoint()).eval().template triangularView<Upper>().toDenseMatrix()); in syrk()
49 …VERIFY_IS_APPROX((m2.template triangularView<Upper>() += s1 * rhs22 * rhs2.adjoint()).nestedExpres… in syrk()
50 … (s1 * rhs22 * rhs2.adjoint()).eval().template triangularView<Upper>().toDenseMatrix()); in syrk()
54 VERIFY_IS_APPROX(m2.template selfadjointView<Lower>().rankUpdate(rhs1.adjoint(),s1)._expression(), in syrk()
55 … (s1 * rhs1.adjoint() * rhs1).eval().template triangularView<Lower>().toDenseMatrix()); in syrk()
57 …VERIFY_IS_APPROX((m2.template triangularView<Lower>() += s1 * rhs11.adjoint() * rhs1).nestedExpres… in syrk()
58 … (s1 * rhs11.adjoint() * rhs1).eval().template triangularView<Lower>().toDenseMatrix()); in syrk()
[all …]
Dmixingtypes.cpp160 VERIFY_IS_APPROX(sd*md.adjoint()*mcd, (sd*md).template cast<CD>().eval().adjoint()*mcd); in mixingtypes()
161 VERIFY_IS_APPROX(sd*mcd.adjoint()*md, sd*mcd.adjoint()*md.template cast<CD>()); in mixingtypes()
162 …VERIFY_IS_APPROX(sd*md.adjoint()*mcd.adjoint(), (sd*md).template cast<CD>().eval().adjoint()*mcd.a… in mixingtypes()
163 …VERIFY_IS_APPROX(sd*mcd.adjoint()*md.adjoint(), sd*mcd.adjoint()*md.template cast<CD>().adjoint()); in mixingtypes()
164 VERIFY_IS_APPROX(sd*md*mcd.adjoint(), (sd*md).template cast<CD>().eval()*mcd.adjoint()); in mixingtypes()
165 VERIFY_IS_APPROX(sd*mcd*md.adjoint(), sd*mcd*md.template cast<CD>().adjoint()); in mixingtypes()
167 VERIFY_IS_APPROX(sf*mf.adjoint()*mcf, (sf*mf).template cast<CF>().eval().adjoint()*mcf); in mixingtypes()
168 VERIFY_IS_APPROX(sf*mcf.adjoint()*mf, sf*mcf.adjoint()*mf.template cast<CF>()); in mixingtypes()
169 …VERIFY_IS_APPROX(sf*mf.adjoint()*mcf.adjoint(), (sf*mf).template cast<CF>().eval().adjoint()*mcf.a… in mixingtypes()
170 …VERIFY_IS_APPROX(sf*mcf.adjoint()*mf.adjoint(), sf*mcf.adjoint()*mf.template cast<CF>().adjoint()); in mixingtypes()
[all …]
Dadjoint.cpp23 VERIFY(test_isApproxWithRef(v1.dot(square * v2), (square.adjoint() * v1).dot(v2), 0)); in run()
58 …)((std::max)(v1.norm(),v2.norm()),(std::max)((square * v2).norm(),(square.adjoint() * v1).norm())); in run()
59 …VERIFY(internal::isMuchSmallerThan(abs(v1.dot(square * v2) - (square.adjoint() * v1).dot(v2)), ref… in run()
67 template<typename MatrixType> void adjoint(const MatrixType& m) in adjoint() function
96 VERIFY_IS_APPROX(m1.transpose().conjugate().adjoint(), m1); in adjoint()
97 VERIFY_IS_APPROX(m1.adjoint().conjugate().transpose(), m1); in adjoint()
100 VERIFY_IS_APPROX((m1.adjoint() * m2).adjoint(), m2.adjoint() * m1); in adjoint()
101 VERIFY_IS_APPROX((s1 * m1).adjoint(), numext::conj(s1) * m1.adjoint()); in adjoint()
115 VERIFY_IS_APPROX(m1.adjoint()(c,r), numext::conj(m1(r,c))); in adjoint()
138 VERIFY_IS_APPROX(m3,m1.adjoint()); in adjoint()
[all …]
Dproduct_notemporary.cpp45 VERIFY_EVALUATION_COUNT( m3 = (m1 * m2.adjoint()), 1); in product_notemporary()
46 VERIFY_EVALUATION_COUNT( m3 = (m1 * m2.adjoint()).transpose(), 1); in product_notemporary()
47 VERIFY_EVALUATION_COUNT( m3.noalias() = m1 * m2.adjoint(), 0); in product_notemporary()
53 VERIFY_EVALUATION_COUNT( m3 = m3 + (m1 * m2.adjoint()), 1); in product_notemporary()
54 VERIFY_EVALUATION_COUNT( m3 = m3 - (m1 * m2.adjoint()), 1); in product_notemporary()
56 VERIFY_EVALUATION_COUNT( m3 = m3 + (m1 * m2.adjoint()).transpose(), 1); in product_notemporary()
64 VERIFY_EVALUATION_COUNT( m3.noalias() = s1 * m1 * s2 * m2.adjoint(), 0); in product_notemporary()
65 VERIFY_EVALUATION_COUNT( m3.noalias() = s1 * m1 * s2 * (m1*s3+m2*s2).adjoint(), 1); in product_notemporary()
66 VERIFY_EVALUATION_COUNT( m3.noalias() = (s1 * m1).adjoint() * s2 * m2, 0); in product_notemporary()
67 VERIFY_EVALUATION_COUNT( m3.noalias() += s1 * (-m1*s3).adjoint() * (s2 * m2 * s3), 0); in product_notemporary()
[all …]
Dproduct_extra.cpp41 VERIFY_IS_APPROX(m3.noalias() = m1 * m2.adjoint(), m1 * m2.adjoint().eval()); in product_extra()
42 …VERIFY_IS_APPROX(m3.noalias() = m1.adjoint() * square.adjoint(), m1.adjoint().eval() * square.ad… in product_extra()
43 VERIFY_IS_APPROX(m3.noalias() = m1.adjoint() * m2, m1.adjoint().eval() * m2); in product_extra()
44 …VERIFY_IS_APPROX(m3.noalias() = (s1 * m1.adjoint()) * m2, (s1 * m1.adjoint()).eval() * m2… in product_extra()
45 …VERIFY_IS_APPROX(m3.noalias() = ((s1 * m1).adjoint()) * m2, (numext::conj(s1) * m1.adjoint(… in product_extra()
46 …VERIFY_IS_APPROX(m3.noalias() = (- m1.adjoint() * s1) * (s3 * m2), (- m1.adjoint() * s1).eval() *… in product_extra()
47 …VERIFY_IS_APPROX(m3.noalias() = (s2 * m1.adjoint() * s1) * m2, (s2 * m1.adjoint() * s1).eval(… in product_extra()
48 …VERIFY_IS_APPROX(m3.noalias() = (-m1*s2) * s1*m2.adjoint(), (-m1*s2).eval() * (s1*m2.adjoin… in product_extra()
51 …VERIFY_IS_APPROX( m1.adjoint() * (s1*m2).conjugate(), (m1.adjoint()).eval() * ((s1*m2).conjugate()… in product_extra()
63 VERIFY_IS_APPROX((s1 * vc2.transpose()) * (-m1.adjoint() * s2), in product_extra()
[all …]
Dproduct_selfadjoint.cpp38 m1 = (m1.adjoint() + m1).eval(); in product_selfadjoint()
43 …VERIFY_IS_APPROX(m2, (m1 + v1 * v2.adjoint()+ v2 * v1.adjoint()).template triangularView<Lower>().… in product_selfadjoint()
47 …VERIFY_IS_APPROX(m2, (m1 + (s3*(-v1)*(s2*v2).adjoint()+numext::conj(s3)*(s2*v2)*(-v1).adjoint())).… in product_selfadjoint()
50 m2.template selfadjointView<Upper>().rankUpdate(-s2*r1.adjoint(),r2.adjoint()*s3,s1); in product_selfadjoint()
51 …OX(m2, (m1 + s1*(-s2*r1.adjoint())*(r2.adjoint()*s3).adjoint() + numext::conj(s1)*(r2.adjoint()*s3… in product_selfadjoint()
58 …ws-1,cols-1) += v1.tail(rows-1) * v2.head(cols-1).adjoint()+ v2.head(cols-1) * v1.tail(rows-1).adj… in product_selfadjoint()
Dnomalloc.cpp48 m2.col(0).noalias() -= m1.adjoint() * m1.col(0); in nomalloc()
49 m2.col(0).noalias() -= m1 * m1.row(0).adjoint(); in nomalloc()
50 m2.col(0).noalias() -= m1.adjoint() * m1.row(0).adjoint(); in nomalloc()
53 m2.row(0).noalias() -= m1.row(0) * m1.adjoint(); in nomalloc()
54 m2.row(0).noalias() -= m1.col(0).adjoint() * m1; in nomalloc()
55 m2.row(0).noalias() -= m1.col(0).adjoint() * m1.adjoint(); in nomalloc()
59 m2.col(0).noalias() -= m1.adjoint().template triangularView<Upper>() * m1.col(0); in nomalloc()
60 m2.col(0).noalias() -= m1.template triangularView<Upper>() * m1.row(0).adjoint(); in nomalloc()
61 m2.col(0).noalias() -= m1.adjoint().template triangularView<Upper>() * m1.row(0).adjoint(); in nomalloc()
64 m2.row(0).noalias() -= m1.row(0) * m1.adjoint().template triangularView<Upper>(); in nomalloc()
[all …]
Dproduct_trmm.cpp54 …VERIFY_IS_APPROX( ge_xs.noalias() = (s1*mat.adjoint()).template triangularView<Mode>() * (s2*ge_le… in trmm()
55 …VERIFY_IS_APPROX( ge_sx.noalias() = ge_right.transpose() * mat.adjoint().template triangularView<M… in trmm()
57 ….noalias() = (s1*mat.adjoint()).template triangularView<Mode>() * (s2*ge_left.adjoint()), s1*triTr… in trmm()
58 …_IS_APPROX( ge_sx.noalias() = ge_right.adjoint() * mat.adjoint().template triangularView<Mode>(), … in trmm()
61 …gate() * (s2*ge_left.adjoint())).eval(), ge_xs.noalias() += (s1*mat.adjoint()).template triangular… in trmm()
64 …ge_sx_save - (ge_right.adjoint() * (-s1 * triTr).conjugate()).eval(), ge_sx.noalias() -= (ge_right… in trmm()
66 … ge_xs = (s1*mat).adjoint().template triangularView<Mode>() * ge_left.adjoint(), numext::conj(s1) … in trmm()
Deigensolver_generalized_real.cpp33 MatrixType spdA = a.adjoint() * a + a1.adjoint() * a1; in generalized_eigensolver_real()
34 MatrixType spdB = b.adjoint() * b + b1.adjoint() * b1; in generalized_eigensolver_real()
75 GeneralizedSelfAdjointEigenSolver<MatrixType> eig1(a.adjoint() * a,b.adjoint() * b); in generalized_eigensolver_real()
76 eig1.compute(a.adjoint() * a,b.adjoint() * b); in generalized_eigensolver_real()
77 GeneralizedEigenSolver<MatrixType> eig2(a.adjoint() * a,b.adjoint() * b); in generalized_eigensolver_real()
78 eig2.compute(a.adjoint() * a,b.adjoint() * b); in generalized_eigensolver_real()
Dproduct_mmtr.cpp48 CHECK_MMTR(matc, Lower, = s*soc*sor.adjoint()); in mmtr()
49 CHECK_MMTR(matc, Upper, = s*(soc*soc.adjoint())); in mmtr()
50 CHECK_MMTR(matr, Lower, = s*soc*soc.adjoint()); in mmtr()
51 CHECK_MMTR(matr, Upper, = soc*(s*sor.adjoint())); in mmtr()
53 CHECK_MMTR(matc, Lower, += s*soc*soc.adjoint()); in mmtr()
55 CHECK_MMTR(matr, Lower, += s*sor*soc.adjoint()); in mmtr()
56 CHECK_MMTR(matr, Upper, += soc*(s*soc.adjoint())); in mmtr()
58 CHECK_MMTR(matc, Lower, -= s*soc*soc.adjoint()); in mmtr()
60 CHECK_MMTR(matr, Lower, -= s*soc*soc.adjoint()); in mmtr()
61 CHECK_MMTR(matr, Upper, -= soc*(s*soc.adjoint())); in mmtr()
[all …]
Dproduct_symm.cpp27 m1 = (m1+m1.adjoint()).eval(); in symm()
54 VERIFY_IS_APPROX(rhs12 = (s1*m2).adjoint().template selfadjointView<Upper>() * (s2*rhs1), in symm()
55 rhs13 = (s1*m1).adjoint() * (s2*rhs1)); in symm()
57 VERIFY_IS_APPROX(rhs12 = (s1*m2).template selfadjointView<Lower>().adjoint() * (s2*rhs1), in symm()
58 rhs13 = (s1*m1).adjoint() * (s2*rhs1)); in symm()
67 VERIFY_IS_APPROX(rhs12 = (s1*m2).template selfadjointView<Lower>() * (s2*rhs2.adjoint()), in symm()
68 rhs13 = (s1*m1) * (s2*rhs2.adjoint())); in symm()
71 VERIFY_IS_APPROX(rhs12 = (s1*m2).template selfadjointView<Upper>() * (s2*rhs2.adjoint()), in symm()
72 rhs13 = (s1*m1) * (s2*rhs2.adjoint())); in symm()
75 …VERIFY_IS_APPROX(rhs12 = (s1*m2.adjoint()).template selfadjointView<Lower>() * (s2*rhs2.adjoint()), in symm()
[all …]
Dproduct_trmv.cpp60 …VERIFY((m3.adjoint() * v1).isApprox(m1.adjoint().template triangularView<Eigen::Lower>() * v1, lar… in trmv()
62 …VERIFY((m3.adjoint() * (s1*v1.conjugate())).isApprox(m1.adjoint().template triangularView<Eigen::U… in trmv()
68 …VERIFY((v1.adjoint() * m3).isApprox(v1.adjoint() * m1.template triangularView<Eigen::Lower>(), lar… in trmv()
69 …VERIFY((v1.adjoint() * m3.adjoint()).isApprox(v1.adjoint() * m1.template triangularView<Eigen::Low… in trmv()
/external/tensorflow/tensorflow/python/ops/linalg/
Dlinear_operator_block_lower_triangular.py394 def matmul(self, x, adjoint=False, adjoint_arg=False, name="matmul"): argument
426 left_operator = self.adjoint() if adjoint else self
427 right_operator = x.adjoint() if adjoint_arg else x
441 block_dimensions = (self._block_range_dimensions() if adjoint
453 op_dimension = (self.range_dimension if adjoint
456 return self._matmul(x, adjoint=adjoint, adjoint_arg=adjoint_arg)
458 def _matmul(self, x, adjoint=False, adjoint_arg=False): argument
460 block_dimensions = (self._block_range_dimensions() if adjoint
476 if adjoint:
481 split_x[index], adjoint=adjoint, adjoint_arg=adjoint_arg)
[all …]
Dlinear_operator.py601 linalg.adjoint(dense),
628 def _matmul(self, x, adjoint=False, adjoint_arg=False): argument
631 def matmul(self, x, adjoint=False, adjoint_arg=False, name="matmul"): argument
661 left_operator = self.adjoint() if adjoint else self
662 right_operator = x.adjoint() if adjoint_arg else x
678 self_dim = -2 if adjoint else -1
684 return self._matmul(x, adjoint=adjoint, adjoint_arg=adjoint_arg)
689 def _matvec(self, x, adjoint=False): argument
691 y_mat = self.matmul(x_mat, adjoint=adjoint)
694 def matvec(self, x, adjoint=False, name="matvec"): argument
[all …]
Dlinear_operator_block_diag.py290 def matmul(self, x, adjoint=False, adjoint_arg=False, name="matmul"): argument
322 left_operator = self.adjoint() if adjoint else self
323 right_operator = x.adjoint() if adjoint_arg else x
337 block_dimensions = (self._block_range_dimensions() if adjoint
349 op_dimension = (self.range_dimension if adjoint
352 return self._matmul(x, adjoint=adjoint, adjoint_arg=adjoint_arg)
354 def _matmul(self, x, adjoint=False, adjoint_arg=False): argument
356 block_dimensions = (self._block_range_dimensions() if adjoint
373 split_x[index], adjoint=adjoint, adjoint_arg=adjoint_arg)]
382 def matvec(self, x, adjoint=False, name="matvec"): argument
[all …]
Dlinear_operator_adjoint.py190 def _matmul(self, x, adjoint=False, adjoint_arg=False): argument
192 x, adjoint=(not adjoint), adjoint_arg=adjoint_arg)
194 def _matvec(self, x, adjoint=False): argument
195 return self.operator.matvec(x, adjoint=(not adjoint))
210 def _solve(self, rhs, adjoint=False, adjoint_arg=False): argument
212 rhs, adjoint=(not adjoint), adjoint_arg=adjoint_arg)
214 def _solvevec(self, rhs, adjoint=False): argument
215 return self.operator.solvevec(rhs, adjoint=(not adjoint))
220 return linalg.adjoint(self.operator.to_dense())
Dlinear_operator_test_util.py156 def make_rhs(self, operator, adjoint, with_batch=True): argument
173 def make_x(self, operator, adjoint, with_batch=True): argument
229 tape.gradient(operator.adjoint().to_dense(), operator.variables))
326 adjoint, argument
341 operator, adjoint=adjoint, with_batch=with_batch)
345 linalg.adjoint(x),
346 adjoint=adjoint,
349 op_matmul = operator.matmul(x, adjoint=adjoint)
350 mat_matmul = math_ops.matmul(mat, x, adjoint_a=adjoint)
364 split_x = [linalg.adjoint(y) for y in split_x]
[all …]
Dlinear_operator_composition.py243 def _matmul(self, x, adjoint=False, adjoint_arg=False): argument
247 if adjoint:
253 x, adjoint=adjoint, adjoint_arg=adjoint_arg)
255 result = operator.matmul(result, adjoint=adjoint)
270 def _solve(self, rhs, adjoint=False, adjoint_arg=False): argument
278 if adjoint:
284 rhs, adjoint=adjoint, adjoint_arg=adjoint_arg)
286 solution = operator.solve(solution, adjoint=adjoint)
Dlinear_operator_low_rank_update.py364 def _matmul(self, x, adjoint=False, adjoint_arg=False): argument
370 leading_term = l.matmul(x, adjoint=adjoint, adjoint_arg=adjoint_arg)
372 if adjoint:
374 d_uh_x = d.matmul(uh_x, adjoint=adjoint)
379 d_vh_x = d.matmul(vh_x, adjoint=adjoint)
417 def _solve(self, rhs, adjoint=False, adjoint_arg=False): argument
432 if adjoint:
440 linv_rhs = l.solve(rhs, adjoint=adjoint, adjoint_arg=adjoint_arg)
449 self._make_capacitance(), vh_linv_rhs, adjoint=adjoint)
453 linv_u_capinv_vh_linv_rhs = l.solve(u_capinv_vh_linv_rhs, adjoint=adjoint)
Dlinear_operator_kronecker.py288 def _matmul(self, x, adjoint=False, adjoint_arg=False): argument
327 x = linalg.adjoint(x)
353 if adjoint:
365 output = operator.matmul(output, adjoint=adjoint, adjoint_arg=False)
378 output = self.operators[-1].matvec(output, adjoint=adjoint)
388 if adjoint:
433 def _solve(self, rhs, adjoint=False, adjoint_arg=False): argument
443 rhs = linalg.adjoint(rhs)
470 if adjoint:
482 output = operator.solve(output, adjoint=adjoint, adjoint_arg=False)
[all …]
/external/tensorflow/tensorflow/compiler/tests/
Dmatrix_solve_op_test.py32 def _verifySolve(self, x, y, adjoint): argument
37 np_ans = np.linalg.solve(np.swapaxes(a, -2, -1) if adjoint else a, b)
40 tf_ans = linalg_ops.matrix_solve(a, b, adjoint=adjoint)
56 def testSolve(self, n, nrhs, batch_dims, rhs_batch_dims, adjoint): argument
59 self._verifySolve(matrix, rhs, adjoint=adjoint)
65 def testConcurrent(self, adjoint): argument
72 s1 = linalg_ops.matrix_solve(lhs1, rhs1, adjoint=adjoint)
73 s2 = linalg_ops.matrix_solve(lhs2, rhs2, adjoint=adjoint)
/external/tensorflow/tensorflow/python/kernel_tests/
Dmatrix_solve_op_test.py47 for adjoint in False, True:
51 a_np = np.transpose(a) if adjoint else a
55 a_np = np.conj(np.transpose(a)) if adjoint else a
65 tf_ans = linalg_ops.matrix_solve(a_ph, b_ph, adjoint=adjoint)
69 tf_ans = linalg_ops.matrix_solve(a, b, adjoint=adjoint)
137 s1 = linalg_ops.matrix_solve(lhs1, rhs1, adjoint=adjoint_)
138 s2 = linalg_ops.matrix_solve(lhs2, rhs2, adjoint=adjoint_)
177 for adjoint in False, True:
185 x = linalg_ops.matrix_solve(matrix, rhs, adjoint=adjoint)
196 adjoint=adjoint))
[all …]
Dmatrix_inverse_op_test.py38 for adjoint in False, True:
42 inv = linalg_ops.matrix_inverse(y, adjoint=adjoint)
43 tf_ans = test_util.matmul_without_tf32(inv, y, adjoint_b=adjoint)
147 inv1 = linalg_ops.matrix_inverse(matrix1, adjoint=adjoint_)
148 inv2 = linalg_ops.matrix_inverse(matrix2, adjoint=adjoint_)
181 for adjoint in False, True:
187 inv = linalg_ops.matrix_inverse(matrix, adjoint=adjoint)
194 shape=shape, adjoint=adjoint))
201 inv = linalg_ops.matrix_inverse(matrix, adjoint=adjoint)
208 shape=shape, adjoint=adjoint))
/external/tensorflow/tensorflow/python/kernel_tests/linalg/
Dlinear_operator_test.py105 def _matmul(self, x, adjoint=False, adjoint_arg=False): argument
108 self._matrix, x, adjoint_a=adjoint, adjoint_b=adjoint_arg)
110 def _solve(self, rhs, adjoint=False, adjoint_arg=False): argument
113 return linalg_ops.matrix_solve(self._matrix, rhs, adjoint=adjoint)
366 adjoint = linalg.adjoint(operator)
367 self.assertIsInstance(adjoint, linalg.LinearOperator)
381 for adjoint in [False, True]:
384 rhs, adjoint=adjoint, adjoint_arg=adjoint_arg)
386 operator, rhs, adjoint_a=adjoint, adjoint_b=adjoint_arg)
390 op_val = operator.solve(rhs, adjoint=adjoint)
[all …]
/external/tensorflow/tensorflow/core/kernels/linalg/
Dcholesky_grad.cc108 C_bar = D.adjoint() in ComputeMatrix()
110 .solve(C_bar.adjoint()) in ComputeMatrix()
111 .adjoint(); in ComputeMatrix()
112 D_bar -= (C_bar.adjoint() * C).template triangularView<Eigen::Lower>(); in ComputeMatrix()
114 R_bar -= C_bar.adjoint() * B; in ComputeMatrix()
116 R_bar -= (D_bar + D_bar.adjoint()) * R; in ComputeMatrix()
150 d_bar -= (c.adjoint() * c_bar) / d; in CholeskyGradUnblocked()
152 r_bar -= d_stack_c_bar.adjoint() * r_stack_B; in CholeskyGradUnblocked()

1234567891011