From 42327643e7c14fef10943dabdb1776a6807d1d0b Mon Sep 17 00:00:00 2001
From: Qianqian Fang <fangqq@gmail.com>
Date: Mon, 1 Jan 2024 20:52:03 -0500
Subject: [PATCH] [mcxlab] add speed comparison between different media formats

---
 .../examples/demo_speed_label_vs_continuous.m | 53 +++++++++++++++++++
 1 file changed, 53 insertions(+)
 create mode 100644 mcxlab/examples/demo_speed_label_vs_continuous.m

diff --git a/mcxlab/examples/demo_speed_label_vs_continuous.m b/mcxlab/examples/demo_speed_label_vs_continuous.m
new file mode 100644
index 00000000..12185e43
--- /dev/null
+++ b/mcxlab/examples/demo_speed_label_vs_continuous.m
@@ -0,0 +1,53 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% MCXLAB - Monte Carlo eXtreme for MATLAB/Octave by Qianqina Fang
+%
+% In this example, we show why using label-based simulation is preferred
+% over per-voxel medium formats whenever possible.
+%
+% This file is part of Monte Carlo eXtreme (MCX) URL:http://mcx.sf.net
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+% only clear cfg to avoid accidentally clearing other useful data
+clear cfg cfgs;
+cfg.nphoton = 1e7;
+cfg.vol = uint8(ones(60, 60, 60));
+cfg.vol(:, :, 4:10) = 2;
+cfg.srcpos = [30 30 1];
+cfg.srcdir = [0 0 1];
+cfg.detpos = [30 33 1 2];
+cfg.gpuid = 1;
+% cfg.gpuid='11'; % use two GPUs together
+cfg.autopilot = 1;
+cfg.prop = [0 0 1 1; 0.05 1 0 1.37; 0.1 1 0 1.37];
+cfg.tstart = 0;
+cfg.tend = 5e-9;
+cfg.tstep = 5e-9;
+
+%% simulate with lable-based medium
+tic;
+flux = mcxlab(cfg);
+toc;
+
+%% simulate with mua_float format
+cfg2 = cfg;
+cfg2.vol = reshape(single(cfg.prop(cfg.vol + 1, 1)), [1 size(cfg.vol)]);
+
+tic;
+flux2 = mcxlab(cfg2);
+toc;
+
+%% simulate with muamus_float format
+cfg3 = cfg;
+cfg3.vol = reshape(single(cfg.prop(cfg.vol + 1, 1:2))', [2 size(cfg.vol)]);
+
+tic;
+flux3 = mcxlab(cfg3);
+toc;
+
+%% simulate with asgn_float format
+cfg4 = cfg;
+cfg4.vol = reshape(single(cfg.prop(cfg.vol + 1, 1:4))', [4 size(cfg.vol)]);
+
+tic;
+flux4 = mcxlab(cfg4);
+toc;