Field access checks

models/electrokinetic/d2q9_poison_boltzmann/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 
 AddDensity( name="g[0]", dx= 0, dy= 0, group="g")
 AddDensity( name="g[1]", dx= 1, dy= 0, group="g")

models/experimental/d3q27_csf/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 
 U = expand.grid(-1:1,-1:1,-1:1)
 

models/flow/auto/Dynamics.R

@@ -57,8 +57,9 @@ AddNodeType(name="Body", group="BODY")
 for (f in fname) AddField(f,dx=0,dy=0,dz=0) # Make f accessible also in present node (not only streamed)
 
 if (Options$part) {
-	AddStage("BaseIteration", "Run", save=Fields$group %in% c("f","Force"), load = DensityAll$group %in% c("f","Force"))
-	AddStage("CalcF", save=Fields$group == "Force", load = DensityAll$group %in% c("f","Force"), particle=TRUE)
+	AddStage("BaseIteration", "Run", save=Fields$group %in% c("f"), load = DensityAll$group %in% c("f","Force"))
+	AddStage("BaseInit", "Init", save=Fields$group %in% c("f"))
+	AddStage("CalcF", save=Fields$group == "Force", load = DensityAll$group %in% c("f"), particle=TRUE)
 	AddAction("Iteration", c("BaseIteration", "CalcF"))
 	AddAction("Init", c("BaseInit", "CalcF"))
 }

models/flow/d2q9_par/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Density - table of variables of LB Node to stream
 #  name - variable name to stream
 #  dx,dy,dz - direction of streaming

models/flow/d3q27_cumulant_part/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 x = c(0,1,-1);
 P = expand.grid(x=0:2,y=0:2,z=0:2)
 U = expand.grid(x,x,x)

models/heat/d3q27_tePSM_per/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 
 AddDensity( name="f[0]", dx= 0, dy= 0, dz= 0, group="f")
 AddDensity( name="f[1]", dx= 1, dy= 0, dz= 0, group="f")

models/multiphase/d2q9_ShanChen/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 
 # Fluid Density Populations
 AddDensity( name="f[0]", dx= 0, dy= 0, group="f")

models/multiphase/d2q9_csf/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Density - table of variables of LB Node to stream
 #  name - variable name to stream
 #  dx,dy,dz - direction of streaming

models/multiphase/d2q9_kuper/Dynamics.R

@@ -13,7 +13,7 @@ AddField("phi",stencil2d=1);
 
 AddStage("BaseIteration", "Run", save=Fields$group == "f", load=DensityAll$group == "f")
 AddStage("CalcPhi", save="phi",load=DensityAll$group == "f")
-AddStage("BaseInit", "Init", save=Fields$group == "f", load=DensityAll$group == "f")
+AddStage("BaseInit", "Init", save=Fields$group == "f")
 
 AddAction("Iteration", c("BaseIteration","CalcPhi"))
 AddAction("Init", c("BaseInit","CalcPhi"))

models/multiphase/d2q9_pf/Dynamics.R

@@ -43,17 +43,16 @@ if (Options$fd) {
     AddField("phi"      ,stencil2d=1 );
 
     AddStage("BaseIteration", "Run", 
-             load=DensityAll$group == "h" | DensityAll$group == "f" | DensityAll$group == "BC",
-             save=Fields$group=="h" | Fields$group=="f"  
-         ) 
+        load=DensityAll$group == "h" | DensityAll$group == "f" | DensityAll$group == "BC",
+        save=Fields$group=="h" | Fields$group=="f" | Fields$group == "BC"
+    ) 
     AddStage("CalcPhi", 
-             save=Fields$name=="phi" ,  
-             load=DensityAll$group == "h"
-         )
+        save=Fields$name=="phi" ,  
+        load=DensityAll$group == "h"
+    )
     AddStage("BaseInit", "Init",  
-	    load=DensityAll$group == "BC",
-        save=Fields$group=="h" | Fields$group == "f" 
-
+	    load=FALSE,
+        save=Fields$group=="h" | Fields$group == "f" | Fields$group == "BC"
     )
     AddAction("Iteration", c("BaseIteration","CalcPhi"))
     AddAction("Init", c("BaseInit","CalcPhi"))

