timedep excitation extra test file and fix

engine/excitation.go

@@ -4,6 +4,7 @@ import (
 	"github.com/mumax/3/cuda"
 	"github.com/mumax/3/data"
 	"github.com/mumax/3/script"
+	"log"
 	"reflect"
 )
 
@@ -62,12 +63,13 @@ func assureGPU(s *data.Slice) *data.Slice {
 }
 
 // user script: has to be 3-vector
-func (e *excitation) SetRegion(region int, value [3]float64) {
-	e.perRegion.setRegion(region, value[:])
+func (e *excitation) SetRegion(region int, f script.VectorFunction) {
+	e.perRegion.SetRegion(region, f)
 }
 
 // for gui (nComp agnostic)
 func (e *excitation) setRegion(region int, value []float64) {
+	log.Println("TODO: time-dep gui")
 	e.perRegion.setRegion(region, value)
 }
 

engine/vectorparam.go

@@ -3,6 +3,7 @@ package engine
 import (
 	"github.com/mumax/3/script"
 	"github.com/mumax/3/util"
+	"log"
 	"reflect"
 )
 
@@ -26,8 +27,10 @@ func (p *VectorParam) SetValue(v interface{}) {
 
 func (p *VectorParam) setRegionsFunc(r1, r2 int, f script.VectorFunction) {
 	if f.Const() {
+		log.Println(p, "[", r1, ":", r2, "]", "is constant")
 		p.setRegions(r1, r2, slice(f.Float3()))
 	} else {
+		log.Println(p, "[", r1, ":", r2, "]", "is not constant")
 		p.setFunc(r1, r2, func() []float64 {
 			return slice(f.Float3())
 		})

test/timedep3Region.txt

@@ -0,0 +1,26 @@
+c:=4e-9
+setgridsize(32, 32, 1)
+setcellsize(c, c, c)
+
+Msat=860e3
+Aex=13e-12
+alpha=0.2
+m=uniform(1, 1, 0)
+setgeom(circle(32*4e-9))
+
+f := 1e9
+A := 0.01
+defRegion(1, xrange(-inf, inf))
+B_ext.setRegion(1, vector(A*sin(2*pi*f*t), A*cos(2*pi*f*t), 0))
+
+tableadd(B_ext)
+tableautosave(1e-12)
+
+run(0.2e-9)
+
+m_ := average(m)
+print(m_)
+
+expect("mx", m_[0], 0.6783973430765087, 1e-7)
+expect("my", m_[1], 0.7192199406365456, 1e-7)
+expect("mz", m_[2], 0.0212674352335812, 1e-7)