fac_asm_build_sort_index_m.F File Reference

Go to the source code of this file.

Modules
module	mumps_build_sort_index_m

Functions/Subroutines
subroutine	mumps_build_sort_index_m::mumps_build_sort_index (myid, inode, n, ioldps, hf, lp, lpok, nfront, nfront_eff, perm, dad, nass1, nass, numstk, numorg, iwposcb, iwpos, ifson, step, pimaster, ptrist, ptraiw, iw, liw, intarr, lintarr, itloc, fils, frere_steps, son_level2, niv1, keep, keep8, iflag, ison_in_place, procnode_steps, slavef, sonrows_per_row, lsonrows_per_row, mumps_tps_arr, l0_omp_mapping)
subroutine	mumps_sort (n, perm, iw, liw)
subroutine	mumps_sorted_merge (n, nass1, perm, itloc, small, lsmall, large, llarge, merge, lmerge)

Function/Subroutine Documentation

mumps_sort()

subroutine mumps_sort	(	integer	n,
		integer, dimension( n )	perm,
		integer, dimension( liw )	iw,
		integer	liw )

Definition at line 626 of file fac_asm_build_sort_index_m.F.

      IMPLICIT NONE
      INTEGER N, LIW
      INTEGER PERM( N ), IW( LIW )
      INTEGER I, SWAP
      LOGICAL DONE
      done = .false.
      DO WHILE ( .NOT. done )
        done = .true.
        DO i = 1, liw - 1
          IF ( perm( iw( i ) ) .GT. perm( iw( i + 1 ) ) ) THEN
            done = .false.
            swap  = iw( i + 1 )
            iw( i + 1 ) = iw( i )
            iw( i ) = swap
          END IF
        END DO
      END DO
      RETURN

mumps_sorted_merge()

subroutine mumps_sorted_merge	(	integer	n,
		integer	nass1,
		integer, dimension( n )	perm,
		integer, dimension( n )	itloc,
		integer, dimension(lsmall)	small,
		integer	lsmall,
		integer, dimension(llarge)	large,
		integer	llarge,
		integer, dimension(lmerge)	merge,
		integer	lmerge )

Definition at line 646 of file fac_asm_build_sort_index_m.F.

      IMPLICIT NONE
      INTEGER N, NASS1, LSMALL, LLARGE, LMERGE
      INTEGER PERM( N ), ITLOC( N ) 
      INTEGER SMALL(LSMALL), LARGE(LLARGE), MERGE(LMERGE)
      INTEGER PSMALL, PLARGE, PMERGE, VSMALL, VLARGE, VMERGE
      psmall = 1
      plarge = 1
      pmerge = 1
      DO WHILE ( psmall .LE. lsmall .or. plarge.LE. llarge )
        IF ( psmall .GT. lsmall ) THEN
          vmerge = large( plarge )
          plarge = plarge + 1
        ELSE IF ( plarge .GT. llarge ) THEN
          vmerge = small( psmall )
          psmall = psmall + 1
        ELSE
          vsmall = small( psmall )
          vlarge = large( plarge )
          IF ( perm( vsmall ) .LT. perm( vlarge ) ) THEN
            vmerge = vsmall
            psmall   = psmall + 1
          ELSE
            vmerge = vlarge
            plarge   = plarge + 1
          END IF
        END IF
        merge( pmerge ) = vmerge
        itloc( vmerge ) = pmerge + nass1
        pmerge = pmerge + 1
      END DO
      pmerge = pmerge - 1
      RETURN