Comparison of the results of simulating OOMMF problem with 'double' and 'float' for different control point specifications


The problem that was taken into contention for simulating had the following specifications:

1) The geometry considered was the following:

2) Initial Magnetization: initmag2.ovf

The problem which was simulated for obtaining the initial magnetization for this problem, had the following specifications:

Starting orienation: Q=89°, F=90°
Applied field: Bx=100 mT to 0 mT in steps of 25 mT,200 iterations at each step
Relaxation magnetization at Bx=0 mT, with max torque=10-4
Precision: long double (default setting in 'app/mmsolve/fft.h'
No time control point specification was given.
The OMF file was converted to an OVF file using 'tclsh oommf.tcl avf2ovf <initmag2.omf >initmag2.ovf'

(3) Field Range: 0.0 0.0 0.0 0.0 0.0 0.0 1 -torque 1e-5 -time 1e-7

This implies that 0 field was applied and the simulation was ended when either the torque at a step was less than 1E-5 or when the elapsed field time for the current field step reaches 1E-7 seconds.

These specification were changed for the different simulations

Calculation of Metric

The comparison of 2 simulations are done on the basis of their OVF files corresponding to a particular step. The OVF file has the magnetization vector values at all the points considered in the geometry of the sample.

The format of the OVF file can be understood from the following sample OVF file:


The metric which we have considered is the Euclidean Norm which is calculated in the following way:

Let one of the vectors be given by *X*=[x 1, x 2,...x M]

Let the other vector be given by *Y*=[y 1, y 2,...y M]

The Euclidean Norm is given by: [(x 1-y 1) 2+(x 2-y 2) 2+...(X M-y M) 2] 1/2


No. Data Type M i (x,y) m x h _ Iterations max( m i - m 0 ) Time spent simulating Figure Comments
1 Double M 0 1e-4 Reference initmag2.ovf
2 Double M 1 10000 1.62e-2 398.05 seconds Figure 1 The largest error is seen approximately at the coordinates (0.0, 1.0) and (1.0, 1.0), where, from the field plot we can see that,a discontinuity arises.
3 Float M 2 10000 4..5e-3 372.54 seconds Figure 2

Figure 1: Magnetization and Euclidian Norm after simulating with double
with different time control points

Figure 2: Magnetization and Euclidean norm after simulating with double and float
with same control specifications.

geometry.png (7.08 KB) Anil Prabhakar, 07/22/2013 11:55 AM

sampleOVF.txt Magnifier (17.3 KB) Anil Prabhakar, 07/22/2013 11:55 AM

float.png (310 KB) Anil Prabhakar, 07/22/2013 11:55 AM

double.png (308 KB) Anil Prabhakar, 07/22/2013 11:55 AM

Redmine Appliance - Powered by TurnKey Linux