|Message: Space varying Magnetic Field Problem||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
Hi folks, I've constructed a model that implements a 3D space varying magnetic field in a particular volume. The field is read in from a file and is a set of values of Bx, By, and Bz. The field values are defined on a 1 cm grid. I wrote the GetFieldValue() function to return the pertinent field components as needed by the stepper. My problem is that I get vastly different steps by changing the field value slightly.
The function GetFieldValue(const G4double point, double *field) that is required from me is written so that each time it is called my code does a tri-linear interpolation on my 3D grid to find the field at the input point from the stepper point,point,point. I'm thinking that since this interpolation is happening on a coarse 1cm x 1cm x 1cm grid then the field values are not smooth enough to give a good smooth track.
I'd like some advice as to the best way to proceed. Should I increase the grid by a factor of 10 so that it is defined every 1mm instead of 1cm or should I use a different stepper. Currently I'm simply using
G4ChordFinder *chordFinder = new G4ChordFinder(magField, stepMinimum);
which I believe uses the RK4 stepper. If I need to change the stepper how do I do that?
|Inline Depth:||Outline Depth:||Add message:|