minor update of sum_seq

When input array is not empty, zero(T) is not required to be defined.
Provides specialized Bidiagonal linear solver and eigensystem solver

Addresses #3688

- New tests of Float16 and BigFloat
- Miscellaneous improvements to Bidiagonal tests
- Clean up bidiagonal methods, particularly to reduct overtyping
- Adds BigFloat(x::Float16)