Computing Magnetic and Potential fields from a Time Dependent Dipole

Example i41: Computing Magnetic and Potential fields from a Time Dependent Dipole

Computes magnetic fields at squid sensors and potential fields at torso surface electrodes from a time dependent moving dipole.

The comfile run by this example is as follows:

$EXPORT=1;
$DIPOLE=1;

$OUT="output";
$FORMAT="binary";

$MAGOFFSET=20000;
$POTOFFSET=20000;

if(!defined $example) 
{  
    $example='.' ;
    $OUT="output";
} else{
    $OUT=".";
}


fem def para;r;$example/torso_dipole
fem def coor 3,1;
fem def base;r;$example/bicubic_bilinear;

fem def node;r;$example/skin;
fem def elem;r;$example/skin_lin;

if($EXPORT)
{
    fem export node;$OUT/torso as torso;
    fem export elem;$OUT/skin  as skin;
}

#
# Group nodes for boundary conditions
#
fem group node 29      as FIXEDNODE;
fem group node 1,290   as SURFENDS;
fem group node 2..289  as SURFNORM;
$COMPNODES="1..290"

fem def equa;r;$example/laplace_nr1 reg all;
fem def mate;r;$example/mate_nr1    reg all;
fem def init;r;$example/noflux_nr1  reg all;
fem def sour;r;$example/dipole_known reg $DIPOLE;

if($EXPORT)
{
    fem exp sour;$OUT/dipole_known as dipole_known;
}

fem def solv;r;$example/static_nr1 reg all;


#
# Electrodes and sensors
#
fem def data;r;$example/electrodes_mag field num_field 3;
fem def xi;c close;
fem exp data;$OUT/electrodes_mag as electrodes_mag offset $MAGOFFSET;

Additional testing commands:


$TSTART="0"
$TEND="3.1"
$TSTEP=1
$FREQ=1

for ($TIME=$TSTART;$TIME<=$TEND;$TIME+=$TSTEP)
{
  printf("*** Dipole solve at time= $TIME\n")
  printf("*********************************\n")
  fem up sour reg $DIPOLE time $TIME
  fem solve at $TIME
  fem up phi tstart $TIME tend $TIME from data phi #data electrodes

  printf("\n\n")
  printf("*** Magnetic solve at time= $TIME\n")
  printf("*********************************\n")
  fem eval solu mag from MFI nss 1 time $TIME frequency $FREQ 
  fem up data field from MFI nss 1 time $TIME frequency $FREQ
  fem export data;$OUT/mfi_known_$TIME as mfi_known field_num 1..3 field_name mag offset $MAGOFFSET
} 

#
# store matrices of potential/magnetic fields
#
fem write matr;$OUT/mfi matri MFI ascii
fem write matr;$OUT/phi matri PHI ascii # data points

Files used by this example are:

Name                     Modified     Size

example_i41.com          05-Dec-2002  1.1k
bicubic_bilinear.ipbase  05-Dec-2002  7.8k
dipole_known.ipsour      05-Dec-2002  1.2k
electrodes_mag.ipdata    05-Dec-2002  522 
laplace_nr1.ipequa       26-May-2003  1.3k
mate_nr1.ipmate          05-Dec-2002  859 
noflux_nr1.ipinit        05-Dec-2002  575 
skin.ipnode              05-Dec-2002  258k
skin_lin.ipelem          05-Dec-2002  118k
static_nr1.ipsolv        13-Apr-2007  1.2k
test_output.com          14-Feb-2003  801 
torso_dipole.ippara      21-Feb-2003  5.9k

Download the entire example:

Name                      Modified     Size

examples_i_i4_i41.tar.gz  18-Jan-2011   73k

Testing status by version:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:03:57 2016	2
cm-debug	Success	Sat Mar 5 00:19:25 2016	10
mips-irix
cm	Success	Sun Aug 19 01:57:40 2007	28
cm-debug	Success	Wed Aug 15 02:32:54 2007	142
cm-debug-clear-malloc	Success	Sat Aug 18 02:53:35 2007	161
cm-debug-clear-malloc7	Success	Mon Aug 20 02:50:47 2007	161
cm64	Success	Sun Aug 19 02:00:21 2007	30
cm64-debug	Success	Tue Aug 21 02:36:01 2007	154
cm64-debug-clear-malloc	Success	Thu Apr 1 12:36:34 2004	60
rs6000-aix
cm	Success	Wed Mar 4 01:10:06 2009	3
cm-debug	Success	Mon Mar 2 01:15:04 2009	32
cm64	Success	Wed Mar 4 01:10:08 2009	3
cm64-debug	Success	Tue Mar 3 01:21:57 2009	35
x86_64-linux
cm	Success	Sun Mar 6 00:01:12 2016	1
cm-debug	Success	Sat Mar 5 00:02:03 2016	5

Testing status by file:

cmiss_test.log

i686-linux
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
mips-irix
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
Success	cm-debug-clear-malloc:	cmiss_test.log.retain.
Success	cm-debug-clear-malloc7:	cmiss_test.log.retain.
Success	cm64:	cmiss_test.log.retain.
Success	cm64-debug:	cmiss_test.log.retain.
Success	cm64-debug-clear-malloc:	cmiss_test.log.retain.
rs6000-aix
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
Success	cm64:	cmiss_test.log.retain.
Success	cm64-debug:	cmiss_test.log.retain.
x86_64-linux
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.

cmiss_test.out

i686-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
mips-irix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc7:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
Success	cm64-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
rs6000-aix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
x86_64-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.

mfi.ipmatr

i686-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
mips-irix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc7:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
Success	cm64-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
rs6000-aix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
x86_64-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.

phi.ipmatr

i686-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
mips-irix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc7:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
Success	cm64-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
rs6000-aix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
x86_64-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.

Html last generated: Sun Mar 6 05:51:39 2016

Input last modified: Mon Jan 17 17:06:30 2011

CMISS Help / Examples / i / i4 / i41