This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Last revision Both sides next revision | ||
howtos:feedback:tracking_method_-_minimizing_root_mean_squared_difference_with_5_point_bracketing [2010/12/24 16:07] shona.weldon |
howtos:feedback:tracking_method_-_minimizing_root_mean_squared_difference_with_5_point_bracketing [2011/11/04 18:46] marcus.williams |
||
---|---|---|---|
Line 1: | Line 1: | ||
- | ====== Minimizing Root Mean Squared Difference with 5 point bracketing ====== | + | ====== Minimizing Root Mean Squared Difference with 5 point bracketing (bi-section variant) ====== |
A typical calibration problem is to get the model results to track on something while tweaking a variable. The approach described here attempts to minimize the difference over time between the model variable and some given tracking variable. | A typical calibration problem is to get the model results to track on something while tweaking a variable. The approach described here attempts to minimize the difference over time between the model variable and some given tracking variable. | ||
Line 142: | Line 142: | ||
integer $numEvalsCompleted | integer $numEvalsCompleted | ||
getobjinfo (numEvalsCompleted[], $numEvalsCompleted; info=dataElement) | getobjinfo (numEvalsCompleted[], $numEvalsCompleted; info=dataElement) | ||
- | say (" number of evaluations completed ",$fbViewIteration ) | + | say (" number of evaluations completed ",$numEvalsCompleted) |
if ($numEvalsCompleted > 1) | if ($numEvalsCompleted > 1) | ||
! Can look at last evaluations results compared to this one then store this one for future comparisons | ! Can look at last evaluations results compared to this one then store this one for future comparisons | ||
Line 150: | Line 150: | ||
!check for stoping | !check for stoping | ||
- | local numberNotConverged[] = sum (boolgt (abs (2 * currentMeshStep[upperDim,auxDims] / currentBestGuess[upperDim,auxDims]), tolerance[])) | + | local numberNotConverged[] = sum (boolgt (abs (2 * currentMeshStep[upperDim,auxDims] / currentBestGuess[upperDim,auxDims]), \ |
+ | tolerance[]); dim=all) | ||
integer $numberNotConverged | integer $numberNotConverged | ||
getobjinfo (numberNotConverged[], $numberNotConverged; info=dataElement, default=1) | getobjinfo (numberNotConverged[], $numberNotConverged; info=dataElement, default=1) |