[Octave-bug-tracker] [bug #47654] odeset does not allow zero value for tolerance

[Francesco Faccio]

Follow-up Comment #1, bug #47654 (project octave):

Dear Carlo, 

I have modified the file ode_struct_value_check.m and now it is allowed to set
RelTol=0 or AbsTol=0 (the changeset is in the file odestructvaluecheck.diff, I
hope to have written it in the correct way). 

While I was testing the code I discovered this 


warning: division by zero
warning: called from
    AbsRel_Norm at line 41 column 12
    integrate_adaptive at line 133 column 9
    ode23 at line 332 column 16


so I have put an error message in the file AbsRel_Norm.m in order to check if
almost one of the elements of the two vectors AbsTol and RelTol is set to
zero. This changeset is in AbsRelNorm.diff

Should I write also some tests in ode23 and ode45?

Best regards,

Francesco Faccio

 

(file #36872, file #36873)
    _______________________________________________________

Additional Item Attachment:

File name: odestructvaluecheck.diff       Size:1 KB
File name: AbsRelNorm.diff                Size:0 KB


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?47654>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/