Louis Lamarche

80 Reputation

6 Badges

14 years, 334 days

MaplePrimes Activity


These are replies submitted by Louis Lamarche

When right clicking, it takes ~5 seconds for the context menu to appear.

@Karen 

Thank you for yous answer !
However it appears to me that E_[1] <> E_[3].

@spalinowy 

I have just installed Maple 2017 and getting this issue in Windows 7.

Copying the 2D output now only copies the text output to the clipboard, not the graphic output.

I am a bit disapointed here

Thanks for any solution except reinstalling Maple 2015.

@acer 

Sorry guy's I replied to Mapleprime directly instead of posting here.

BK0(x) is a fast version off the BesselK(0,x) 100 times faster, based on the
Ooura's algorythm

aK0:=hfarray(1..16, [...14 digits numbers....]];
bK0:=hfarray(1..112, [...14 digits numbers....]];
cK0:=hfarray(1..135, [...14 digits numbers....]];
dK0:=hfarray(1..40, [...14 digits numbers....]];
BK0:= proc(X)
local x, k::integer, t, y;
#global aK1, bK1, cK1, dK1;x:=evalf(X);
if (x < 0.86) then
  t := x * x;
  y := ((((((aK0[0+1] * t + aK0[1+1]) * t +
      aK0[2+1]) * t + aK0[3+1]) * t + aK0[4+1]) * t +
      aK0[5+1]) * t + aK0[6+1]) * t + aK0[7+1];
  y := ((((((aK0[8+1] * t + aK0[9+1]) * t +
      aK0[10+1]) * t + aK0[11+1]) * t + aK0[12+1]) * t +
      aK0[13+1]) * t + aK0[14+1]) * t + aK0[15+1] - y * ln(x);
elif (x < 4.15) then
...
return y;
end proc;

In my case rho1= 0.323   and a = 0.01785

I will keep the piecewise since integral from 0 to kmin must also be evaluated, not
causing any memory consumpion issue however.

With dz=7000, kmin=0.224399475256414e-3, kmax= 91.5632886852015
number of cycles 102009, The result is -0.109403613457678e-3 and is given
in ~50 seconds on my dual core 2.8 GHz CPU. One evaluation consumes
1.1 Gbytes of memory. gc() does not appear to reclaim any memory.

Maybe there is no way to send this memory to garbage and reclaim it?
That was my question.







1 2 Page 2 of 2