Re: [eigen] NEON remaining failing tests

[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]


For cholesky_2:

Confirmed on x86-64 linux...

Backtrace:

Program received signal SIGABRT, Aborted.
0x00007ffff6cf4035 in raise () from /lib/libc.so.6
(gdb) bt
#0  0x00007ffff6cf4035 in raise () from /lib/libc.so.6
#1  0x00007ffff6cf5460 in abort () from /lib/libc.so.6
#2  0x00007ffff6d2e3f7 in __libc_message () from /lib/libc.so.6
#3  0x00007ffff6d33876 in malloc_printerr () from /lib/libc.so.6
#4  0x00007ffff6d385ac in free () from /lib/libc.so.6
#5  0x0000000000483efb in Eigen::ei_aligned_free (ptr=0x7db390) at
/home/bjacob/eigen/Eigen/src/Core/util/Memory.h:202
#6  0x0000000000493f31 in Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run (size=184, depth=16, _mat=0x78d200,
    matStride=200, res=0x793600, resStride=200, alpha=-1)
    at /home/bjacob/eigen/Eigen/src/Core/products/SelfadjointProduct.h:116
#7  0x00000000004900e8 in
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32> > (
    this=0x7fffffffddd0, u=..., alpha=-1) at
/home/bjacob/eigen/Eigen/src/Core/products/SelfadjointProduct.h:137
#8  0x000000000048ce5c in
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> > (m=
    Eigen::Matrix<double,33331,33331> = {...}) at
/home/bjacob/eigen/Eigen/src/Cholesky/LLT.h:204
#9  0x000000000048a39d in Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>, 1>::inplace_decomposition
(m=Eigen::Matrix<double,33331,33331> = {...}) at
/home/bjacob/eigen/Eigen/src/Cholesky/LLT.h:233
#10 0x0000000000487eaa in Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute (
    this=0x7fffffffdf90, a=Eigen::Matrix<double,33331,33331> = {...})
at /home/bjacob/eigen/Eigen/src/Cholesky/LLT.h:259
#11 0x0000000000486341 in LLT (this=0x7fffffffdf90,
matrix=Eigen::Matrix<double,33331,33331> = {...})
    at /home/bjacob/eigen/Eigen/src/Cholesky/LLT.h:84
#12 0x000000000048475e in cholesky<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> > (m=
    Eigen::Matrix<double,33331,33331> = {...}) at
/home/bjacob/eigen/test/cholesky.cpp:97
#13 0x00000000004804ab in test_cholesky () at
/home/bjacob/eigen/test/cholesky.cpp:158
#14 0x00000000004802e1 in main (argc=1, argv=0x7fffffffe468) at
/home/bjacob/eigen/test/main.h:539





and now the valgrind output. Notice that the first error is an invalid
write in ei_gebp_kernel().


