Wrong outcome from kron() with small numbers

After spending a day, I found a bug in kron() unless I miss something.

From the following code,

#include "blaze/Blaze.h"
#include <iostream>

using namespace blaze;

int main()
{
    DynamicVector<double> P1{0.073512780078877, 0.564448879998923, 0.350576679157777, 0.011448881656057, 0.000012778790294, 0.000000000316875, 0.000000000001194, 0.000000000000004, 0.000000000000000};

    DynamicVector<double> P2{1.0};

    std::cout << std::fixed << std::setprecision(15);
    std::cout << P1 << "\n";
    std::cout << P2 << "\n";
    std::cout << kron(P1, P2) << "\n";
};

I got the following outcome.

Please see that the last four numbers become zero.

( 0.073512780078877 )
( 0.564448879998923 )
( 0.350576679157777 )
( 0.011448881656057 )
( 0.000012778790294 )
( 0.000000000316875 )
( 0.000000000001194 )
( 0.000000000000004 )
( 0.000000000000000 )

( 1.000000000000000 )

( 0.073512780078877 )
( 0.564448879998923 )
( 0.350576679157777 )
( 0.011448881656057 )
( 0.000012778790294 )
( 0.000000000000000 )
( 0.000000000000000 )
( 0.000000000000000 )
( 0.000000000000000 )

‌

Comments (5)