/*****************************************************************  

  Description: 

   This mumax3 script template checks if the skyrmion gets through
   a VCMA gate for a given current density and a given increase 
   (percentage) of the uniaxial anisotropy due to the VCMA effect.

   This template contains three place holders:
        JX     current density along the x direction 
        VCMA   precentage increase of the anisotropy due to the VCMA
        TCUT   cut-off time

   This file is written for the mumax3 workshop organised for the 
   spintronic seminar series in the summer of 2020. For more 
   information go to https://www.mumax.ugent.be/mumax3-workshop
  
  Authors:

   Dr. Jonathan Leliaert (jonathan.leliaert@ugent.be)
   Dr. Jeroen Mulkers    (jeroen.mulkers@ugent.be)

******************************************************************

setgridsize(126,64,1)
setcellsize(1e-9,1e-9,1e-9) 
setpbc(4,0,0)

Msat  = 580e3
Aex   = 15e-12
Dind  = 3.0e-3
Ku1   = 0.8e6
AnisU = vector(0,0,1)
alpha = 0.1

m = neelskyrmion(-1, 1).transl(-40e-9,0,0)
minimize()

defRegion(1,rect(50e-9,inf))
Ku1.setRegion(1, (1+$VCMA/100.0)*Ku1.GetRegion(0) )

Pol = 0.4
xi  = 0.2
j   = vector(-$JX,0,0)

runwhile( ext_bubblepos.Get().X() < 0 && t < $TCUT )
    
if ext_bubblepos.Get().X() > 0 {
    print("PASSED!")
}  else {
    print("STOPPED!")
}