==5300== Memcheck, a memory error detector
==5300== Copyright (C) 2002-2009, and GNU GPL'd, by Julian Seward et al.
==5300== Using Valgrind-3.5.0 and LibVEX; rerun with -h for copyright info
==5300== Command: ./test/cholesky_2
==5300==
Initializing random number generator with seed 1267903652
Repeating each test 10 times
==5300== Invalid write of size 8
==5300==    at 0x4A0B5D: void Eigen::ei_pstore<double,
double>(double*, double const&) (GenericPacketMath.h:166)
==5300==    by 0x4997BC: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:65)
==5300==    by 0x49C352: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:184)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217c0 is 0 bytes after a block of size 23,552 alloc'd
==5300==    at 0x4C23E03: malloc (in
/usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==5300==    by 0x483EA2: Eigen::ei_aligned_malloc(unsigned long) (Memory.h:178)
==5300==    by 0x493C1D: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:76)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x499E66: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:223)
==5300==    by 0x49C352: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:184)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217c0 is 0 bytes after a block of size 23,552 alloc'd
==5300==    at 0x4C23E03: malloc (in
/usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==5300==    by 0x483EA2: Eigen::ei_aligned_malloc(unsigned long) (Memory.h:178)
==5300==    by 0x493C1D: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:76)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x499F2F: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:226)
==5300==    by 0x49C352: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:184)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217c8 is 8 bytes after a block of size 23,552 alloc'd
==5300==    at 0x4C23E03: malloc (in
/usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==5300==    by 0x483EA2: Eigen::ei_aligned_malloc(unsigned long) (Memory.h:178)
==5300==    by 0x493C1D: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:76)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x499FF8: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:229)
==5300==    by 0x49C352: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:184)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217d0 is not stack'd, malloc'd or (recently) free'd
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x49A0F7: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:234)
==5300==    by 0x49C352: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:184)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217d8 is not stack'd, malloc'd or (recently) free'd
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x49A1C0: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:238)
==5300==    by 0x49C352: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:184)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217e0 is not stack'd, malloc'd or (recently) free'd
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x49A289: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:241)
==5300==    by 0x49C352: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:184)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217e8 is not stack'd, malloc'd or (recently) free'd
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x49A352: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:244)
==5300==    by 0x49C352: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:184)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217f0 is 16 bytes before a block of size 23,552 alloc'd
==5300==    at 0x4C23E03: malloc (in
/usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==5300==    by 0x483EA2: Eigen::ei_aligned_malloc(unsigned long) (Memory.h:178)
==5300==    by 0x493CB1: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:77)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x49A451: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:249)
==5300==    by 0x49C352: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:184)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217f8 is 8 bytes before a block of size 23,552 alloc'd
==5300==    at 0x4C23E03: malloc (in
/usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==5300==    by 0x483EA2: Eigen::ei_aligned_malloc(unsigned long) (Memory.h:178)
==5300==    by 0x493CB1: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:77)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==
==5300== Invalid write of size 8
==5300==    at 0x4A0B5D: void Eigen::ei_pstore<double,
double>(double*, double const&) (GenericPacketMath.h:166)
==5300==    by 0x4997BC: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:65)
==5300==    by 0x49C4B2: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:198)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217c0 is 0 bytes after a block of size 23,552 alloc'd
==5300==    at 0x4C23E03: malloc (in
/usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==5300==    by 0x483EA2: Eigen::ei_aligned_malloc(unsigned long) (Memory.h:178)
==5300==    by 0x493C1D: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:76)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x499E66: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:223)
==5300==    by 0x49C4B2: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:198)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217c0 is 0 bytes after a block of size 23,552 alloc'd
==5300==    at 0x4C23E03: malloc (in
/usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==5300==    by 0x483EA2: Eigen::ei_aligned_malloc(unsigned long) (Memory.h:178)
==5300==    by 0x493C1D: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:76)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x499F2F: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:226)
==5300==    by 0x49C4B2: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:198)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217c8 is 8 bytes after a block of size 23,552 alloc'd
==5300==    at 0x4C23E03: malloc (in
/usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==5300==    by 0x483EA2: Eigen::ei_aligned_malloc(unsigned long) (Memory.h:178)
==5300==    by 0x493C1D: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:76)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x499FF8: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:229)
==5300==    by 0x49C4B2: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:198)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217d0 is not stack'd, malloc'd or (recently) free'd
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x49A0F7: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:234)
==5300==    by 0x49C4B2: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:198)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217d8 is not stack'd, malloc'd or (recently) free'd
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x49A1C0: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:238)
==5300==    by 0x49C4B2: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:198)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217e0 is not stack'd, malloc'd or (recently) free'd
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x49A289: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:241)
==5300==    by 0x49C4B2: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:198)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217e8 is not stack'd, malloc'd or (recently) free'd
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x49A352: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:244)
==5300==    by 0x49C4B2: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:198)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217f0 is 16 bytes before a block of size 23,552 alloc'd
==5300==    at 0x4C23E03: malloc (in
/usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==5300==    by 0x483EA2: Eigen::ei_aligned_malloc(unsigned long) (Memory.h:178)
==5300==    by 0x493CB1: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:77)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==
==5300== Invalid read of size 8
==5300==    at 0x4A0B6E: Eigen::ei_packet_traits<double>::type
Eigen::ei_pload<double>(double const*) (GenericPacketMath.h:150)
==5300==    by 0x49A451: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:249)
==5300==    by 0x49C4B2: Eigen::ei_sybb_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false>, 1>::operator()(double*, int,
double const*, double const*, int, int) (SelfadjointProduct.h:198)
==5300==    by 0x493EB7: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:106)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==  Address 0x63217f8 is 8 bytes before a block of size 23,552 alloc'd
==5300==    at 0x4C23E03: malloc (in
/usr/lib/valgrind/vgpreload_memcheck-amd64-linux.so)
==5300==    by 0x483EA2: Eigen::ei_aligned_malloc(unsigned long) (Memory.h:178)
==5300==    by 0x493CB1: Eigen::ei_selfadjoint_product<double, 0, 0,
true, 1>::run(int, int, double const*, int, double*, int, double)
(SelfadjointProduct.h:77)
==5300==    by 0x4900E7:
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>, 1u>&
Eigen::SelfAdjointView<Eigen::Block<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 33331, 33331, 32>,
1u>::rankUpdate<Eigen::Block<Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>, 33331, 33331, 32>
>(Eigen::MatrixBase<Eigen::Block<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331>, 33331, 33331, 32> > const&, double)
(SelfadjointProduct.h:137)
==5300==    by 0x48CE5B: bool
Eigen::ei_llt_inplace<1>::blocked<Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0, 33331,
33331>&) (LLT.h:204)
==5300==    by 0x48A39C: Eigen::LLT_Traits<Eigen::Matrix<double,
33331, 33331, 0, 33331, 33331>,
1>::inplace_decomposition(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331>&) (LLT.h:233)
==5300==    by 0x487EA9: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::compute(Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> const&) (LLT.h:259)
==5300==    by 0x486340: Eigen::LLT<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331>, 1>::LLT(Eigen::Matrix<double, 33331, 33331,
0, 33331, 33331> const&) (LLT.h:84)
==5300==    by 0x48475D: void cholesky<Eigen::Matrix<double, 33331,
33331, 0, 33331, 33331> >(Eigen::Matrix<double, 33331, 33331, 0,
33331, 33331> const&) (cholesky.cpp:97)
==5300==    by 0x4804AA: test_cholesky() (cholesky.cpp:158)
==5300==    by 0x4802E0: main (main.h:539)
==5300==
==5300== Invalid write of size 8
==5300==    at 0x4A0B5D: void Eigen::ei_pstore<double,
double>(double*, double const&) (GenericPacketMath.h:166)
==5300==    by 0x4997BC: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:65)
==5300==    by 0x3FE1B73113CF679F: ???
==5300==    by 0xBFD30321E277961A: ???
==5300==    by 0xBFBC48C821DF473A: ???
==5300==    by 0x3F78D88C95F99AB6: ???
==5300==    by 0x3FC913A601A38E92: ???
==5300==    by 0xBFB0803A98AD7C7D: ???
==5300==    by 0xBFED5438EC96B55E: ???
==5300==    by 0x3FC461E37238A370: ???
==5300==  Address 0x3fc461e37238a409 is not stack'd, malloc'd or
(recently) free'd
==5300==
==5300==
==5300== Process terminating with default action of signal 11 (SIGSEGV)
==5300==  General Protection Fault
==5300==    at 0x4A0B5D: void Eigen::ei_pstore<double,
double>(double*, double const&) (GenericPacketMath.h:166)
==5300==    by 0x4997BC: Eigen::ei_gebp_kernel<double, 2, 4,
Eigen::ei_conj_helper<false, false> >::operator()(double*, int, double
const*, double const*, int, int, int, int, int, int, int, double*)
(GeneralBlockPanelKernel.h:65)
==5300==    by 0x3FE1B73113CF679F: ???
==5300==    by 0xBFD30321E277961A: ???
==5300==    by 0xBFBC48C821DF473A: ???
==5300==    by 0x3F78D88C95F99AB6: ???
==5300==    by 0x3FC913A601A38E92: ???
==5300==    by 0xBFB0803A98AD7C7D: ???
==5300==    by 0xBFED5438EC96B55E: ???
==5300==    by 0x3FC461E37238A370: ???
==5300==
==5300== HEAP SUMMARY:
==5300==     in use at exit: 2,563,301 bytes in 13 blocks
==5300==   total heap usage: 70 allocs, 57 frees, 7,169,077 bytes allocated
==5300==
==5300== LEAK SUMMARY:
==5300==    definitely lost: 0 bytes in 0 blocks
==5300==    indirectly lost: 0 bytes in 0 blocks
==5300==      possibly lost: 85 bytes in 2 blocks
==5300==    still reachable: 2,563,216 bytes in 11 blocks
==5300==         suppressed: 0 bytes in 0 blocks
==5300== Rerun with --leak-check=full to see details of leaked memory
==5300==
==5300== For counts of detected and suppressed errors, rerun with: -v
==5300== ERROR SUMMARY: 1441 errors from 19 contexts (suppressed: 7 from 7)
Segmentation fault



