Defining differential operator that acts like curl

Loading [MathJax]/jax/output/HTML-CSS/jax.js

First time here? Check out the FAQ!

Hello,

As a newbie in using SAGE (experience with Python and Numpy), I was wondering how to define a differential operator that acts like curl (or $\nabla \times \vec F$ ). The curl for a vector field $\vec F=(F_x,F_y, F_z)$ is defined as a determinant $\mathrm{det} (\nabla,\vec F)$ . I want to define such operators that act on $F_x, F_y, F_z$ without calculating the determinant a priori. In other words if I was to multiply $\frac {\partial }{\partial x}$ with $F_y$ , I would expect to get $\frac{\partial F_y}{\partial x}$ EDIT:

class DiffOpp(SageObject):

    def __init__(self, dep_var):
        self.dep_var = dep_var

    def __mul__(self, f):
        return diff(f, self.dep_var)

add a comment

1 Answer

Sort by » oldest newest most voted

answered 10 years ago

calc314
4191 ●22 ●48 ●113

The previous question here links to a worksheet which gives the following code:

var ('x y z')
def divergence(F):
    assert(len(F) == 3)
    return diff(F[0],x) + diff(F[1],y) + diff(F[2],z) 
def curl(F):
    assert(len(F) == 3)
    return vector([diff(F[2],y)-diff(F[1],z), diff(F[0],z)-diff(F[2],x), diff(F[1],x)-diff(F[0],y)])

There is also a ticket here on this issue that has a link to some code as well.

link

Comments

However $\nabla \cdot \vec F$ work in n-dimensions since it's a dot product
I do not want to just calculate the derivative but I want the differential operator to act as a multiplier

For example the hamiltonian operator $\hat H = \frac {\hbar}{-2m} \nabla^2 + V$ can act on a wavefunction $\psi$ as $\hat H | \psi \rangle = - \frac{\hbar}{2m} \nabla^2 \psi + V \psi$

Marios Papachristou ( 10 years ago )

add a comment

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower

Stats

Asked: 10 years ago

Seen: 499 times

Last updated: Dec 30 '14

Defining differential operator that acts like curl savecancel

1 Answer

Comments