November 2009 - Hedge - tiker.net Mailing Lists

Problem running Maxwell 2D example, and question about creating custom operators

by David Powell

Dear hedge developers, I have played around with hedge a bit, and have done a bit of reading on DG methods. It certainly looks like a very promising tool. I recently updated the git repository to the latest changes, because I would like to be able to import gmsh meshes (although it seems that the mesh.reader package needs to be added to the install scripts by hand). Now when I try running the maxwell-2d example, I get the following errors: Traceback (most recent call last): File "maxwell-2d.py", line 142, in <module> main() File "maxwell-2d.py", line 129, in main fields = stepper(fields, t, dt, rhs) File "/home/dap124/src/hedge/hedge/timestep/rk4.py", line 101, in __call__ self.residual = 0*rhs(t, y) File "/home/dap124/src/hedge/hedge/models/em.py", line 252, in rhs j = self.current.volume_interpolant(t, discr)[e_indices] File "/home/dap124/src/hedge/hedge/data.py", line 247, in volume_interpolant return self.gf.volume_interpolant(t, discr) TypeError: volume_interpolant() takes exactly 2 arguments (3 given) Is this a new bug, or have I somehow corrupted my build of hedge (I did repotool start-over to try and ensure that everything gets updated). Are there instructions somewhere on how to create your own operator, or detailed documentation on how the existing operators work? For my problems I would need to modify the Maxwell operator to allow for different (piecewise constant) values of epsilon and mu in different regions (preferably defined using physical groups in gmsh). This would seem to be possible, as this is already done for other equations eg. VariableVelocityStrongWaveOperator. The overall idea of the operator template seems clear enough, but understanding all the specific details of the code is extremely difficult. So far all my attempts just resulted in incomprehensible errors. regards David

10 years