<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"><html xmlns="http://www.w3.org/1999/xhtml"><head>

    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />

 </head><body style="">

  <div>

   Dear Getdp-users/developers,

   <br />I am afraid I do not completely understand the meaning of VolRectShell.

   <br />I assume it transforms a shell (constructed as set difference of two co-centric axi-parallel cubes, the smaller one contained in the interior of the larger one) into the (infinitely large) complement of a cube.

   <br />

   <br />The only hint I found so far for the rectangular shell transformation is from the paper of Abert et. al.: "magnum.fe: A micromagnetic finite-element simulation code based on FEniCS" (http://arxiv.org/abs/1301.3061).

   <br />

   <br />I have attached a geometry file CuboidCyl.geo and a problem description CuboidCyl.pro. Both of them include the attached parameter file CuboidCyl.par. Note, that I am exploiting symmetry in all three space directions.

   <br />

   <br />In GuboidCyl.pro I have tried to apply VolRectShell without luck.

   <br />When I run getdp I get the error

   <br />Error   : GetDP - Bad parameters for transformation Jacobian: 0.0163227 not in [0.023,0.047571]

   <br />

   <br />Perhaps my misunderstanding of the radius parameters of VolRectShell is the origin of this error. I already tried several variants in this respect. One of those variants is implemented in CuboidCyl.pro.

   <br />

   <br />Could you please tell me what I am doing wrong?

   <br />

   <br />The problem description works (of course with a worse approximation) if I replace VolRectShell by Vol in domInf.

   <br />

   <br />With best regards,

   <br />Tobias Zawada

   <br />

   <br />

  </div>

</body></html>