Group Abstract Group Abstract

Message Boards Message Boards

0
|
4.2K Views
|
4 Replies
|
3 Total Likes
View groups...
Share
Share this post:
GROUPS:
Mathematica

Problem with Interpolation

Ranjan Laha
Posted 10 years ago

Hello,

I am having a problem with Interpolation function in Mathematica. It will be great if I can get some help on this issue.

f = Interpolation[{{{10000., 1000.}, 7.092814429178659*^-35}, {{10000., 10000.}, 0.}, {{1.*^6, 1000.}, 1.2388687117821762*^-33}, {{1.^6, 10000.}, 1.2349871178216793^-33}, {{1.*^6, 100000.}, 1.1951273497985842*^-33}, {{1.^6, 1.*^6}, 0.}, {{1.*^8, 1000.}, 8.865974897582733*^-33}, {{1.^8, 10000.}, 8.86574897582677^-33}, {{1.*^8, 100000.}, 8.863489758267155*^-33}, {{1.^8, 1.*^6}, 8.840897247916136^-33}, {{1.*^8, 1.^7}, 8.607630769651254*^-33}, {{1.^8, 1.*^8}, 0.}, {{1.*^10, 1000.}, 4.7849998780149045*^-32}, {{1.^10, 10000.}, 4.784998780148429^-32}, {{1.*^10, 100000.}, 4.784987801442968*^-32}, {{1.^10, 1.*^6}, 4.78487801031734^-32}, {{1.*^10, 1.^7}, 4.783779691857841*^-32}, {{1.^10, 1.*^8}, 4.7727557829322513^-32}, {{1.*^10, 1.^9}, 4.6599274090650744*^-32}, {{1.^10, 1.*^10}, 0.}, {{1.*^12, 1000.}, 2.944829033896313*^-31}, {{1.^12, 10000.}, 2.97448290275792246^-31}, {{1.*^12, 100000.}, 3.`*^-31}}, InterpolationOrder -> 1]

When I try this in Mathematica, the kernel is quitting. Can I get some help with this?

Thanks and best regards, Ranjan Laha.
POSTED BY: Ranjan Laha
4 Replies
Sort By:
Ranjan Laha
Ranjan Laha
Posted 10 years ago

Thank you for the prompt and helpful response.
POSTED BY: Ranjan Laha

Thanks, now I was able to reproduce it in 9.0.1 but the issue appears to be resolved in version 10.0.0 and later.
POSTED BY: Ilian Gachevski
Ranjan Laha
Ranjan Laha
Posted 10 years ago

Thank you very much for the message. The related notebook is very long. I attach a code, with changed inputs for simplicity, which reproduce the crash.

f = Interpolation[ Flatten[ParallelTable[{{10^Enuindex, 10^ Etauindex}, (NIntegrate[ 10^Enuindex yy, {yy, 0, 1. - 10^Etauindex/10^Enuindex}] + NIntegrate[ 10^Enuindex yy, {yy, 0, 1. - 10^Etauindex/10^Enuindex}]) }, {Enuindex, 4., 12., 2.}, {Etauindex, 3., Enuindex, 1.}], 1]]

This is the code that is reproducing the crash. I am using Mathematica 9.0.1.0 in Mac OS X.

Thanks and best regards, Ranjan Laha.
POSTED BY: Ranjan Laha

I cannot reproduce the crash with the input given. Since the data has duplicate abscissas, Interpolation will return unevaluated with an error message.

Which version and operating system are you using, and could you attach a notebook with input in case it was misformatted?
POSTED BY: Ilian Gachevski
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Reply Preview
Attachments
Remove
or Discard