The multiplication of a triangular matrix with a scalar value is wrongly implemented in triangular.jl, lines 155-172
The lines
for i = UpLo == :L ? j:n : 1:j
should say
for i = UpLo == :U ? j:n : 1:j
The method should probably also work on a copy of A instead of manipulating A itself
Regards,
Tom
—
Reply to this email directly, view it on GitHub, or unsubscribe.
You are receiving this because you were mentioned.![]()