|
@@ -77,8 +77,10 @@ switch lower(mode)
|
|
|
if k == 1 % allocate dose matrix in the first loop
|
|
|
dose = zeros(beamlet.x_count, beamlet.y_count, beamlet.z_count);
|
|
|
end
|
|
|
+% dose(beamlet.non_zero_indices) = dose(beamlet.non_zero_indices) + ...
|
|
|
+% beamlet.non_zero_values * weights(beamlet.num);
|
|
|
dose(beamlet.non_zero_indices) = dose(beamlet.non_zero_indices) + ...
|
|
|
- beamlet.non_zero_values * weights(beamlet.num);
|
|
|
+ beamlet.non_zero_values * weights(k);
|
|
|
end
|
|
|
|
|
|
results = beamlets;
|
|
@@ -123,7 +125,7 @@ switch lower(mode)
|
|
|
xpg = xprogressbar();
|
|
|
for k=1:nBmax%Nbeamlets
|
|
|
xpg.update(k/Nbeamlets);
|
|
|
- beamlet.num = fread(fid,1,'int');
|
|
|
+ beamlet.num = fread(fid,1,'int')
|
|
|
beamlet.x_count = fread(fid,1,'int');
|
|
|
beamlet.y_count = fread(fid,1,'int');
|
|
|
beamlet.z_count = fread(fid,1,'int');
|