Ok

Here I copy a modified version of my input file. I launch it via ssh as was suggested to me in the forum

nohup maple <input.mpl> out.out &

and then I renice it +19. The system is a Linux Fedora 11 64bit, Maple 13.

assume(eta,'positive');

assume(phi,'positive');

assume(t,'real');

assume(Enn,'positive');

assume(M,'positive');

assume(f,'positive');

assume(Emm,'positive');

with(LinearAlgebra):

parameters:= [Enn,t,phi,Emm,eta]:

param_dim:= nops(parameters):

f_low:=50.0:

f_lso:= (6^(3/2)*Pi*M_tot)^(-1):

f_high:=f_lso:

M_test:=7*10^(-6):

m1:= M_test*1.4 : #########change it

m2:= M_test*1.4 : ######### change it

M_tot:=m1+m2:

eta_eff:=(m1*m2)/M_tot^2:

v_lso:=(Pi*M_tot*f_lso)^(1/3):

variables:=[f]:

variab_dim:=nops(variables):

Digits:=15:

p:= proc()

local Psi,i,psi,alpha_par,v,theta,euler_gamma,lambda;

global Enn,PNorder,param_dim,f,eta,M,M_tot,f_lso,v_lso,Emm;

M:=Emm*eta^(-3/5);

v:=(Pi*M*f)^(1/3);

theta:=-1.28;

lambda:=-0.6451;

euler_gamma:=0.577215664901532860606512090082:

alpha_par:=[

1

,0

,3*eta

,-16*Pi];

Psi:= (2*Pi*f*t-phi-Pi/4+ 3/(128*eta*v^5)*add(alpha_par[i+1]*v^i, i=0..3));

Enn*f^(-7/6)*exp(I*Psi):

end proc:

Digits:=10:

h1:=Array(1..param_dim, i->

'eval'(

(diff(p(), parameters[i]))

,[Emm=(M_tot*(eta_eff)^(3/5)),eta=eta_eff] )

):

scalar_pr:= proc(A,B)

global variables, parameters,f_low, f_high,S_h,M_tot,eta_eff;

local num;

num:=

simplify(

A*conjugate(B)+conjugate(A)*B

)

;

evalf(

IntegrationTools:-Expand(

2*Int(num/S_h(f), f=f_low..f_high)

)

)

;

#simplify(%);

end proc:

Digits:=10:

AMatrix:=proc()

global low,f_low,f_high,param_dim,M_tot,eta_eff:

local gg,a_1,b_1:

gg:=(a_1,b_1)->

scalar_pr(h1[a_1],h1[b_1])

;

low:= Matrix(param_dim,param_dim, shape=symmetric,gg);

end proc:

AMatrix():

save low, "./low.txt":

Inverse:=proc()

global up,M_tot,param_dim:

local a_3,b_3,gg;

gg:=(a_3,b_3)-> (MatrixInverse(low))[a_3,b_3];

up:= Matrix(param_dim,param_dim,shape=symmetric,gg);

end proc:

Inverse():

h2:=Array(1..param_dim, 1..param_dim, (i,j)->

'eval'(

(diff(p(), parameters[i],parameters[j]))

,[Emm=(M_tot*(eta_eff)^(3/5)),eta=eta_eff])

):

h3:=Array(1..param_dim, 1..param_dim,1..param_dim, (i,j,k)->

'eval'(

(diff(p(), parameters[i],parameters[j],parameters[k]))

,[Emm=(M_tot*(eta_eff)^(3/5)),eta=eta_eff])

):

try:=proc(i)

local temp,MMM,NNN,PPP,QQQ,TTT,ZZZ:

global up,low,V,param_dim;

return

(

+ add(add(add(add(

up(i,MMM)*up(i,NNN)*up(PPP,QQQ)*

(

+ V4[NNN,MMM,PPP,QQQ]

+ 3*scalar_pr(h2[NNN,QQQ],h2[PPP,MMM])

+ 2*V31[NNN,MMM,PPP,QQQ]

+ V31[MMM,PPP,QQQ,NNN]

)

,MMM=1..param_dim),NNN=1..param_dim),PPP=1..param_dim),QQQ=1..param_dim)

+

add(add(add(add(add(add(

up(i,MMM)*up(i,NNN)*up(PPP,ZZZ)*up(QQQ,TTT)*

(

+ V3[NNN,PPP,MMM]*V3[QQQ,ZZZ,TTT]

+ 5/2*(-scalar_pr(h2[NNN,PPP],h1[QQQ])*V3[MMM,ZZZ,TTT]-scalar_pr(h2[QQQ,PPP],h1[NNN])*V3[MMM,ZZZ,TTT])

+5/2*( 2*scalar_pr(h2[NNN,QQQ],h1[PPP])*scalar_pr(h2[MMM,ZZZ],h1[TTT]) + scalar_pr(h2[NNN,QQQ],h1[PPP])*scalar_pr(h2[MMM,TTT],h1[ZZZ]) +scalar_pr(h2[NNN,QQQ],h1[PPP])*scalar_pr(h2[TTT,ZZZ],h1[MMM]) )

+2*V21[QQQ,ZZZ,NNN]*V3[MMM,PPP,TTT]

+2* V3[NNN,MMM,TTT]*V21[QQQ,PPP,ZZZ]

+6* V21[NNN,TTT,ZZZ]*V3[MMM,PPP,QQQ]

+V21[MMM,TTT,NNN]*V3[PPP,QQQ,ZZZ]

+2*V21[NNN,QQQ,ZZZ]*V21[PPP,TTT,MMM]

+2*V21[NNN,QQQ,MMM]*V21[PPP,TTT,ZZZ]

)

,MMM=1..param_dim),NNN=1..param_dim),PPP=1..param_dim),QQQ=1..param_dim),TTT=1..param_dim),ZZZ=1..param_dim)

)

end proc:

prova2:= try(2):

prova2:= simplify(Re(try2));#assuming Enn::positive;

save prova2, "./try2.txt":

where the v's are array of this kind

v3:=proc()

global h1,h2,h3,param_dim, V3;

V3:= Array(1..param_dim,1..param_dim,1..param_dim, (ii,jj,kk)->

- scalar_pr(h2[ii,jj],h1[kk])

):

end proc:

v3():

and S_h is a rel valued function. (I made changes with respect to the true code, for privacy reasons, so if something is unclear please ask)

I launched the imput 5 times, in the same machine, and the file try2.txt contained always a different value:

prova2 := .2050331085e-93/Enn^4;

prova2 := .7892899913e-93/Enn^4;

prova2 := .3927438368e-93/Enn^4;

prova2 := -.2472605212e-93/Enn^4;

prova2 := .6520808383e-93/Enn^4;

prova2 := -.1950615287e-93/Enn^4;

This means that my work of the last three weeks is garbage. I'd love at least to discover why.

Thanks

S.