dc_scaledsec::operator(*) Interface Reference

Public Member Functions
type(dc_scaled_sec) function	dcscaledsec_mul_ss (sclsec1, sclsec2)
type(dc_scaled_sec) function	dcscaledsec_mul_si (sclsec, factor)
type(dc_scaled_sec) function	dcscaledsec_mul_is (factor, sclsec)
type(dc_scaled_sec) function	dcscaledsec_mul_sd (sclsec, factor)
type(dc_scaled_sec) function	dcscaledsec_mul_ds (factor, sclsec)
type(dc_scaled_sec) function	dcscaledsec_mul_sr (sclsec, factor)
type(dc_scaled_sec) function	dcscaledsec_mul_rs (factor, sclsec)

Detailed Description

Definition at line 204 of file dc_scaledsec.f90.

Member Function/Subroutine Documentation

dcscaledsec_mul_ds()

type(dc_scaled_sec) function dc_scaledsec::operator(*)::dcscaledsec_mul_ds	(	real(dp), intent(in)	factor,
		type(dc_scaled_sec), intent(in)	sclsec )

Multiply double precision by DC_SCALED_SEC

Parameters

[in]	factor	Double precision multiplier
[in]	sclsec	DC_SCALED_SEC variable

Returns

Product of double and DC_SCALED_SEC

Definition at line 1789 of file dc_scaledsec.f90.

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

dcscaledsec_mul_is()

type(dc_scaled_sec) function dc_scaledsec::operator(*)::dcscaledsec_mul_is	(	integer, intent(in)	factor,
		type(dc_scaled_sec), intent(in)	sclsec )

Multiply integer by DC_SCALED_SEC

Parameters

[in]	factor	Integer multiplier
[in]	sclsec	DC_SCALED_SEC variable

Returns

Product of integer and DC_SCALED_SEC

Definition at line 1740 of file dc_scaledsec.f90.

    implicit none
    integer, intent(in):: factor
    type(DC_SCALED_SEC), intent(in):: sclsec
  continue
    result = sclsec * factor

dcscaledsec_mul_rs()

type(dc_scaled_sec) function dc_scaledsec::operator(*)::dcscaledsec_mul_rs	(	real, intent(in)	factor,
		type(dc_scaled_sec), intent(in)	sclsec )

Multiply single precision real by DC_SCALED_SEC

Parameters

[in]	factor	Single precision multiplier
[in]	sclsec	DC_SCALED_SEC variable

Returns

Product of real and DC_SCALED_SEC

Definition at line 1839 of file dc_scaledsec.f90.

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

dcscaledsec_mul_sd()

type(dc_scaled_sec) function dc_scaledsec::operator(*)::dcscaledsec_mul_sd	(	type(dc_scaled_sec), intent(in)	sclsec,
		real(dp), intent(in)	factor )

Multiply DC_SCALED_SEC by double precision

Parameters

[in]	sclsec	DC_SCALED_SEC variable
[in]	factor	Double precision multiplier

Returns

Product of DC_SCALED_SEC and double

Definition at line 1763 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 = sclsec * factor_scl

dcscaledsec_mul_si()

type(dc_scaled_sec) function dc_scaledsec::operator(*)::dcscaledsec_mul_si	(	type(dc_scaled_sec), intent(in)	sclsec,
		integer, intent(in)	factor )

Multiply DC_SCALED_SEC by integer

Optimized implementation without using mul_ss for performance.

Parameters

[in]	sclsec	DC_SCALED_SEC variable
[in]	factor	Integer multiplier

Returns

Product of DC_SCALED_SEC and integer

