Art Kalb

252 Reputation

12 Badges

17 years, 260 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