Update matRad_fluenceOptimization.m

Some fields necessary for RBExD optimization have to be written if initial weights are passed in
HomolkaN · Dec 1, 2022 · 74f4f5d · 74f4f5d
1 parent 093eb96
commit 74f4f5d
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/matRad_fluenceOptimization.m b/matRad_fluenceOptimization.m
@@ -103,6 +103,15 @@
     %do nothing as wInit was passed to the function
     matRad_cfg.dispInfo('chosen provided wInit!\n');
 
+    % Write ixDose which is needed for the optimizer
+    if pln.bioParam.bioOpt
+        dij.ixDose  = dij.bx~=0;
+
+        %pre-calculations
+        dij.gamma             = zeros(dij.doseGrid.numOfVoxels,dij.numOfScenarios);
+        dij.gamma(dij.ixDose) = dij.ax(dij.ixDose)./(2*dij.bx(dij.ixDose));
+    end
+
 elseif strcmp(pln.bioParam.model,'constRBE') && strcmp(pln.radiationMode,'protons')
     % check if a constant RBE is defined - if not use 1.1
     if ~isfield(dij,'RBE')