Dear Sympy list,
I'm trying to evaluate a numerical dispersion relation for a finite
element method we are working on, which requires construction of a 4x4
Hermitian matrix A which is a function of l,k,w, and evaluation of
it's determinant. This gives a quartic polynomial which can be solved
for
On Jun 29, 3:31 pm, Colin wrote:
> Dear Sympy list,
> I'm trying to evaluate a numerical dispersion relation for a finite
> element method we are working on, which requires construction of a 4x4
> Hermitian matrix A which is a function of l,k,w, and evaluation of
> it's determinant. This gives
Dear Sympy list (sorry for deluge!),
I'm getting some unexpected (to me) behaviour from evalf() - please
see the bottom of my email. Is this a) incorrect expectations, b)
incorrect usage, c) known bug, or d) something else?
all the best
--cjc
In [108]: k = Symbol('k')
In [109]: w = Symbol('w')
On Tue, Jun 29, 2010 at 12:31 PM, Colin wrote:
> Dear Sympy list,
> I'm trying to evaluate a numerical dispersion relation for a finite
> element method we are working on, which requires construction of a 4x4
> Hermitian matrix A which is a function of l,k,w, and evaluation of
> it's determinant.
On 29 June, 12:40, Colin wrote:
> Dear Sympy list (sorry for deluge!),
> I'm getting some unexpected (to me) behaviour from evalf() - please
> see the bottom of my email. Is this a) incorrect expectations, b)
> incorrect usage, c) known bug, or d) something else?
Ah sorry, I went to the mailing
On Tue, Jun 29, 2010 at 1:43 PM, Colin wrote:
> On 29 June, 12:40, Colin wrote:
>> Dear Sympy list (sorry for deluge!),
>> I'm getting some unexpected (to me) behaviour from evalf() - please
>> see the bottom of my email. Is this a) incorrect expectations, b)
>> incorrect usage, c) known bug, o
> So it finished in about 11s. How long does it take for you?
About the same, I'm pretty impressed with the performance on this
calculation.
> But I don't know how to simplify this. In my experience, if you need
> to invert some matrices, different methods give different results, I
> think reason
There are different methods for finding the determinant too. See the docstring
of the .det() method.
Aaron Meurer
On Jun 29, 2010, at 7:49 AM, Colin wrote:
>> So it finished in about 11s. How long does it take for you?
>
> About the same, I'm pretty impressed with the performance on this
> cal
I'm making a manual equation solver, where my students have to
manually apply algebra rules in order to simplify and solve problems.
The problem is that, for some subjects, I want that the expression
x1 = Symbol('x1')
b = x1 + x1
remainx x1 + x1 and not became 2*x1...
is it possibile in sympy?
With a little trickery, yes.
A normal evaluate=False does not help in this case, otherwise
Mul(x,x,evaluate=False) might work. Perhaps issue 1966 could fix this.
But in the meantime, realizing that variables with distinct
assumptions are distinct variables you can create two variables with
the sam
10 matches
Mail list logo