Currently cannot do something like B["ij"]=A["ij"]B["ij"] with sparse A, B, consider instead using summation with a special operator, i.e Transform<>([](double a, double b){ b=a; })(A,B);. However, there is no such simple workaround for computations like C["ijk"]=A["ijl"]*B["ljk"], so this is a high priority enhancement. For now please do this as dense and expect a fix for the sparse case soon.
Currently cannot do something like B["ij"]=A["ij"]B["ij"] with sparse A, B, consider instead using summation with a special operator, i.e Transform<>([](double a, double b){ b=a; })(A,B);. However, there is no such simple workaround for computations like C["ijk"]=A["ijl"]*B["ljk"], so this is a high priority enhancement. For now please do this as dense and expect a fix for the sparse case soon.