Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Intel Community
- Software
- Software Development SDKs and Libraries
- Intel® oneAPI Math Kernel Library
- Question regarding step size control in ?jacobix

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

Jim_G_1

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

10-02-2013
12:07 PM

24 Views

Question regarding step size control in ?jacobix

Hi Everyone,

I am using djacobix to differentiate a function which can only take positve inputs. As my program proceeds jacobix the point of interest approaches zero and eventually jacobix tries to evaluate the function with negative arguments.

It seems to me that at some point the way jacobix is using the "eps" parameter changes; for large point of interest x the step size used in (what i assume is) finite differencing is

dx = (plus and minus) eps*x

which remains positive as x approaches zero. However at some point jacobix starts using

dx = (plus and minus) eps

resulting in negative arguments to the target function.

Can someone explain this behaviour so I can suppress it. I'm happy with a cluge that modifies the value of eps to avoid negative values if that is necessary - in which case knowing when jacobix changes behaviour is enough. If it's complicated I will probably write my own routine.

Cheers!

Jim

Link Copied

1 Reply

mecej4

Black Belt

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

10-03-2013
08:58 PM

24 Views

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

For more complete information about compiler optimizations, see our Optimization Notice.