RepeatCurrentTimeStep_AdditionalComputationsAfterDataTransfer
additional computations on original pointcloud after data transfer is finished
As the repeating of the current time step happens after the original time step is already finished,
we need a means of computing additional values on the original pointcloud. WE MUST NOT use
CODI_eq (see
CODI), as this function is processed
DURING the executions of the original time step.
RepeatCurrentTimeStep(n) = (
%RepeatCurrentTimeStep_BasedOnSamePointCloud% )
# repeatine based on the same pointcloud
RepeatCurrentTimeStep_ChangeCVconfiguration(n) = ( "ord_laplace=2" ,0, "ord_gradient=2" )
# set the approximation order (temporarily to 2 for the repeated time step)
RepeatCurrentTimeStep_SaveVariables(n) = ( %indU_v(1)%,
%ind_v(1)%,
# define data trensfer from temporary to original point cloud
%indU_v(2)%,
%ind_v(2)%,
%indU_v(3)%,
%ind_v(3)%,
%indU_p_corr%,
%ind_p_dyn%,
%indU_c%,
%ind_c% )
# perform additional computations on original pointcloud, based on the data transferred
RepeatCurrentTimeStep_AdditionalComputationsAfterDataTransfer(n) = ( %indU_Dv%, [ sqrt( (Y%indU_v(1)%-Y
%ind_v(1)%)^2 + (Y%indU_v(2)%-Y
%ind_v(2)%)^2 + (Y%indU_v(3)%-Y
%ind_v(3)%)^2 )] ,
%indU_DpCorr%, [ abs(Y%indU_p_corr%-Y
%ind_p_dyn%) ] ,
%indU_Dc%, [ abs(Y%indU_c%-Y
%ind_c%) ] )
In the example above,we compute the difference between the velocity, dynamic pressure, and correction pressure solutions between the
original time step and the additionally performed time step. The results of the computations are written to the index
variables %indU_Dv%, %indU_DpCorr%, and %indU_Dc%, respectively.
ATTENTION!!!!!!!
The additional computations are executed regardless of the order as they appear in the brackets, i.e. dependent solution cannot be produced. The following example
is constructed wrongly, as we presume a dependence of %indU_C% on %indU_A% and %indU_B% which
cannot be provided by RepeatCurrentTimeStep_AdditionalComputationsAfterDataTransfer .