geomean is poorly implemented
Reported by Michael BAUDIN
geomean-MB.sci (3.12 KB, text/x-csrc)
A robust implementation of the geometric mean.
geomean is poorly implemented :
* the arguments are not checked e.g. the inputs can be negative
* the implementation is not robust against large inputs, producing INF when the result is actually a double
* there is no test
A sample implementation is provided in attachment.