[Octave-bug-tracker] [bug #59403] integral2 and quad2d functions do not accurately document form of integrand function F

[David Bateman]

Follow-up Comment #6, bug #59403 (project octave):

Matlab uses a G3,K7 Guass-Konrad rule and Octave a G7,K15 Guass-Konrad rule in
quad2d. So the only difference is the matlab should call the integrand with a
7x7 matrix of scalar arguments and Octave with a 15x15 matrix of scalar
arguments. So matlab is detecting the issue somehow. Perhaps they do something
like


try
   f(ones(2,1),ones(2,1)
catch
   error('ERROR HERE')
endcatch


To detect matrix cases before the rest of quad2d. But if they are then calling
the integrand better not have any side effects

D.



    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?59403>

_______________________________________________
  Message posté via Savannah
  https://savannah.gnu.org/