matRad IPOPT callback: objective function for direct aperture optimization call f = matRad_objectiveFunction(optiProb,apertureInfoVect,dij,cst) input optiProb: option struct defining the type of optimization apertureInfoVect: aperture info in form of vector dij: matRad dij struct as generated by bixel-based dose calculation cst: matRad cst struct output f: objective function value References [1] http://dx.doi.org/10.1118/1.4914863 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0001 function f = matRad_objectiveFunction(optiProb,apertureInfoVect,dij,cst) 0002 % matRad IPOPT callback: objective function for direct aperture optimization 0003 % 0004 % call 0005 % f = matRad_objectiveFunction(optiProb,apertureInfoVect,dij,cst) 0006 % 0007 % input 0008 % optiProb: option struct defining the type of optimization 0009 % apertureInfoVect: aperture info in form of vector 0010 % dij: matRad dij struct as generated by bixel-based dose calculation 0011 % cst: matRad cst struct 0012 % 0013 % output 0014 % f: objective function value 0015 % 0016 % References 0017 % [1] http://dx.doi.org/10.1118/1.4914863 0018 % 0019 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0020 0021 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0022 % 0023 % Copyright 2015 the matRad development team. 0024 % 0025 % This file is part of the matRad project. It is subject to the license 0026 % terms in the LICENSE file found in the top-level directory of this 0027 % distribution and at https://github.com/e0404/matRad/LICENSES.txt. No part 0028 % of the matRad project, including this file, may be copied, modified, 0029 % propagated, or distributed except according to the terms contained in the 0030 % LICENSE file. 0031 % 0032 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 0033 0034 apertureInfo = optiProb.apertureInfo; 0035 0036 % update apertureInfo, bixel weight vector an mapping of leafes to bixels 0037 apertureInfo = optiProb.matRad_daoVec2ApertureInfo(apertureInfo,apertureInfoVect); 0038 0039 % bixel based objective function calculation 0040 f = matRad_objectiveFunction@matRad_OptimizationProblem(optiProb,apertureInfo.bixelWeights,dij,cst);