Art Kalb

252 Reputation

12 Badges

18 years, 7 days

MaplePrimes Activity


These are questions asked by Art Kalb

Hi,

I'm trying to figure out how to get the HermitianTranpose vector to distribute over a sum of vectors. Maple seems to want to ignore the linearity of the conjugate operation.

See the example code attached.

Thanks.
 

alp := [a, b, c, d, e]

[a, b, c, d, e]

(1)

u[a] := Vector(5, proc (m) options operator, arrow; u[alp[m], a] end proc)

u[b] := Vector(5, proc (m) options operator, arrow; u[alp[m], b] end proc)

u[c] := Vector(5, proc (m) options operator, arrow; u[alp[m], c] end proc)

u[d] := Vector(5, proc (m) options operator, arrow; u[alp[m], d] end proc)

u[e] := Vector(5, proc (m) options operator, arrow; u[alp[m], e] end proc)

NULL

NULL

NULL

N := r[a]*u[a]+r[b]*u[b]+r[c]*u[c]+r[d]*u[d]+r[e]*u[e]

Vector[column](%id = 36893490335086041316)

(2)

``

LinearAlgebra:-HermitianTranspose(N)

Vector[row](%id = 36893490335105454316)

(3)

NULL


 

Download DistributeHermitian.mw

By default, CodeGeneration generates temporary names during optimize of t####. For compatibility with other code, I would like to have flexibility in choosing the preamble to the temporary variable name ("t" is default). Is there a way to do this?

I don't see anything in the documentation, but it seems like something one should be able to do.

(My target language is Python - although I'm not really using Python, the syntax is compatible).

I am trying to get Maple to recognize that

diff(x^n,x) does not equal n*x^(n-1), but rather 0 if n=0, or n*x^(n-1) otherwise.

This comes up when differentiating an infinite sum (power series). The constant term gets transformed into n/x instead of becoming zero. Maybe this is really a bug/lack of feature in how sum/Sum works.

For example: diff(Sum(x^n, n = 0 .. infinity), x) yields

 

Does anybody have a fix to get the differentiation right (other than expanding some terms of the series before taking the derivative)? 

(edited) Please note I am not trying to get an answer to this specific question. It is just illustrates a simplified example of behavior that leads to other bizarre results. Please see the attached worksheet for more weirdness: Problems_with_summation_Mapleprimes.mw

Hi,

I just upgraded to Maple 2021 and started getting strange results - similar calculations seemed to work in Maple 2020.

I am looking at correlated bivariate normal distributions. When doing a fully symbolic verification of the normalization, I get an incorrect result of infinity. I also get strange results when integrating over a circular region, but the normalization seemed to be a very elementary calculation.

Maybe someone can check this make sure I'm not losing my marbles. I've attached a worksheet.

Thanks.

Bivariate_Gaussian.mw

Hi,

I am looking to have Maple compute the following:

simplify(diff(diff(Sum(epsilon^k*apply(index(f, k), t)/k!, k = 1 .. infinity), t), epsilon, epsilon))

The return value is:

 

Is there any way to get this to simplify appropriately? The problem is with the case k=1, which should just be 0.

What is the best way to expand this if I want to see some of the first few terms?

 

1 2 3 4 5 6 7 Page 2 of 8