Hello all,
I've been toying with the idea of an extended precision accumulator
for the dot product written in numpy/core/src/multiarraymodule.c
Once the modification is being performed, there is no reason not to
allow the specification of an output array.
The functions that exist now:
Th
like they would
entail modifying the source file "multiarrayobject.c". The existing
code is:
/*NUMPY_API
Numeric.innerproduct(a,v)
*/
static PyObject *
PyArray_InnerProduct(PyObject *op1, PyObject *op2)
{
If I make changes, can I get them inserted into the numpy source base
(af
Hello all,
I'd like accumulate the summation in extended precision, "double"
sum for float inputs, "long double" sum for "double" inputs.
A way of doing this is to add an optional third argument to dot - to
specify the summation type.
I've looked into the code, and the source file multiarr
