From bf54fa9a0f5fc70c9d6afef60d0f50e5ca32909c Mon Sep 17 00:00:00 2001
From: Qianqian Fang <fangqq@gmail.com>
Date: Fri, 13 Sep 2024 11:53:14 -0400
Subject: [PATCH] [bug] fix the potential typo in Custo et al for CSF mua,
 fangq/mcx#232

---
 examples/colin27/prop_brain.dat   | 2 +-
 mmclab/example/demo_head_atlas.m  | 2 +-
 mmclab/example/demo_mmcl_b2_b2d.m | 4 ++--
 mmclab/example/demo_mmcl_b3.m     | 2 +-
 mmclab/example/demo_sfdi_2layer.m | 2 +-
 5 files changed, 6 insertions(+), 6 deletions(-)

diff --git a/examples/colin27/prop_brain.dat b/examples/colin27/prop_brain.dat
index 3fc54785..d9ae0d38 100644
--- a/examples/colin27/prop_brain.dat
+++ b/examples/colin27/prop_brain.dat
@@ -1,5 +1,5 @@
 1 4
 1 0.019 7.8182   0.89 1.37
-2 0.004 0.00909 0.89 1.37
+2 0.0004 0.00909 0.89 1.37
 3 0.02  9.0   0.89 1.37
 4 0.08  40.9   0.84 1.37
diff --git a/mmclab/example/demo_head_atlas.m b/mmclab/example/demo_head_atlas.m
index 482a7391..842e30f8 100644
--- a/mmclab/example/demo_head_atlas.m
+++ b/mmclab/example/demo_head_atlas.m
@@ -29,7 +29,7 @@
 cfg.node = double(node);
 cfg.elem = double(elem);
 cfg.elemprop = double(prop);
-cfg.prop = [0, 0, 1, 1; 0, 0, 1, 1; 0.019 7.8 0.89 1.37; 0.019 7.8 0.89 1.37; 0.004 0.009 0.89 1.37; 0.02 9.0 0.89 1.37; 0.08 40.9 0.84 1.37];
+cfg.prop = [0, 0, 1, 1; 0, 0, 1, 1; 0.019 7.8 0.89 1.37; 0.019 7.8 0.89 1.37; 0.0004 0.009 0.89 1.37; 0.02 9.0 0.89 1.37; 0.08 40.9 0.84 1.37];
 
 % light source
 cfg.srctype = 'pencil';
diff --git a/mmclab/example/demo_mmcl_b2_b2d.m b/mmclab/example/demo_mmcl_b2_b2d.m
index 3c9aa06c..ffbb8ed8 100644
--- a/mmclab/example/demo_mmcl_b2_b2d.m
+++ b/mmclab/example/demo_mmcl_b2_b2d.m
@@ -27,7 +27,7 @@
 cfg.tstart = 0;
 cfg.tend = 5e-9;
 cfg.tstep = 5e-10;
-cfg.prop = [0, 0, 1, 1; 0.02 7.0 0.89 1.37; 0.004 0.009 0.89 1.37; 0.02 9.0 0.89 1.37; 0.05 0.0 1.0 1.37];
+cfg.prop = [0, 0, 1, 1; 0.02 7.0 0.89 1.37; 0.0004 0.009 0.89 1.37; 0.02 9.0 0.89 1.37; 0.05 0.0 1.0 1.37];
 cfg.debuglevel = 'TP';
 cfg.isreflect = 1;
 
@@ -112,7 +112,7 @@
 cfg_mcx.issrcfrom0 = 1;
 
 % format: [mua(1/mm) mus(1/mm) g n]
-cfg_mcx.prop = [0, 0, 1, 1; 0.02 7.0 0.89 1.37; 0.004 0.009 0.89 1.37; 0.02 9.0 0.89 1.37; 0.05 0.0 1.0 1.37];
+cfg_mcx.prop = [0, 0, 1, 1; 0.02 7.0 0.89 1.37; 0.0004 0.009 0.89 1.37; 0.02 9.0 0.89 1.37; 0.05 0.0 1.0 1.37];
 
 % time-domain simulation parameters
 cfg_mcx.tstart = 0;
diff --git a/mmclab/example/demo_mmcl_b3.m b/mmclab/example/demo_mmcl_b3.m
index 81fcfc52..e3743dbd 100644
--- a/mmclab/example/demo_mmcl_b3.m
+++ b/mmclab/example/demo_mmcl_b3.m
@@ -35,7 +35,7 @@
 cfg.prop = [0         0    1.0000    1.0000 % background/air
             0.0190    7.8182    0.8900    1.3700 % scalp & skull are both labeled using 1
             %     0.0190    7.8182    0.8900    1.3700 % skull
-            0.0040    0.0090    0.8900    1.3700 % csf
+            0.0004    0.0090    0.8900    1.3700 % csf
             0.0200    9.0000    0.8900    1.3700 % gray matters
             0.0800   40.9000    0.8400    1.3700 % white matters
             0         0    1.0000    1.0000]; % air pockets
diff --git a/mmclab/example/demo_sfdi_2layer.m b/mmclab/example/demo_sfdi_2layer.m
index 80d4e0b7..f5ec60cd 100644
--- a/mmclab/example/demo_sfdi_2layer.m
+++ b/mmclab/example/demo_sfdi_2layer.m
@@ -29,7 +29,7 @@
 else
     [node, face, c0] = latticegrid([0 60], [0 60], [0 20 25 30]); % if you like a 3-layer model
     c0(:, 4) = [2; 2; 1];
-    cfg.prop = [0 0 1 1; 0.02 9.0, 0.89 1.37; 0.004 0.009, 0.89 1.37; 0.019 7.8 0.89 1.37];
+    cfg.prop = [0 0 1 1; 0.02 9.0, 0.89 1.37; 0.0004 0.009, 0.89 1.37; 0.019 7.8 0.89 1.37];
 end
 [cfg.node, cfg.elem] = surf2mesh(node, face, [], [], 1, [], c0);