Benoit


2010/3/6 Konstantinos Margaritis <markos@xxxxxxxx>:
> On Saturday 06 March 2010 13:53:00 Konstantinos Margaritis wrote:
>> make in test/ fails in cholesky:
>>
>> [ 12%] Building CXX object test/CMakeFiles/cholesky_1.dir/cholesky.cpp.o
>> In file included from /home/oem/eigen/Eigen/Core:272,
>>                  from /home/oem/eigen/Eigen/QR:4,
>>                  from /home/oem/eigen/test/main.h:151,
>>                  from /home/oem/eigen/test/cholesky.cpp:29:
>> /home/oem/eigen/Eigen/src/Core/Product.h: In instantiation of
>> 'Eigen::ei_product_type<Eigen::Block<Eigen::Matrix<double, 1, 1, 0, 1, 1>,
>> 33331, 33331, 32>, Eigen::Block<Eigen::Block<Eigen::Matrix<double, 1, 1, 0,
>>  1, 1>, 1, 1, 32>, 1, 1, 32> >':
>> /home/oem/eigen/Eigen/src/Cholesky/LDLT.h:249:   instantiated from
>> 'Eigen::LDLT<MatrixType>& Eigen::LDLT<MatrixType>::compute(const
>>  _MatrixType&) [with _MatrixType = Eigen::Matrix<double, 1, 1, 0, 1, 1>]'
>> /home/oem/eigen/Eigen/src/Cholesky/LDLT.h:86:   instantiated from
>> 'Eigen::LDLT<MatrixType>::LDLT(const _MatrixType&) [with _MatrixType =
>> Eigen::Matrix<double, 1, 1, 0, 1, 1>]'
>> /home/oem/eigen/test/cholesky.cpp:121:   instantiated from 'void
>> cholesky(const MatrixType&) [with MatrixType = Eigen::Matrix<double, 1, 1,
>>  0, 1, 1>]'
>> /home/oem/eigen/test/cholesky.cpp:153:   instantiated from here
>> /home/oem/eigen/Eigen/src/Core/Product.h:78: error: incomplete type
>> 'Eigen::ei_product_type_selector<2, 1, 1>' used in nested name specifier
>
> Uhm I just reproduced this sans-NEON, and also on x86 and sse2 builds.
> Also, cholesky_2 builds fine but I get the following on ARM(with/without NEON):
>
> $ ./cholesky_2
> Initializing random number generator with seed 1267893603
> Repeating each test 10 times
> *** glibc detected *** ./cholesky_2: double free or corruption (!prev): 0x00132750 ***
> Aborted
>
> and on x86(both with/without sse):
>
> $ ./cholesky_2
> Initializing random number generator with seed 1267894057
> Repeating each test 10 times
> *** glibc detected *** ./cholesky_2: double free or corruption (!prev): 0x000000000072d390 ***
> ======= Backtrace: =========
> /lib64/libc.so.6[0x7f09f4e0cc76]
> /lib64/libc.so.6(cfree+0x6c)[0x7f09f4e1196c]
> ./cholesky_2(_ZN5Eigen22ei_selfadjoint_productIdLi0ELi0ELb1ELi1EE3runEiiPKdiPdid+0x6c7)[0x422207]
> ./cholesky_2(_ZN5Eigen14ei_llt_inplaceILi1EE7blockedINS_6MatrixIdLi33331ELi33331ELi0ELi33331ELi33331EEEEEbRT_+0x160)
> [0x429b20]
> ./cholesky_2(_ZN5Eigen3LLTINS_6MatrixIdLi33331ELi33331ELi0ELi33331ELi33331EEELi1EE7computeERKS2_+0x13b)
> [0x429d5b]
> ./cholesky_2(_Z8choleskyIN5Eigen6MatrixIdLi33331ELi33331ELi0ELi33331ELi33331EEEEvRKT_+0xc06)[0x42a9f6]
> ./cholesky_2(_Z13test_choleskyv+0x17d)[0x4180cd]
> ./cholesky_2(main+0x1f0)[0x4185d0]
> /lib64/libc.so.6(__libc_start_main+0xfd)[0x7f09f4db8a7d]
> ./cholesky_2[0x40da69]
> ======= Memory map: ========
> 00400000-00441000 r-xp 00000000 08:06 2113301                            /home/markos/Development/eigen-build-
> plain/test/cholesky_2
> 00640000-00641000 r--p 00040000 08:06 2113301                            /home/markos/Development/eigen-build-
> plain/test/cholesky_2
> 00641000-00642000 rw-p 00041000 08:06 2113301                            /home/markos/Development/eigen-build-
> plain/test/cholesky_2
> 00642000-0074e000 rw-p 00000000 00:00 0                                  [heap]
> 7f09f0000000-7f09f0021000 rw-p 00000000 00:00 0
> 7f09f0021000-7f09f4000000 ---p 00000000 00:00 0
> 7f09f4d4b000-7f09f4d9a000 rw-p 00000000 00:00 0
> 7f09f4d9a000-7f09f4eeb000 r-xp 00000000 08:05 2490380                    /lib64/libc-2.10.1.so
> 7f09f4eeb000-7f09f50eb000 ---p 00151000 08:05 2490380                    /lib64/libc-2.10.1.so
> 7f09f50eb000-7f09f50ef000 r--p 00151000 08:05 2490380                    /lib64/libc-2.10.1.so
> 7f09f50ef000-7f09f50f0000 rw-p 00155000 08:05 2490380                    /lib64/libc-2.10.1.so
> 7f09f50f0000-7f09f50f5000 rw-p 00000000 00:00 0
> 7f09f50f5000-7f09f510b000 r-xp 00000000 08:05 2490422                    /lib64/libgcc_s.so.1
> 7f09f510b000-7f09f530a000 ---p 00016000 08:05 2490422                    /lib64/libgcc_s.so.1
> 7f09f530a000-7f09f530b000 r--p 00015000 08:05 2490422                    /lib64/libgcc_s.so.1
> 7f09f530b000-7f09f530c000 rw-p 00016000 08:05 2490422                    /lib64/libgcc_s.so.1
> 7f09f530c000-7f09f5360000 r-xp 00000000 08:05 2490412                    /lib64/libm-2.10.1.so
> 7f09f5360000-7f09f555f000 ---p 00054000 08:05 2490412                    /lib64/libm-2.10.1.so
> 7f09f555f000-7f09f5560000 r--p 00053000 08:05 2490412                    /lib64/libm-2.10.1.so
> 7f09f5560000-7f09f5561000 rw-p 00054000 08:05 2490412                    /lib64/libm-2.10.1.so
> 7f09f5561000-7f09f5651000 r-xp 00000000 08:05 2403626                    /usr/lib64/libstdc++.so.6.0.12
> 7f09f5651000-7f09f5851000 ---p 000f0000 08:05 2403626                    /usr/lib64/libstdc++.so.6.0.12
> 7f09f5851000-7f09f5858000 r--p 000f0000 08:05 2403626                    /usr/lib64/libstdc++.so.6.0.12
> 7f09f5858000-7f09f585a000 rw-p 000f7000 08:05 2403626                    /usr/lib64/libstdc++.so.6.0.12
> 7f09f585a000-7f09f586f000 rw-p 00000000 00:00 0
> 7f09f586f000-7f09f588d000 r-xp 00000000 08:05 2490635                    /lib64/ld-2.10.1.so
> 7f09f590c000-7f09f5a4b000 rw-p 00000000 00:00 0
> 7f09f5a89000-7f09f5a8c000 rw-p 00000000 00:00 0
> 7f09f5a8c000-7f09f5a8d000 r--p 0001d000 08:05 2490635                    /lib64/ld-2.10.1.so
> 7f09f5a8d000-7f09f5a8e000 rw-p 0001e000 08:05 2490635                    /lib64/ld-2.10.1.so
> 7fff95f66000-7fff95f7b000 rw-p 00000000 00:00 0                          [stack]
> 7fff95fff000-7fff96000000 r-xp 00000000 00:00 0                          [vdso]
> ffffffffff600000-ffffffffff601000 r-xp 00000000 00:00 0                  [vsyscall]
> Ακυρώθηκε
>



Mail converted by MHonArc 2.6.19+ http://listengine.tuxfamily.org/