Definition at line 1681 of file dc_scaledsec.f90.

    use dc_message, only: messagenotify
    implicit none
    type(DC_SCALED_SEC), intent(in):: sclsec
    integer, intent(in):: factor
    integer:: factor_abs
    type(DC_SCALED_SEC):: zero_sec
    real(DP):: sec_ary_dp(imin:imax)
    integer:: i, move_up
  continue
    if ( sclsec == zero_sec .or. factor == 0 ) then
      result = zero_sec
      return
    end if
 
    if ( sclsec % flag_negative ) then
      result % flag_negative = .not. factor < 0
    else
      result % flag_negative = factor < 0
    end if
    factor_abs = abs(factor)
 
    move_up = 0
    sec_ary_dp(:) = 0.0_dp
    do i = imin, imax
      sec_ary_dp(i) = sclsec % sec_ary(i) * factor_abs + move_up
 
      if ( .not. sec_ary_dp(i) < scale_factor ) then
        move_up = int( sec_ary_dp(i) / scale_factor )
        sec_ary_dp(i) = sec_ary_dp(i) - move_up * scale_factor
      else
        move_up = 0
      end if
    end do
 
    if ( move_up /= 0 ) then
      call messagenotify( 'E', 'dc_scaledsec#operator(*)', &
        & 'DC_SCALED_SEC must be smaller than 10^24' )
    end if
 
    result % sec_ary(imin:imax) = nint( sec_ary_dp(imin:imax) )
 

dcscaledsec_mul_sr()

type(dc_scaled_sec) function dc_scaledsec::operator(*)::dcscaledsec_mul_sr	(	type(dc_scaled_sec), intent(in)	sclsec,
		real, intent(in)	factor )

Multiply DC_SCALED_SEC by single precision real

Parameters

[in]	sclsec	DC_SCALED_SEC variable
[in]	factor	Single precision multiplier

Returns

Product of DC_SCALED_SEC and real

Definition at line 1813 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 = sclsec * factor_scl

dcscaledsec_mul_ss()

type(dc_scaled_sec) function dc_scaledsec::operator(*)::dcscaledsec_mul_ss	(	type(dc_scaled_sec), intent(in), target	sclsec1,
		type(dc_scaled_sec), intent(in), target	sclsec2 )

Multiply two DC_SCALED_SEC variables

Parameters

[in]	sclsec1	First DC_SCALED_SEC variable
[in]	sclsec2	Second DC_SCALED_SEC variable

Returns

Product of two DC_SCALED_SEC variables

Definition at line 1598 of file dc_scaledsec.f90.

    use dc_message, only: messagenotify
    implicit none
    type(DC_SCALED_SEC), intent(in), target:: sclsec1, sclsec2
    integer:: sec_ary_int(imin:imax,imin:imax)
!    real(DP):: sec_ary_int(imin:imax,imin:imax)
    integer:: i, j, move_up
    type(DC_SCALED_SEC):: zero_sec
  continue
    if ( sclsec1 == zero_sec .or. sclsec2 == zero_sec ) then
      result = zero_sec
      return
    end if
 
    if ( sclsec1 % flag_negative ) then
      result % flag_negative = .not. sclsec2 % flag_negative
    else
      result % flag_negative = sclsec2 % flag_negative
    end if
 
    move_up = 0
    sec_ary_int(:,:) = 0
    do i = imin, imax
      do j = imin, imax
        sec_ary_int(i,j) = &
          & sclsec1 % sec_ary(j) * sclsec2 % sec_ary(i) + move_up
        if ( i + j > imax .and. sec_ary_int(i,j) /= 0 ) then
          call messagenotify( 'E', 'dc_scaledsec#operator(*)', &
            & 'DC_SCALED_SEC must be smaller than 10^24' )
        end if
        if ( .not. sec_ary_int(i,j) < scale_factor ) then
          move_up = int( sec_ary_int(i,j) / scale_factor_int )
          sec_ary_int(i,j) = sec_ary_int(i,j) - move_up * scale_factor_int
        else
          move_up = 0
        end if
      end do
    end do
 
    result % sec_ary = 0 
    do i = imin, imax
      do j = imin, imax
        if ( i + j < imin ) cycle
        if ( i + j > imax ) cycle
        result % sec_ary(i+j) = result % sec_ary(i+j) + sec_ary_int(i,j)
      end do
    end do
 
    move_up = 0
    do i = imin, imax
      result % sec_ary(i) = result % sec_ary(i) + move_up
      move_up = 0
      do while ( .not. result % sec_ary(i) < scale_factor_int )
        if ( i == imax ) then
          call messagenotify( 'E', 'dc_scaledsec#operator(*)', &
            & 'DC_SCALED_SEC must be smaller than 10^24' )
        end if
        result % sec_ary(i) = result % sec_ary(i) - scale_factor_int
        move_up = move_up + 1
      end do
    end do
 

---

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

• dc_utils/dc_scaledsec.f90