models/multiphase/d2q9_pf_pressureEvolution/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Density - table of variables of LB Node to stream
 #  name - variable name to stream
 #  dx,dy,dz - direction of streaming

models/multiphase/d2q9_pf_velocity/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # 	Density - table of variables of LB Node to stream
 #	Pressure Evolution:
 AddDensity( name="g[0]", dx= 0, dy= 0, group="g")

models/multiphase/d2q9_scmp/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 
 AddDensity( name="f[0]", dx= 0, dy= 0, group="f")
 AddDensity( name="f[1]", dx= 1, dy= 0, group="f")

models/multiphase/d3q19_kuper/Dynamics.R

@@ -38,7 +38,7 @@ AddField("phi",stencil3d=1);
 
 AddStage("BaseIteration", "Run", save=Fields$group == "f", load=DensityAll$group == "f")
 AddStage("CalcPhi", save="phi",load=DensityAll$group == "f")
-AddStage("BaseInit", "Init", save=Fields$group == "f", load=DensityAll$group == "f")
+AddStage("BaseInit", "Init", save=Fields$group == "f")
 
 AddAction("Iteration", c("BaseIteration","CalcPhi"))
 AddAction("Init", c("BaseInit","CalcPhi"))

models/multiphase/d3q27_pf_velocity/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Density - table of variables of LB Node to stream
 #	Velocity-based Evolution d3q27:
 source("lattice.R")

models/multiphase/experimental/d2q9_pp_LBL/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Fluid Density Populations
 AddDensity( name="f[0]", dx= 0, dy= 0, group="f")
 AddDensity( name="f[1]", dx= 1, dy= 0, group="f")

models/multiphase/experimental/d2q9_pp_MCMP/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 # Density - table of variables of LB Node to stream
 AddDensity( name="f[0]", dx= 0, dy= 0, group="f")
 AddDensity( name="f[1]", dx= 1, dy= 0, group="f")

models/reaction/d2q9_reaction_diffusion_system/Dynamics.R

@@ -1,3 +1,10 @@
+# Setting permissive access policy.
+#  * This skips checks of fields being overwritten or read prematurely.
+#  * Otherwise the model compilation was failing.
+#  * This should be removed if the issue is fixed
+SetOptions(permissive.access=TRUE)  ### WARNING
+
+
 if (Options$AllenCahn) {
  	Qname = 'Allen-Cahn'
 	DREs <- ('PHI')

src/LatticeAccess.inc.cpp.Rt

@@ -204,8 +204,8 @@ CudaDeviceFunction flag_t LatticeContainer::getType(int x, int y, int z) const
 	)
  ))
  all_stages = lapply(all_stages, function(s) {
-  if (is.null(s$load.densities))  s$load.densities = DensityAll[,s$tag];
-  if (is.null(s$save.fields))  s$save.fields = Fields[,s$tag];
+  if (is.null(s$load.densities))  s$load.densities = DensityAll[,s$loadtag];
+  if (is.null(s$save.fields))  s$save.fields = Fields[,s$savetag];
   if (is.null(s$suffix))  s$suffix = paste("_",s$name,sep="")
   s
  })
@@ -225,6 +225,7 @@ CudaDeviceFunction flag_t LatticeContainer::getType(int x, int y, int z) const
   real_to_storage = function(val,f) storage_convert("real_to_storage",val,f)
 
 
+
 resolve.symmetries = function(D) {
   if (nrow(D) < 1) return();
   dp = D[,c("dx","dy","dz"),drop=FALSE]

src/SUMMARY.Rt

@@ -1,4 +1,5 @@
 <?R
+  plot.access=TRUE
   source("conf.R")
   save.image(file="SUMMARY.Rdata")
   options(width=4000)
@@ -45,4 +46,6 @@
     cat(sprintf("%-35s [%7s] | %1s | %s",s$name, s$unit, ifelse(s$adjoint, "X"," "), s$comment),"\n");
   }
   cat("---------------------------------------------------------------------------------------\n");
+  cat("\n");
+  cat("Field access plot: ", paste0("file://",getwd(),"/field_access.pdf"),"\n")
 ?>