dc_scaledsec::mod Interface Reference

Public Member Functions
type(dc_scaled_sec) function	dcscaledsec_mod_si (sclsec, factor)
type(dc_scaled_sec) function	dcscaledsec_mod_sr (sclsec, factor)
type(dc_scaled_sec) function	dcscaledsec_mod_sd (sclsec, factor)
type(dc_scaled_sec) function	dcscaledsec_mod_ss (sclsec, factor)

Detailed Description

Definition at line 221 of file dc_scaledsec.f90.

Member Function/Subroutine Documentation

dcscaledsec_mod_sd()

type(dc_scaled_sec) function dc_scaledsec::mod::dcscaledsec_mod_sd	(	type(dc_scaled_sec), intent(in)	sclsec,
		real(dp), intent(in)	factor
	)

Calculate remainder of DC_SCALED_SEC divided by double

Parameters

[in]	sclsec	DC_SCALED_SEC variable (dividend)
[in]	factor	Double precision divisor

Returns

Remainder of division

Definition at line 2190 of file dc_scaledsec.f90.

    use dc_message, only: messagenotify
    implicit none
    type(DC_SCALED_SEC), intent(in):: sclsec
    real(DP), intent(in):: factor
    type(DC_SCALED_SEC):: factor_scl
 
  continue
    factor_scl = factor
    result = mod( sclsec, factor_scl )

dcscaledsec_mod_si()

type(dc_scaled_sec) function dc_scaledsec::mod::dcscaledsec_mod_si	(	type(dc_scaled_sec), intent(in)	sclsec,
		integer, intent(in)	factor
	)

Calculate remainder of DC_SCALED_SEC divided by integer

Parameters

[in]	sclsec	DC_SCALED_SEC variable (dividend)
[in]	factor	Integer divisor

Returns

Remainder of division

Definition at line 2136 of file dc_scaledsec.f90.

    use dc_message, only: messagenotify
    implicit none
    type(DC_SCALED_SEC), intent(in):: sclsec
    integer, intent(in):: factor
    type(DC_SCALED_SEC):: factor_scl
 
  continue
    factor_scl = factor
    result = mod( sclsec, factor_scl )

dcscaledsec_mod_sr()

type(dc_scaled_sec) function dc_scaledsec::mod::dcscaledsec_mod_sr	(	type(dc_scaled_sec), intent(in)	sclsec,
		real, intent(in)	factor
	)

Calculate remainder of DC_SCALED_SEC divided by real

Parameters

[in]	sclsec	DC_SCALED_SEC variable (dividend)
[in]	factor	Single precision divisor

Returns

Remainder of division

Definition at line 2163 of file dc_scaledsec.f90.

    use dc_message, only: messagenotify
    implicit none
    type(DC_SCALED_SEC), intent(in):: sclsec
    real, intent(in):: factor
    type(DC_SCALED_SEC):: factor_scl
 
  continue
    factor_scl = factor
    result = mod( sclsec, factor_scl )

dcscaledsec_mod_ss()

type(dc_scaled_sec) function dc_scaledsec::mod::dcscaledsec_mod_ss	(	type(dc_scaled_sec), intent(in)	sclsec,
		type(dc_scaled_sec), intent(in)	factor
	)

Calculate remainder of DC_SCALED_SEC division

Note: Due to precision limitations in Fortran compilers, factor must be smaller than 10^12.

Parameters

[in]	sclsec	DC_SCALED_SEC variable (dividend)
[in]	factor	DC_SCALED_SEC variable (divisor)

Returns

Remainder of division

Definition at line 2032 of file dc_scaledsec.f90.

    use dc_message, only: messagenotify
    use dc_types, only: dp_eps
    implicit none
    type(DC_SCALED_SEC), intent(in):: sclsec, factor
 
    type(DC_SCALED_SEC):: factor_scl
    real(DP):: sec_ary_mod(imin+imin:imax)
    integer:: i, move_down_index, sf_idx
    real(DP):: move_down
    real(DP):: factor_dp
    type(DC_SCALED_SEC):: zero_sec
    logical:: done
  continue
 
    ! factor must be smaller than 10^12 due to compiler precision limitations.
    !
    if ( .not. all( factor % sec_ary (imax-4:imax) == (/ 0, 0, 0, 0, 0 /) ) ) then
      call messagenotify( 'E', 'dc_scaledsec#mod', &
        & 'factor must be smaller than 10^12' )
    end if
 
    if ( sclsec == factor ) then
      result = zero_sec
      return
    end if
 
    factor_scl % sec_ary(imin:-1) = 0
    factor_scl % sec_ary(imin-imin:imax) = factor % sec_ary(imin:imax+imin)
    factor_scl % flag_negative = factor % flag_negative
 
    factor_dp = factor_scl
 
    done = .false.
    move_down = 0.0_dp
    do i = imax, imin, -1
      move_down_index = i
      if ( abs(move_down) > dp_eps ) then
        sf_idx = i - imin
        if ( sf_idx > ubound(scale_factor_xx, 1) ) sf_idx = ubound(scale_factor_xx, 1)
        if ( sf_idx < lbound(scale_factor_xx, 1) ) sf_idx = lbound(scale_factor_xx, 1)
        if ( abs(factor_dp) > ( move_down + scale_factor ) * scale_factor_xx( sf_idx ) ) then
          done = .true.
          exit
        end if
      end if
 
      sec_ary_mod(i) = &
        & mod( ( sclsec % sec_ary(i) + move_down ), factor_dp )
 
      if ( abs(sec_ary_mod(i)) > dp_eps ) then
        move_down = sec_ary_mod(i) * scale_factor
      else
        move_down = 0.0_dp
      end if
    end do
 
    if ( .not. done ) then
      do i = imin - 1, imin + imin, -1
      move_down_index = i
      if ( abs(move_down) > dp_eps ) then
        sf_idx = i - imin
        if ( sf_idx > ubound(scale_factor_xx, 1) ) sf_idx = ubound(scale_factor_xx, 1)
        if ( sf_idx < lbound(scale_factor_xx, 1) ) sf_idx = lbound(scale_factor_xx, 1)
        if ( abs(factor_dp) > ( move_down + scale_factor ) * scale_factor_xx( sf_idx ) ) then
          done = .true.
          exit
        end if
      end if
 
      sec_ary_mod(i) = mod( move_down, factor_dp )
 
      if ( abs(sec_ary_mod(i)) > dp_eps ) then
        move_down = sec_ary_mod(i) * scale_factor
      else
        move_down = 0.0_dp
      end if
      end do
    end if
 
    result = move_down * scale_factor_xx(move_down_index)
    if ( move_down_index > imin - 1 ) then
      result % sec_ary(imin:move_down_index) = sclsec % sec_ary(imin:move_down_index)
    end if
 
    result % flag_negative = sclsec % flag_negative
 

References dc_types::dp_eps.

---

The documentation for this interface was generated from the following file:

• dc_utils/dc_scaledsec.f90