|
On 10/16/2017 03:17 PM, Nicholas
Jankowski wrote:
Ok, thanks. Obviously this function can do
whatever it wants since there is not a Matlab implementation,
but because we use it as the underlying integrator for those
other functions, they will no longer produce compatible output
for single inputs. We could easily just add a line to those
other functions to recast output as single since there is
already an issingle check. Just didn't know if quadcc behave
like the other quads in this regard
Yes, for the moment I think we should just work around quadcc in the
m-files by casting to single. The tolerances will already be
correctly reduced so, even though the quadcc function is using
doubles, it won't be wasting time calculating unnecessary precision
only to have it discarded.
--Rik
|