MathGroup Archive 2009

[Date Index] [Thread Index] [Author Index]

Search the Archive

Interpolation -> segfault?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg98128] Interpolation -> segfault?
  • From: oshaughn <oshaughn at gravity.psu.edu>
  • Date: Tue, 31 Mar 2009 04:16:24 -0500 (EST)

Hello,

I generated and saved the following interpolating function (and
wrapper).  Whenever I evaluate it, I immediately get a segfault in V7
and V6.

Could someone let me know why this is, why the code that generated it
doesn't segfault, and what I should do to avoid segfault issues in the
future?

fn = If[#1 > 0.1 && #1 <
      10^2, (1/(#1 Log[10])) InterpolatingFunction[{{-1., 2.}}, {3, 3,
         0, {136}, {4}, 0, 0, 0,
        0}, {{-1., -0.75, -0.5, -0.25, -0.125, 0., 0.125, 0.15625,
         0.1875, 0.21875, 0.234375, 0.25, 0.265625, 0.273438, 0.28125,
          0.289063, 0.296875, 0.304688, 0.3125, 0.320313, 0.328125,
         0.335938, 0.34375, 0.351563, 0.359375, 0.367188, 0.375,
         0.382813, 0.390625, 0.398438, 0.40625, 0.414063, 0.421875,
         0.429688, 0.4375, 0.445313, 0.453125, 0.460938, 0.46875,
         0.476563, 0.484375, 0.492188, 0.5, 0.507813, 0.515625,
         0.523438, 0.53125, 0.539063, 0.546875, 0.554688, 0.5625,
         0.570313, 0.578125, 0.585938, 0.59375, 0.601563, 0.609375,
         0.617188, 0.625, 0.632813, 0.640625, 0.648438, 0.65625,
         0.664063, 0.671875, 0.679688, 0.6875, 0.695313, 0.703125,
         0.710938, 0.71875, 0.726563, 0.734375, 0.742188, 0.75,
         0.757813, 0.765625, 0.773438, 0.78125, 0.789063, 0.796875,
         0.804688, 0.8125, 0.820313, 0.828125, 0.835938, 0.84375,
         0.851563, 0.859375, 0.867188, 0.875, 0.882813, 0.890625,
         0.898438, 0.90625, 0.914063, 0.921875, 0.929688, 0.9375,
         0.945313, 0.953125, 0.960938, 0.96875, 0.976563, 0.984375,
         0.992188, 1., 1.00781, 1.01563, 1.02344, 1.03125, 1.03906,
         1.04688, 1.05469, 1.0625, 1.07031, 1.07813, 1.08594, 1.09375,
          1.10156, 1.10938, 1.11719, 1.125, 1.13281, 1.14063, 1.14844,
          1.15625, 1.17188, 1.1875, 1.20313, 1.21875, 1.23438, 1.25,
         1.5, 1.75, 2.}}, {Developer`PackedArrayForm, {0, 1, 2, 3, 4,
         5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20,
         21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35,
         36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50,
         51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65,
         66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80,
         81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95,
         96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108,
         109, 110, 111, 112, 113, 114, 115, 116, 117, 118, 119, 120,
         121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132,
         133, 134, 135, 136}, {0., 0., 1.839556450000000*10^-316,
         1.21549*10^-159, 2.31335*10^-101, 1.67253*10^-56,
         4.5936*10^-25, 2.66022*10^-19, 2.23311*10^-14,
         2.71726*10^-10, 1.4528*10^-8, 4.79291*10^-7, 9.75747*10^-6,
         0.0000367374, 0.000122606, 0.000362738, 0.000951569,
         0.00221421, 0.00457332, 0.00839564, 0.0137329, 0.0201077,
         0.0265815, 0.0322166, 0.0367286, 0.0408717, 0.0462177,
         0.0544092, 0.0663, 0.0813791, 0.0976722, 0.112228, 0.122235,
         0.126523, 0.126732, 0.127295, 0.133834, 0.150431, 0.177017,
         0.208434, 0.236308, 0.253376, 0.257998, 0.255774, 0.256824,
         0.270362, 0.300085, 0.342984, 0.391402, 0.436275, 0.469696,
         0.486402, 0.484727, 0.467353, 0.441479, 0.41775, 0.407822,
         0.421313, 0.4633, 0.533418, 0.626979, 0.737778, 0.861279,
         0.996289, 1.1438, 1.30339, 1.46955, 1.63063, 1.77232, 1.8848,
          1.96966, 2.04154, 2.1219, 2.2277, 2.36187, 2.51091, 2.65028,
          2.7539, 2.80324, 2.79347, 2.7351, 2.65062, 2.56727, 2.50909,
          2.49206, 2.52399, 2.60708, 2.73955, 2.91432, 3.11604,
         3.32016, 3.49681, 3.61939, 3.67374, 3.66173, 3.59611,
         3.48923, 3.34377, 3.15256, 2.90785, 2.61315, 2.28878,
         1.96732, 1.68224, 1.45678, 1.29843, 1.20019, 1.14535,
         1.11267, 1.08029, 1.02917, 0.946603, 0.829484, 0.685201,
         0.529164, 0.379515, 0.25137, 0.15301, 0.0852144, 0.0432423,
         0.0199185, 0.00829953, 0.00311858, 0.00105388, 0.000319547,
         0.0000867615, 0.0000210592, 8.81696*10^-7, 2.32524*10^-8,
         3.8403*10^-10, 3.95627*10^-12, 2.53536*10^-14,
         1.00878*10^-16, 1.67729*10^-83, 7.17379*10^-204,
         0.}}, {Automatic}][Log[10, #1]], 0] &;


  • Prev by Date: Re: Exporting a discrete Manipulate[] sequence?
  • Next by Date: Re: simpler way to get a particular banded matrix?
  • Previous by thread: Graphics and DisplayFunction