logical res = Gravity_bhMAC( real(in) :: node(:), real(in) :: ndSize2, real(in) :: dr(MDIM+2), integer(in) :: blockno, integer(in) :: point(MDIM), integer(in) :: blkLimits(2,MDIM), real,pointer :: solnData(:,:,:,:) )
Multipole Acceptance Criterion. Determines whether the contribution of the node to the potential at the point of calculation will have a nacessary accuracy.
node : array of the node tested ndSize2 : square of the physical size of the node (the largest extent of the node) dr : (1:MDIM) - position vector from the point-of-calculation to the node (MDIM+1) - square of the magnitude of the position vector (MDIM+2) - inverted magnitude of the position vector blockno : number of block into which the point-of-calculation belongs point : indeces of the point-of-calculation in the block blkLimits : limits of indeces in the block solnData : solution data from the grid
Returns TRUE if the node is accepted for the calculation of the potential. Otherwise returns FALSE.