Louis Lamarche

105 Reputation

7 Badges

16 years, 45 days
Hydro-Quebec

MaplePrimes Activity


These are questions asked by Louis Lamarche

Converting a VectorCalculus vector to a PhysicsVectors vector and reconverting it back
to VectorCalculus vector do not return the same result !
is this a bug ? Any clue ?

with(plots):
with(VectorCalculus):
with(Physics[Vectors]):
with(inttrans):
Setup(mathematicalnotation = true):

SetCoordinates('spherical'[r, phi, theta]);
E_[1]:=Vector(<r,phi,theta>);
'VectorCalculus(E_[1])'=type(E_,Vector);
'PhysicsVectors(E_[1])'=type(E_,PhysicsVectors);

Ev_[1]:=convert(E_[1],PhysicsVectors);
'VectorCalculus(Ev_[1])'=type(Ev_,Vector);
'PhysicsVectors(Ev_[1])'=type(Ev_,PhysicsVectors);

E_[2]:=convert(Ev_[1],VectorCalculus);
'E_[2]'<>'E_[1] !';

spherical[r, phi, theta]

 

`#mover(mi("E"),mo("&rarr;"))`[1] := Vector(3, {(1) = r, (2) = phi, (3) = theta}, attributes = [coords = spherical[r, phi, theta]])

 

VectorCalculus(E_[1]) = true

 

PhysicsVectors(E_[1]) = false

 

_phi*phi+_r*r+_theta*theta

 

VectorCalculus(Ev_[1]) = false

 

PhysicsVectors(Ev_[1]) = true

 

`#mover(mi("E"),mo("&rarr;"))`[2] := Vector[row](3, {(1) = r, (2) = theta, (3) = phi}, attributes = [vectorfield, coords = spherical[r, phi, theta]])

 

E_[2] <> factorial(E_[1])

(1)

 


 

Download BugConversion.mw

In previous version of Maple (e.g. Maple 2015) I was able to copy
"2D Output" equations directly from a Maple Worksheet to MS Word.

These equations could be edited for fine element placement. In the 2017.3 version
it is does not seem to be possible anymore. When I copy the "2D Output" the image
is not copied only the text.

Is this a feature that was removed or a configuration problem ?

Thank you

LL 


I am evaluating a difficult integral of an oscillating function.

f:=k->evalhf(piecewise(k=0,1,k>0,cos(k*(dz))*BK0(k*rho1)/BK0(k*a)));

evalf(Int(f, kmin..kmax, method = _d01akc,epsilon=eps,maxintervals=1000000))

with dz>0 and possibly very large
kmin is Pi/(2*dz)
the number of cycles is nc = ceil((1/2)*kl*dz/Pi) with kl = 91.563
kmax=  kmin+nc*2*Pi/dz

Each evaluation of this integral using a call to the NAG routine (_d01akc ) eats up 

 

The Physics package appears broken for the tensor product of Bras and Kets.
It is not possible to describe new operators creatde by the tensorial
product of operators on two Hilbert spaces. ie H12 = H1 X H2

Download BraKetTensorProduc.mw

Am I missing something here.

Thank you for your help

1 2 Page 2 of 2