Advanced Search

Louis Lamarche

105 Reputation

7 Badges

19 years, 293 days
Hydro-Quebec
Contact Louis Lamarche

MaplePrimes Activity


These are questions asked by Louis Lamarche

Curious behavior the convert( V, VectorC...

Louis Lamarche 105 Maple 2017
vectorcalculus physics
October 18 2017
1 3

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

Copying "2D Output" equations from Mapl...

Louis Lamarche 105 Maple 2017
export msword 2d-output
October 16 2017
1 2

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 

Any way to reclaim NAG (_d01akc) allocat...

Louis Lamarche 105 Maple 15
integral memory nag
December 09 2011
9 3


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 

Problem with tensor product of Bras and ...

Louis Lamarche 105 Maple
package
November 01 2011
0 0

 

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