[Getdp] problem to compute magnetic energy using potential vector and current density

[Frederic Trillaud]

Hi,

Following up on my previous emails, I am having trouble to compute the
magnetic energy using an alternative method based on the potential
vector and the current density (see attached *.doc).

The results are very different and could not find my mistake. I am
attaching the code.

Best,

Frederic
-------------- next part --------------
A non-text attachment was scrubbed...
Name: axisymmetryMagneticEnergy.pdf
Type: application/pdf
Size: 51859 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20130521/813b8b0f/attachment.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: definitionOfMagnetostaticProblem.pro
Type: text/x-csrc
Size: 3836 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20130521/813b8b0f/attachment.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: integrationMethod.pro
Type: text/x-csrc
Size: 1687 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20130521/813b8b0f/attachment-0001.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: jacobianMethod.pro
Type: text/x-csrc
Size: 936 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20130521/813b8b0f/attachment-0002.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: methodForMagnetostaticProblem.pro
Type: text/x-csrc
Size: 4302 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20130521/813b8b0f/attachment-0003.pro>