An Algebraic Programming Style for Numerical Software and its Optimization