dc_scaledsec::modulo Interface Reference

Private Member Functions
type(dc_scaled_sec) function	dcscaledsec_modulo_si (sclsec, factor)
type(dc_scaled_sec) function	dcscaledsec_modulo_sr (sclsec, factor)
type(dc_scaled_sec) function	dcscaledsec_modulo_sd (sclsec, factor)
type(dc_scaled_sec) function	dcscaledsec_modulo_ss (sclsec, factor)

Detailed Description

Definition at line 165 of file dc_scaledsec.f90.

Member Function/Subroutine Documentation

dcscaledsec_modulo_sd()

type(dc_scaled_sec) function dc_scaledsec::modulo::dcscaledsec_modulo_sd ( type(dc_scaled_sec), intent(in)  sclsec, real(dp), intent(in)  factor  )

private

Definition at line 1872 of file dc_scaledsec.f90.

    !
    ! DC_SCALED_SEC 型変数を割った際の剰余を計算. 
    ! 
    ! Calculate of remainder by division of a "DC_SCALED_SEC" variable
    !
    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 = modulo( sclsec, factor_scl )

dcscaledsec_modulo_si()

type(dc_scaled_sec) function dc_scaledsec::modulo::dcscaledsec_modulo_si ( type(dc_scaled_sec), intent(in)  sclsec, integer, intent(in)  factor  )

private

Definition at line 1834 of file dc_scaledsec.f90.

    !
    ! DC_SCALED_SEC 型変数を割った際の剰余を計算. 
    ! 
    ! Calculate of remainder by division of a "DC_SCALED_SEC" variable
    !
    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 = modulo( sclsec, factor_scl )

dcscaledsec_modulo_sr()

type(dc_scaled_sec) function dc_scaledsec::modulo::dcscaledsec_modulo_sr ( type(dc_scaled_sec), intent(in)  sclsec, real, intent(in)  factor  )

private

Definition at line 1853 of file dc_scaledsec.f90.

    !
    ! DC_SCALED_SEC 型変数を割った際の剰余を計算. 
    ! 
    ! Calculate of remainder by division of a "DC_SCALED_SEC" variable
    !
    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 = modulo( sclsec, factor_scl )

dcscaledsec_modulo_ss()

type(dc_scaled_sec) function dc_scaledsec::modulo::dcscaledsec_modulo_ss ( type(dc_scaled_sec), intent(in)  sclsec, type(dc_scaled_sec), intent(in)  factor  )

private

Definition at line 1747 of file dc_scaledsec.f90.

    !
    ! DC_SCALED_SEC 型変数を割った際の剰余を計算. 
    ! 
    ! Calculate of remainder by division of a "DC_SCALED_SEC" variable
    !
    use dc_message, only: messagenotify
    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
    real(DP):: move_down
    real(DP):: factor_dp
    type(dc_scaled_sec):: zero_sec
  continue
  
    ! frt, ifort などでは, 1.0e+23 などの実数とすると, 
    ! 9.9999e+22 などとなってしまうため, 
    ! factor として指定するものは 10e+12 までとする. (うーむ, 汚い対応だな....)
    ! (morikawa 2008/09/01) 
    !
    if ( .not. all( factor % sec_ary (imax-4:imax) == (/ 0, 0, 0, 0, 0 /) ) ) then
      call messagenotify( 'E', 'dc_scaledsec#modulo', &
        & '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

    move_down = 0.0_dp
    do i = imax, imin + imin, -1
      move_down_index = i
      if ( move_down /= 0.0_dp ) then
        if ( abs(factor_dp) > ( move_down + scale_factor ) * scale_factor_xx( i - imin ) ) exit
      end if

      if ( i > imin - 1 ) then
        sec_ary_mod(i) = &
          & mod( ( sclsec % sec_ary(i) + move_down ), factor_dp )
      else
        sec_ary_mod(i) = mod( move_down, factor_dp )
      end if

      if ( sec_ary_mod(i) /= 0.0_dp ) then
        move_down = sec_ary_mod(i) * scale_factor
      else
        move_down = 0.0_dp
      end if

    end do

    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 = .false.

    if ( .not. result == zero_sec ) then
      if ( .not. sclsec % flag_negative .and. factor % flag_negative ) then
        result = - factor - result
        result % flag_negative = .not. sclsec % flag_negative
        
      elseif ( sclsec % flag_negative .and. .not. factor % flag_negative ) then
        result = factor - result
        result % flag_negative = .not. sclsec % flag_negative

      else
        result % flag_negative = sclsec % flag_negative
        
      end if
    end if

---

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

• dc_utils/dc_scaledsec.f90