compact_real_string Subroutine

  1. penf.F90
  2. penf_stringify
  3. compact_real_string

private pure subroutine compact_real_string(string)

Arguments

Type IntentOptional AttributesName
character(len=*), intent(inout) :: string

string representation of a real number.

Description

Compact a string representing a real number, so that the same value is displayed with fewer characters.

Called By

proc~~compact_real_string~~CalledByGraph proc~compact_real_string compact_real_string proc~str_r4p str_R4P proc~str_r4p->proc~compact_real_string proc~str_r8p str_R8P proc~str_r8p->proc~compact_real_string proc~str_r16p str_R16P proc~str_r16p->proc~compact_real_string proc~str_a_r4p str_a_R4P proc~str_a_r4p->proc~str_r4p interface~str str interface~str->proc~str_r4p interface~str->proc~str_r8p interface~str->proc~str_a_r4p proc~str_a_r8p str_a_R8P interface~str->proc~str_a_r8p proc~bctor_r8p bctor_R8P proc~bctor_r8p->interface~str proc~string_assign_integer_i8p string_assign_integer_I8P proc~string_assign_integer_i8p->interface~str proc~bctoi_i2p bctoi_I2P proc~bctoi_i2p->interface~str proc~string_assign_integer_i4p string_assign_integer_I4P proc~string_assign_integer_i4p->interface~str proc~bctor_r4p bctor_R4P proc~bctor_r4p->interface~str proc~bctoi_i1p bctoi_I1P proc~bctoi_i1p->interface~str proc~autotest autotest proc~autotest->interface~str proc~string_assign_integer_i1p string_assign_integer_I1P proc~string_assign_integer_i1p->interface~str proc~bctoi_i8p bctoi_I8P proc~bctoi_i8p->interface~str proc~string_assign_real_r4p string_assign_real_R4P proc~string_assign_real_r4p->interface~str proc~bctoi_i4p bctoi_I4P proc~bctoi_i4p->interface~str proc~string_assign_integer_i2p string_assign_integer_I2P proc~string_assign_integer_i2p->interface~str proc~bctor_r16p bctor_R16P proc~bctor_r16p->interface~str proc~string_assign_real_r16p string_assign_real_R16P proc~string_assign_real_r16p->interface~str proc~string_assign_real_r8p string_assign_real_R8P proc~string_assign_real_r8p->interface~str interface~bcton bcton interface~bcton->proc~bctor_r8p interface~bcton->proc~bctoi_i2p interface~bcton->proc~bctor_r4p interface~bcton->proc~bctoi_i1p interface~bcton->proc~bctoi_i8p interface~bcton->proc~bctoi_i4p proc~str_a_r8p->proc~str_r8p proc~str_a_r16p str_a_R16P proc~str_a_r16p->proc~str_r16p
Help

Variables

TypeVisibility AttributesNameInitial
character(len=len(string)), public :: significand

Significand characters.
character(len=len(string)), public :: expnt

Exponent characters.
character(len=2), public :: separator

Separator characters.
integer(kind=I4P), public :: exp_start

Start position of exponent.
integer(kind=I4P), public :: decimal_pos

Decimal positions.
integer(kind=I4P), public :: sig_trim

Signature trim.
integer(kind=I4P), public :: exp_trim

Exponent trim.
integer(kind=I4P), public :: i

counter

Source Code

  pure subroutine compact_real_string(string)
  !---------------------------------------------------------------------------------------------------------------------------------
  !< author: Izaak Beekman
  !< date: 02/24/2015
  !<
  !< Compact a string representing a real number, so that the same value is displayed with fewer characters.
  !---------------------------------------------------------------------------------------------------------------------------------
  character(len=*),intent(inout) :: string      !< string representation of a real number.
  character(len=len(string))     :: significand !< Significand characters.
  character(len=len(string))     :: expnt       !< Exponent characters.
  character(len=2)               :: separator   !< Separator characters.
  integer(I4P)                   :: exp_start   !< Start position of exponent.
  integer(I4P)                   :: decimal_pos !< Decimal positions.
  integer(I4P)                   :: sig_trim    !< Signature trim.
  integer(I4P)                   :: exp_trim    !< Exponent trim.
  integer(I4P)                   :: i           !< counter
  !---------------------------------------------------------------------------------------------------------------------------------

  !---------------------------------------------------------------------------------------------------------------------------------
  string = adjustl(string)
  exp_start = scan(string, 'eEdD')
  if (exp_start == 0) exp_start = scan(string, '-+', back=.true.)
  decimal_pos = scan(string, '.')
  if (exp_start /= 0) separator = string(exp_start:exp_start)
  if ( exp_start < decimal_pos ) then ! possibly signed, exponent-less float
    significand = string
    sig_trim = len(trim(significand))
    do i = len(trim(significand)), decimal_pos+2, -1 ! look from right to left at 0s, but save one after the decimal place
      if (significand(i:i) == '0') then
        sig_trim = i-1
      else
        exit
      endif
    enddo
    string = trim(significand(1:sig_trim))
  elseif (exp_start > decimal_pos) then ! float has exponent
    significand = string(1:exp_start-1)
    sig_trim = len(trim(significand))
    do i = len(trim(significand)),decimal_pos+2,-1 ! look from right to left at 0s
      if (significand(i:i) == '0') then
        sig_trim = i-1
      else
        exit
      endif
    enddo
    expnt = adjustl(string(exp_start+1:))
    if (expnt(1:1) == '+' .or. expnt(1:1) == '-') then
      separator = trim(adjustl(separator))//expnt(1:1)
      exp_start = exp_start + 1
      expnt     = adjustl(string(exp_start+1:))
    endif
    exp_trim = 1
    do i = 1,(len(trim(expnt))-1) ! look at exponent leading zeros saving last
      if (expnt(i:i) == '0') then
        exp_trim = i+1
      else
        exit
      endif
    enddo
    string = trim(adjustl(significand(1:sig_trim)))// &
             trim(adjustl(separator))// &
             trim(adjustl(expnt(exp_trim:)))
  !else ! mal-formed real, BUT this code should be unreachable
  endif
  !---------------------------------------------------------------------------------------------------------------------------------
  endsubroutine compact_real_string