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