| 218 | } |
| 219 | |
| 220 | template<typename ArrayType> void array_complex(const ArrayType& m) |
| 221 | { |
| 222 | typedef typename ArrayType::Index Index; |
| 223 | |
| 224 | Index rows = m.rows(); |
| 225 | Index cols = m.cols(); |
| 226 | |
| 227 | ArrayType m1 = ArrayType::Random(rows, cols), |
| 228 | m2(rows, cols); |
| 229 | |
| 230 | for (Index i = 0; i < m.rows(); ++i) |
| 231 | for (Index j = 0; j < m.cols(); ++j) |
| 232 | m2(i,j) = sqrt(m1(i,j)); |
| 233 | |
| 234 | VERIFY_IS_APPROX(m1.sqrt(), m2); |
| 235 | VERIFY_IS_APPROX(m1.sqrt(), Eigen::sqrt(m1)); |
| 236 | } |
| 237 | |
| 238 | template<typename ArrayType> void min_max(const ArrayType& m) |
| 239 | { |
no test coverage detected