3 from pybind11_tests
import ConstructorStats
5 np = pytest.importorskip(
"numpy")
6 m = pytest.importorskip(
"pybind11_tests.eigen")
11 [0.0, 3, 0, 0, 0, 11],
12 [22, 0, 0, 0, 17, 11],
21 np.testing.assert_array_equal(mat, ref)
47 ref2 = np.array([[0.0, 1, 2, 3], [4, 5, 6, 7], [8, 9, 10, 11], [12, 13, 14, 15]])
48 np.testing.assert_array_equal(m.partial_copy_four_rm_r(ref2), ref2)
49 np.testing.assert_array_equal(m.partial_copy_four_rm_c(ref2), ref2)
50 np.testing.assert_array_equal(m.partial_copy_four_rm_r(ref2[:, 1]), ref2[:, [1]])
51 np.testing.assert_array_equal(m.partial_copy_four_rm_c(ref2[0, :]), ref2[[0], :])
52 np.testing.assert_array_equal(
53 m.partial_copy_four_rm_r(ref2[:, (0, 2)]), ref2[:, (0, 2)]
55 np.testing.assert_array_equal(
56 m.partial_copy_four_rm_c(ref2[(3, 1, 2), :]), ref2[(3, 1, 2), :]
59 np.testing.assert_array_equal(m.partial_copy_four_cm_r(ref2), ref2)
60 np.testing.assert_array_equal(m.partial_copy_four_cm_c(ref2), ref2)
61 np.testing.assert_array_equal(m.partial_copy_four_cm_r(ref2[:, 1]), ref2[:, [1]])
62 np.testing.assert_array_equal(m.partial_copy_four_cm_c(ref2[0, :]), ref2[[0], :])
63 np.testing.assert_array_equal(
64 m.partial_copy_four_cm_r(ref2[:, (0, 2)]), ref2[:, (0, 2)]
66 np.testing.assert_array_equal(
67 m.partial_copy_four_cm_c(ref2[(3, 1, 2), :]), ref2[(3, 1, 2), :]
72 m.partial_copy_four_rm_r,
73 m.partial_copy_four_rm_c,
74 m.partial_copy_four_cm_r,
75 m.partial_copy_four_cm_c,
77 matrix_with_wrong_shape = [[1, 2], [3, 4]]
79 with pytest.raises(TypeError)
as excinfo:
80 f(matrix_with_wrong_shape)
81 assert "incompatible function arguments" in str(excinfo.value)
85 zr = np.arange(30, dtype=
"float32").
reshape(5, 6)
86 zc = zr.reshape(6, 5).transpose()
92 with pytest.raises(TypeError)
as excinfo:
95 "(arg0: numpy.ndarray[numpy.float32[5, 6],"
96 " flags.writeable, flags.c_contiguous]) -> None" in str(excinfo.value)
98 with pytest.raises(TypeError)
as excinfo:
101 "(arg0: numpy.ndarray[numpy.float32[5, 6],"
102 " flags.writeable, flags.f_contiguous]) -> None" in str(excinfo.value)
104 with pytest.raises(TypeError)
as excinfo:
105 m.fixed_mutator_a(np.array([[1, 2], [3, 4]], dtype=
"float32"))
106 assert "(arg0: numpy.ndarray[numpy.float32[5, 6], flags.writeable]) -> None" in str(
109 zr.flags.writeable =
False
110 with pytest.raises(TypeError):
111 m.fixed_mutator_r(zr)
112 with pytest.raises(TypeError):
113 m.fixed_mutator_a(zr)
117 assert m.cpp_copy(m.fixed_r()) == 22.0
118 assert m.cpp_copy(m.fixed_c()) == 22.0
119 z = np.array([[5.0, 6], [7, 8]])
120 assert m.cpp_copy(z) == 7.0
121 assert m.cpp_copy(m.get_cm_ref()) == 21.0
122 assert m.cpp_copy(m.get_rm_ref()) == 21.0
123 assert m.cpp_ref_c(m.get_cm_ref()) == 21.0
124 assert m.cpp_ref_r(m.get_rm_ref()) == 21.0
125 with pytest.raises(RuntimeError)
as excinfo:
127 m.cpp_ref_any(m.fixed_c())
128 assert "Unable to cast Python instance" in str(excinfo.value)
129 with pytest.raises(RuntimeError)
as excinfo:
131 m.cpp_ref_any(m.fixed_r())
132 assert "Unable to cast Python instance" in str(excinfo.value)
133 assert m.cpp_ref_any(m.ReturnTester.create()) == 1.0
135 assert m.cpp_ref_any(m.get_cm_ref()) == 21.0
136 assert m.cpp_ref_any(m.get_cm_ref()) == 21.0
140 z = np.full((5, 6), 42.0)
141 z.flags.writeable =
False
142 np.testing.assert_array_equal(z, m.fixed_copy_r(z))
143 np.testing.assert_array_equal(m.fixed_r_const(), m.fixed_r())
144 assert not m.fixed_r_const().flags.writeable
145 np.testing.assert_array_equal(m.fixed_copy_r(m.fixed_r_const()), m.fixed_r_const())
149 counting_mat = np.arange(9.0, dtype=np.float32).
reshape((3, 3))
150 second_row = counting_mat[1, :]
151 second_col = counting_mat[:, 1]
152 np.testing.assert_array_equal(m.double_row(second_row), 2.0 * second_row)
153 np.testing.assert_array_equal(m.double_col(second_row), 2.0 * second_row)
154 np.testing.assert_array_equal(m.double_complex(second_row), 2.0 * second_row)
155 np.testing.assert_array_equal(m.double_row(second_col), 2.0 * second_col)
156 np.testing.assert_array_equal(m.double_col(second_col), 2.0 * second_col)
157 np.testing.assert_array_equal(m.double_complex(second_col), 2.0 * second_col)
159 counting_3d = np.arange(27.0, dtype=np.float32).
reshape((3, 3, 3))
160 slices = [counting_3d[0, :, :], counting_3d[:, 0, :], counting_3d[:, :, 0]]
161 for ref_mat
in slices:
162 np.testing.assert_array_equal(m.double_mat_cm(ref_mat), 2.0 * ref_mat)
163 np.testing.assert_array_equal(m.double_mat_rm(ref_mat), 2.0 * ref_mat)
166 m.double_threer(second_row)
167 m.double_threec(second_col)
168 np.testing.assert_array_equal(counting_mat, [[0.0, 2, 2], [6, 16, 10], [6, 14, 8]])
172 """Eigen doesn't support (as of yet) negative strides. When a function takes an Eigen matrix by
173 copy or const reference, we can pass a numpy array that has negative strides. Otherwise, an
174 exception will be thrown as Eigen will not be able to map the numpy array."""
176 counting_mat = np.arange(9.0, dtype=np.float32).
reshape((3, 3))
177 counting_mat = counting_mat[::-1, ::-1]
178 second_row = counting_mat[1, :]
179 second_col = counting_mat[:, 1]
180 np.testing.assert_array_equal(m.double_row(second_row), 2.0 * second_row)
181 np.testing.assert_array_equal(m.double_col(second_row), 2.0 * second_row)
182 np.testing.assert_array_equal(m.double_complex(second_row), 2.0 * second_row)
183 np.testing.assert_array_equal(m.double_row(second_col), 2.0 * second_col)
184 np.testing.assert_array_equal(m.double_col(second_col), 2.0 * second_col)
185 np.testing.assert_array_equal(m.double_complex(second_col), 2.0 * second_col)
187 counting_3d = np.arange(27.0, dtype=np.float32).
reshape((3, 3, 3))
188 counting_3d = counting_3d[::-1, ::-1, ::-1]
189 slices = [counting_3d[0, :, :], counting_3d[:, 0, :], counting_3d[:, :, 0]]
190 for ref_mat
in slices:
191 np.testing.assert_array_equal(m.double_mat_cm(ref_mat), 2.0 * ref_mat)
192 np.testing.assert_array_equal(m.double_mat_rm(ref_mat), 2.0 * ref_mat)
195 with pytest.raises(TypeError)
as excinfo:
196 m.double_threer(second_row)
200 double_threer(): incompatible function arguments. The following argument types are supported:
201 1. (arg0: numpy.ndarray[numpy.float32[1, 3], flags.writeable]) -> None
204 +
repr(np.array([5.0, 4.0, 3.0], dtype=
"float32"))
207 with pytest.raises(TypeError)
as excinfo:
208 m.double_threec(second_col)
212 double_threec(): incompatible function arguments. The following argument types are supported:
213 1. (arg0: numpy.ndarray[numpy.float32[3, 1], flags.writeable]) -> None
216 +
repr(np.array([7.0, 4.0, 1.0], dtype=
"float32"))
221 assert np.all(m.diagonal(ref) == ref.diagonal())
222 assert np.all(m.diagonal_1(ref) == ref.diagonal(1))
223 for i
in range(-5, 7):
224 assert np.all(m.diagonal_n(ref, i) == ref.diagonal(i)), f
"m.diagonal_n({i})"
226 assert np.all(m.block(ref, 2, 1, 3, 3) == ref[2:5, 1:4])
227 assert np.all(m.block(ref, 1, 4, 4, 2) == ref[1:, 4:])
228 assert np.all(m.block(ref, 1, 4, 3, 2) == ref[1:4, 4:])
232 chols = [m.cholesky1, m.cholesky2, m.cholesky3, m.cholesky4]
233 for i, chol
in enumerate(chols, start=1):
234 mymat = chol(np.array([[1.0, 2, 4], [2, 13, 23], [4, 23, 77]]))
236 mymat == np.array([[1, 0, 0], [2, 3, 0], [4, 5, 6]])
246 z = np.array(a, copy=
True)
252 """Tests various ways of returning references and non-referencing copies"""
254 master = np.ones((10, 10))
257 assert not a_get1.flags.owndata
and a_get1.flags.writeable
260 assert not a_get2.flags.owndata
and a_get2.flags.writeable
264 assert not a_view1.flags.owndata
and not a_view1.flags.writeable
265 with pytest.raises(ValueError):
267 a_view2 = a.view_ptr()
268 assert not a_view2.flags.owndata
and not a_view2.flags.writeable
269 with pytest.raises(ValueError):
272 a_copy1 = a.copy_get()
273 assert a_copy1.flags.owndata
and a_copy1.flags.writeable
274 np.testing.assert_array_equal(a_copy1, master)
277 a_copy2 = a.copy_view()
278 assert a_copy2.flags.owndata
and a_copy2.flags.writeable
279 np.testing.assert_array_equal(a_copy2, master)
284 assert not a_ref1.flags.owndata
and a_ref1.flags.writeable
286 a_ref2 = a.ref_const()
287 assert not a_ref2.flags.owndata
and not a_ref2.flags.writeable
288 with pytest.raises(ValueError):
290 a_ref3 = a.ref_safe()
291 assert not a_ref3.flags.owndata
and a_ref3.flags.writeable
293 a_ref4 = a.ref_const_safe()
294 assert not a_ref4.flags.owndata
and not a_ref4.flags.writeable
295 with pytest.raises(ValueError):
296 a_ref4[7, 0] = 987654321
298 a_copy3 = a.copy_ref()
299 assert a_copy3.flags.owndata
and a_copy3.flags.writeable
300 np.testing.assert_array_equal(a_copy3, master)
303 a_copy4 = a.copy_ref_const()
304 assert a_copy4.flags.owndata
and a_copy4.flags.writeable
305 np.testing.assert_array_equal(a_copy4, master)
309 a_block1 = a.block(3, 3, 2, 2)
310 assert not a_block1.flags.owndata
and a_block1.flags.writeable
313 a_block2 = a.block_safe(2, 2, 3, 2)
314 assert not a_block2.flags.owndata
and a_block2.flags.writeable
315 a_block2[2, 1] = -123
317 a_block3 = a.block_const(6, 7, 4, 3)
318 assert not a_block3.flags.owndata
and not a_block3.flags.writeable
319 with pytest.raises(ValueError):
320 a_block3[2, 2] = -44444
322 a_copy5 = a.copy_block(2, 2, 2, 3)
323 assert a_copy5.flags.owndata
and a_copy5.flags.writeable
324 np.testing.assert_array_equal(a_copy5, master[2:4, 2:5])
328 a_corn1 = a.corners()
329 assert not a_corn1.flags.owndata
and a_corn1.flags.writeable
336 a_corn2 = a.corners_const()
337 assert not a_corn2.flags.owndata
and not a_corn2.flags.writeable
338 with pytest.raises(ValueError):
343 np.testing.assert_array_equal(a_get1, master)
344 np.testing.assert_array_equal(a_get2, master)
345 np.testing.assert_array_equal(a_view1, master)
346 np.testing.assert_array_equal(a_view2, master)
347 np.testing.assert_array_equal(a_ref1, master)
348 np.testing.assert_array_equal(a_ref2, master)
349 np.testing.assert_array_equal(a_ref3, master)
350 np.testing.assert_array_equal(a_ref4, master)
351 np.testing.assert_array_equal(a_block1, master[3:5, 3:5])
352 np.testing.assert_array_equal(a_block2, master[2:5, 2:4])
353 np.testing.assert_array_equal(a_block3, master[6:10, 7:10])
354 np.testing.assert_array_equal(
355 a_corn1, master[0 :: master.shape[0] - 1, 0 :: master.shape[1] - 1]
357 np.testing.assert_array_equal(
358 a_corn2, master[0 :: master.shape[0] - 1, 0 :: master.shape[1] - 1]
361 np.testing.assert_array_equal(a_copy1, c1want)
362 np.testing.assert_array_equal(a_copy2, c2want)
363 np.testing.assert_array_equal(a_copy3, c3want)
364 np.testing.assert_array_equal(a_copy4, c4want)
365 np.testing.assert_array_equal(a_copy5, c5want)
370 start_with = cstats.alive()
372 assert cstats.alive() == start_with + 1
374 assert cstats.alive() == start_with + 1
377 assert cstats.alive() == start_with + 1
380 assert cstats.alive() == start_with
386 assert cstats.alive() == 1
387 unsafe = [a.ref(), a.ref_const(), a.block(1, 2, 3, 4)]
393 a.copy_block(4, 3, 2, 1),
396 assert cstats.alive() == 0
402 m.ReturnTester.get_ptr,
404 m.ReturnTester.view_ptr,
405 m.ReturnTester.ref_safe,
406 m.ReturnTester.ref_const_safe,
407 m.ReturnTester.corners,
408 m.ReturnTester.corners_const,
412 for meth
in [m.ReturnTester.block_safe, m.ReturnTester.block_const]:
417 """Tests Eigen's ability to mutate numpy values"""
419 orig = np.array([[1.0, 2, 3], [4, 5, 6], [7, 8, 9]])
421 zc = np.array(orig, order=
"F")
422 m.add_rm(zr, 1, 0, 100)
423 assert np.all(zr == np.array([[1.0, 2, 3], [104, 5, 6], [7, 8, 9]]))
424 m.add_cm(zc, 1, 0, 200)
425 assert np.all(zc == np.array([[1.0, 2, 3], [204, 5, 6], [7, 8, 9]]))
427 m.add_any(zr, 1, 0, 20)
428 assert np.all(zr == np.array([[1.0, 2, 3], [124, 5, 6], [7, 8, 9]]))
429 m.add_any(zc, 1, 0, 10)
430 assert np.all(zc == np.array([[1.0, 2, 3], [214, 5, 6], [7, 8, 9]]))
433 with pytest.raises(TypeError):
434 m.add_rm(zc, 1, 0, 1)
435 with pytest.raises(TypeError):
436 m.add_cm(zr, 1, 0, 1)
439 m.add1(zr, 1, 0, -100)
440 m.add2(zr, 1, 0, -20)
441 assert np.all(zr == orig)
442 m.add1(zc, 1, 0, -200)
443 m.add2(zc, 1, 0, -10)
444 assert np.all(zc == orig)
448 cornersr = zr[0::2, 0::2]
449 cornersc = zc[0::2, 0::2]
451 assert np.all(cornersr == np.array([[1.0, 3], [7, 9]]))
452 assert np.all(cornersc == np.array([[1.0, 3], [7, 9]]))
454 with pytest.raises(TypeError):
455 m.add_rm(cornersr, 0, 1, 25)
456 with pytest.raises(TypeError):
457 m.add_cm(cornersr, 0, 1, 25)
458 with pytest.raises(TypeError):
459 m.add_rm(cornersc, 0, 1, 25)
460 with pytest.raises(TypeError):
461 m.add_cm(cornersc, 0, 1, 25)
462 m.add_any(cornersr, 0, 1, 25)
463 m.add_any(cornersc, 0, 1, 44)
464 assert np.all(zr == np.array([[1.0, 2, 28], [4, 5, 6], [7, 8, 9]]))
465 assert np.all(zc == np.array([[1.0, 2, 47], [4, 5, 6], [7, 8, 9]]))
469 zro.flags.writeable =
False
470 with pytest.raises(TypeError):
471 m.add_rm(zro, 0, 0, 0)
472 with pytest.raises(TypeError):
473 m.add_any(zro, 0, 0, 0)
474 with pytest.raises(TypeError):
476 with pytest.raises(TypeError):
480 zi = np.array([[1, 2], [3, 4]])
481 with pytest.raises(TypeError):
486 """Tests numpy mutating Eigen matrices (for returned Eigen::Ref<...>s)"""
491 zcro = m.get_cm_const_ref()
493 zrro = m.get_rm_const_ref()
495 assert [zc[1, 2], zcro[1, 2], zr[1, 2], zrro[1, 2]] == [23] * 4
497 assert not zc.flags.owndata
and zc.flags.writeable
498 assert not zr.flags.owndata
and zr.flags.writeable
499 assert not zcro.flags.owndata
and not zcro.flags.writeable
500 assert not zrro.flags.owndata
and not zrro.flags.writeable
503 expect = np.array([[11.0, 12, 13], [21, 22, 99], [31, 32, 33]])
505 assert np.all(zc == expect)
506 assert np.all(zcro == expect)
507 assert np.all(m.get_cm_ref() == expect)
510 assert np.all(zr == expect)
511 assert np.all(zrro == expect)
512 assert np.all(m.get_rm_ref() == expect)
515 with pytest.raises(ValueError):
517 with pytest.raises(ValueError):
522 y1 = np.array(m.get_cm_const_ref())
524 assert y1.flags.owndata
and y1.flags.writeable
526 assert y1[1, 2] == 99
528 assert y1[1, 2] == 111
529 assert zc[1, 2] == 99
533 """Tests a complex chain of nested eigen/numpy references"""
539 z2 = m.incr_matrix(z, 1)
541 z3 = m.incr_matrix(z, 2)
543 z4 = m.incr_matrix(z, 3)
545 z5 = m.incr_matrix(z, 4)
547 assert np.all(z == z2)
548 assert np.all(z == z3)
549 assert np.all(z == z4)
550 assert np.all(z == z5)
551 expect = np.array([[0.0, 22, 20], [31, 37, 33], [41, 42, 38]])
552 assert np.all(z == expect)
554 y = np.array(
range(100), dtype=
"float64").
reshape(10, 10)
555 y2 = m.incr_matrix_any(y, 10)
556 y3 = m.incr_matrix_any(
561 y6 = m.incr_matrix_any(y5, 1000)
564 yexpect = np.array(
range(100), dtype=
"float64").
reshape(10, 10)
566 yexpect[0::2, 0::2] -= 33
567 yexpect[0::4, 0::4] += 1000
568 assert np.all(y6 == yexpect[0::4, 0::4])
569 assert np.all(y5 == yexpect[0::4, 0::4])
570 assert np.all(y4 == yexpect[0::4, 0::2])
571 assert np.all(y3 == yexpect[0::2, 0::2])
572 assert np.all(y2 == yexpect)
573 assert np.all(y == yexpect)
579 int_matrix_colmajor = np.array([[1, 2, 3], [4, 5, 6], [7, 8, 9]], order=
"F")
580 dbl_matrix_colmajor = np.array(
581 int_matrix_colmajor, dtype=
"double", order=
"F", copy=
True
583 int_matrix_rowmajor = np.array(int_matrix_colmajor, order=
"C", copy=
True)
584 dbl_matrix_rowmajor = np.array(
585 int_matrix_rowmajor, dtype=
"double", order=
"C", copy=
True
589 assert m.get_elem(int_matrix_colmajor) == 8
590 assert m.get_elem(dbl_matrix_colmajor) == 8
591 assert m.get_elem(int_matrix_rowmajor) == 8
592 assert m.get_elem(dbl_matrix_rowmajor) == 8
595 with pytest.raises(TypeError)
as excinfo:
596 m.get_elem_nocopy(int_matrix_colmajor)
597 assert "get_elem_nocopy(): incompatible function arguments." in str(
599 )
and ", flags.f_contiguous" in str(excinfo.value)
600 assert m.get_elem_nocopy(dbl_matrix_colmajor) == 8
601 with pytest.raises(TypeError)
as excinfo:
602 m.get_elem_nocopy(int_matrix_rowmajor)
603 assert "get_elem_nocopy(): incompatible function arguments." in str(
605 )
and ", flags.f_contiguous" in str(excinfo.value)
606 with pytest.raises(TypeError)
as excinfo:
607 m.get_elem_nocopy(dbl_matrix_rowmajor)
608 assert "get_elem_nocopy(): incompatible function arguments." in str(
610 )
and ", flags.f_contiguous" in str(excinfo.value)
613 with pytest.raises(TypeError)
as excinfo:
614 m.get_elem_rm_nocopy(int_matrix_colmajor)
615 assert "get_elem_rm_nocopy(): incompatible function arguments." in str(
617 )
and ", flags.c_contiguous" in str(excinfo.value)
618 with pytest.raises(TypeError)
as excinfo:
619 m.get_elem_rm_nocopy(dbl_matrix_colmajor)
620 assert "get_elem_rm_nocopy(): incompatible function arguments." in str(
622 )
and ", flags.c_contiguous" in str(excinfo.value)
623 assert m.get_elem_rm_nocopy(int_matrix_rowmajor) == 8
624 with pytest.raises(TypeError)
as excinfo:
625 m.get_elem_rm_nocopy(dbl_matrix_rowmajor)
626 assert "get_elem_rm_nocopy(): incompatible function arguments." in str(
628 )
and ", flags.c_contiguous" in str(excinfo.value)
632 """Ensure the lifetime of temporary arrays created by the `Ref` caster
634 The `Ref` caster sometimes creates a copy which needs to stay alive. This needs to
635 happen both for directs casts (just the array) or indirectly (e.g. list of arrays).
638 a = np.full(shape=10, fill_value=8, dtype=np.int8)
639 assert m.get_elem_direct(a) == 8
642 assert m.get_elem_indirect(list_of_a) == 8
646 assert np.all(m.incr_diag(7) == np.diag([1.0, 2, 3, 4, 5, 6, 7]))
648 asymm = np.array([[1.0, 2, 3, 4], [5, 6, 7, 8], [9, 10, 11, 12], [13, 14, 15, 16]])
649 symm_lower = np.array(asymm)
650 symm_upper = np.array(asymm)
652 for j
in range(i + 1, 4):
653 symm_lower[i, j] = symm_lower[j, i]
654 symm_upper[j, i] = symm_upper[i, j]
656 assert np.all(m.symmetric_lower(asymm) == symm_lower)
657 assert np.all(m.symmetric_upper(asymm) == symm_upper)
664 double_col(arg0: numpy.ndarray[numpy.float32[m, 1]]) -> numpy.ndarray[numpy.float32[m, 1]]
670 double_row(arg0: numpy.ndarray[numpy.float32[1, n]]) -> numpy.ndarray[numpy.float32[1, n]]
673 assert doc(m.double_complex) == (
675 double_complex(arg0: numpy.ndarray[numpy.complex64[m, 1]])"""
676 """ -> numpy.ndarray[numpy.complex64[m, 1]]
679 assert doc(m.double_mat_rm) == (
681 double_mat_rm(arg0: numpy.ndarray[numpy.float32[m, n]])"""
682 """ -> numpy.ndarray[numpy.float32[m, n]]
688 a = np.array([[1.0, 2], [3, 4], [5, 6]])
691 assert np.all(m.matrix_multiply(a, b) == np.array([[3.0], [7], [11]]))
692 assert np.all(m.matrix_multiply(A=a, B=b) == np.array([[3.0], [7], [11]]))
693 assert np.all(m.matrix_multiply(B=b, A=a) == np.array([[3.0], [7], [11]]))
695 with pytest.raises(ValueError)
as excinfo:
696 m.matrix_multiply(b, a)
697 assert str(excinfo.value) ==
"Nonconformable matrices!"
699 with pytest.raises(ValueError)
as excinfo:
700 m.matrix_multiply(A=b, B=a)
701 assert str(excinfo.value) ==
"Nonconformable matrices!"
703 with pytest.raises(ValueError)
as excinfo:
704 m.matrix_multiply(B=a, A=b)
705 assert str(excinfo.value) ==
"Nonconformable matrices!"
709 pytest.importorskip(
"scipy")
719 pytest.importorskip(
"scipy")
723 sparse_copy_r(arg0: scipy.sparse.csr_matrix[numpy.float32]) -> scipy.sparse.csr_matrix[numpy.float32]
729 sparse_copy_c(arg0: scipy.sparse.csc_matrix[numpy.float32]) -> scipy.sparse.csc_matrix[numpy.float32]
735 """Ignore strides on a length-1 dimension (even if they would be incompatible length > 1)"""
736 assert np.all(m.iss738_f1(np.array([[1.0, 2, 3]])) == np.array([[1.0, 102, 203]]))
738 m.iss738_f1(np.array([[1.0], [2], [3]])) == np.array([[1.0], [12], [23]])
741 assert np.all(m.iss738_f2(np.array([[1.0, 2, 3]])) == np.array([[1.0, 102, 203]]))
743 m.iss738_f2(np.array([[1.0], [2], [3]])) == np.array([[1.0], [12], [23]])
747 @pytest.mark.parametrize(
"func", [m.iss738_f1, m.iss738_f2])
748 @pytest.mark.parametrize(
"sizes", [(0, 2), (2, 0)])
750 """Ignore strides on a length-0 dimension (even if they would be incompatible length > 1)"""
751 assert np.all(
func(np.zeros(sizes)) == np.zeros(sizes))
755 """Issue 1105: 1xN or Nx1 input arrays weren't accepted for eigen
756 compile-time row vectors or column vector"""
757 assert m.iss1105_row(np.ones((1, 7)))
758 assert m.iss1105_col(np.ones((7, 1)))
761 with pytest.raises(TypeError)
as excinfo:
762 m.iss1105_row(np.ones((7, 1)))
763 assert "incompatible function arguments" in str(excinfo.value)
764 with pytest.raises(TypeError)
as excinfo:
765 m.iss1105_col(np.ones((1, 7)))
766 assert "incompatible function arguments" in str(excinfo.value)
770 """Using Eigen types as member variables requires a class-specific
771 operator new with proper alignment"""
773 o = m.CustomOperatorNew()
774 np.testing.assert_allclose(o.a, 0.0)
775 np.testing.assert_allclose(o.b.diagonal(), 1.0)