sasomao

602 Reputation

7 Badges

12 years, 85 days
Phd.
Paris, France

MaplePrimes Activity


These are replies submitted by sasomao

Hi

I discovered that it works well if I break the command in two parts:

storage:= simplify( A B* +  A* B) assuming f::positive:
evalf(IntegrationTools:-Expand(2*Int( storage, f=f_low..f_high))) :
 
only I really can't understand why Maple doesn't simplify automatically that expression, if it knows that f is positive (as the two extreme are real positive float numbers)!!!
And, if I remember well, an expression like A B* + A* B should be real whitout regards on what the scalars functions A and B are (reals or complex)
I feel there is something important I don't understand about this Maple's behaviour.

Thank for any remarks.

S.V.

thanks acer and Doug

it works!

S.

thanks acer and Doug

it works!

S.

Hi

I tried  

nohup /usr/bin/maple < /home/Salvatore/Desktop/aaa.mw > aaa.txt
 

with aaa.mw being :

eval( 0.475632267551*10^(-47)/Enn^2,
Enn = rho/sqrt( 0.187440162848*10^46)):
sqrt(%); evalf(eval(%,rho=10));

and aaa.txt contains

    |\^/|     Maple 13 (IBM INTEL LINUX)
._|\|   |/|_. Copyright (c) Maplesoft, a division of Waterloo Maple Inc. 2009
 \  MAPLE  /  All rights reserved. Maple is a trademark of
 <____ ____>  Waterloo Maple Inc.
      |       Type ? for help.
> <?xml version="1.0" encoding="UTF-8"?>
on line 1, syntax error, character `?` unexpected:
<?xml version="1.0" encoding="UTF-8"?>
 ^
> quit
memory used=0.6MB, alloc=0.8MB, time=0.02

Something seems to be wrong!

S.V.

Hi

I tried  

nohup /usr/bin/maple < /home/Salvatore/Desktop/aaa.mw > aaa.txt
 

with aaa.mw being :

eval( 0.475632267551*10^(-47)/Enn^2,
Enn = rho/sqrt( 0.187440162848*10^46)):
sqrt(%); evalf(eval(%,rho=10));

and aaa.txt contains

    |\^/|     Maple 13 (IBM INTEL LINUX)
._|\|   |/|_. Copyright (c) Maplesoft, a division of Waterloo Maple Inc. 2009
 \  MAPLE  /  All rights reserved. Maple is a trademark of
 <____ ____>  Waterloo Maple Inc.
      |       Type ? for help.
> <?xml version="1.0" encoding="UTF-8"?>
on line 1, syntax error, character `?` unexpected:
<?xml version="1.0" encoding="UTF-8"?>
 ^
> quit
memory used=0.6MB, alloc=0.8MB, time=0.02

Something seems to be wrong!

S.V.

Thanks guys.

This reply of acer, in particular, is very intresting.

Just two questions:

1) what if I want to run this command line (nohup maple <input> output & ) with a niceness 19?
2) With this command, all the lines in the maple input code that end with a ; will have their results printend in output-file?

Thanks.

S.
 

Thanks guys.

This reply of acer, in particular, is very intresting.

Just two questions:

1) what if I want to run this command line (nohup maple <input> output & ) with a niceness 19?
2) With this command, all the lines in the maple input code that end with a ; will have their results printend in output-file?

Thanks.

S.
 

thank everybody!

It's good to know about IntegrationTools, I'll try it. Before acer replied I made this procedure to hande the situation

integral_cutter:=proc(A)
global parameters,f_low,f_high;
local number,i,temp,storage,final,B;
if (has(A,Int)) then

B:=op(1,select(has,A,Int));
number:=nops(expand(B));
final:=0;
storage:=0;

for i from 1 to numer do
temp:=op(i,expand(B));
if has(temp,parameters) then
storage:=select(has,temp,parameters)*Int(remove(has,temp,parameters),f=f_low..f_high);
else
storage:=Int(op(i,expand(B)),f=f_low..f_high);
end if;
final:=final+storage;
end do;

else
return A;
end if
end proc:
 

And it seems to work.

You were really nice. I spent all the week-end thinking how to solve this problem!

thanks!

S.V.

thank everybody!

It's good to know about IntegrationTools, I'll try it. Before acer replied I made this procedure to hande the situation

integral_cutter:=proc(A)
global parameters,f_low,f_high;
local number,i,temp,storage,final,B;
if (has(A,Int)) then

B:=op(1,select(has,A,Int));
number:=nops(expand(B));
final:=0;
storage:=0;

for i from 1 to numer do
temp:=op(i,expand(B));
if has(temp,parameters) then
storage:=select(has,temp,parameters)*Int(remove(has,temp,parameters),f=f_low..f_high);
else
storage:=Int(op(i,expand(B)),f=f_low..f_high);
end if;
final:=final+storage;
end do;

else
return A;
end if
end proc:
 

And it seems to work.

You were really nice. I spent all the week-end thinking how to solve this problem!

thanks!

S.V.

thanks!

 

it works as you said.
Maybe it woud be useful to write this hint in the overload help page, it doesn't seem evident to me that one shoud use the $ marker.

 

S.

thanks!

 

it works as you said.
Maybe it woud be useful to write this hint in the overload help page, it doesn't seem evident to me that one shoud use the $ marker.

 

S.

Hi

 

but a full evaluation should'n have removed the quotes, like it happens in what follows:

 

a:=2;

eval('a');

2

 

The full evaluation performed in the head of a procedure is not of the same kind of that of "eval"?

 

Tkyou

 

S.V.

Hi

 

but a full evaluation should'n have removed the quotes, like it happens in what follows:

 

a:=2;

eval('a');

2

 

The full evaluation performed in the head of a procedure is not of the same kind of that of "eval"?

 

Tkyou

 

S.V.

In all cases, what I want is a formal result for the riemann tensor. That is, I would know the ijkm component without fix the actual value for i j k m (after all the possible simplification). I don't understand if Maple can do it. Can you make clearer this point? Tk Salvatore
In all cases, what I want is a formal result for the riemann tensor. That is, I would know the ijkm component without fix the actual value for i j k m (after all the possible simplification). I don't understand if Maple can do it. Can you make clearer this point? Tk Salvatore
First 7 8 9 10 Page 9 of 10