estim_flops.F File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	mumps_estim_flops (inode, n, procnode_steps, keep199, nd, fils, frere_steps, step, pimaster, keep28, keep50, keep253, flop1, iw, liw, xsize)
subroutine	mumps_update_flops_root (opeliw, keep50, nfront, npiv, nprow, npcol, myid)
subroutine	mumps_get_flops_cost (nfront, npiv, nass, keep50, level, cost)
subroutine	mumps_print_still_active (myid, keep, dkeep17, opeliw, oplast_printed, mpa)

Function/Subroutine Documentation

mumps_estim_flops()

subroutine mumps_estim_flops	(	integer	inode,
		integer	n,
		integer, dimension(keep28)	procnode_steps,
		integer	keep199,
		integer, dimension(keep28)	nd,
		integer, dimension(n)	fils,
		integer, dimension(keep28)	frere_steps,
		integer, dimension(n)	step,
		integer, dimension(keep28)	pimaster,
		integer	keep28,
		integer	keep50,
		integer	keep253,
		double precision	flop1,
		integer, dimension( liw )	iw,
		integer	liw,
		integer	xsize )

Definition at line 14 of file estim_flops.F.

      IMPLICIT NONE
      INTEGER INODE, N, KEEP50, LIW, KEEP199, KEEP28, KEEP253
      INTEGER PROCNODE_STEPS(KEEP28), ND(KEEP28),
     &        FILS(N), FRERE_STEPS(KEEP28),
     &        STEP(N), 
     & PIMASTER(KEEP28),
     &  IW( LIW )
      INTEGER XSIZE
      DOUBLE PRECISION FLOP1
      INTEGER NUMORG, IN, NASS, IFSON, NUMSTK, NFRONT, NPIV, NCB, 
     &        LEVEL, ISON
      LOGICAL MUMPS_IN_OR_ROOT_SSARBR
      INTEGER MUMPS_TYPENODE
      EXTERNAL mumps_in_or_root_ssarbr, mumps_typenode
      include 'mumps_headers.h'
      FLOP1 = 0.0D0
      IF (MUMPS_IN_OR_ROOT_SSARBR(PROCNODE_STEPS(STEP(INODE)),
     &                KEEP199) ) RETURN
      IN     = INODE
      NUMORG = 0
   10 NUMORG = NUMORG + 1
      IN = FILS(IN)
.GT.      IF (IN  0) GOTO 10
      NUMSTK = 0
      NASS = 0
      IFSON = -IN
      ISON = IFSON
.EQ.      IF (ISON  0) GOTO 30
   20 NUMSTK = NUMSTK + 1
      NASS = NASS + IW(PIMASTER(STEP(ISON)) + 1 +XSIZE)
      ISON = FRERE_STEPS(STEP(ISON))
.GT.      IF (ISON  0) GOTO 20
   30 NFRONT = ND(STEP(INODE)) + NASS + KEEP253
      NPIV  = NASS + NUMORG
      NCB   = NFRONT - NPIV
      LEVEL = MUMPS_TYPENODE(PROCNODE_STEPS(STEP(INODE)),KEEP199)
      CALL MUMPS_GET_FLOPS_COST(NFRONT,NPIV,NPIV,KEEP50,LEVEL,FLOP1)
      RETURN

mumps_get_flops_cost()

subroutine mumps_get_flops_cost	(	integer, intent(in)	nfront,
		integer, intent(in)	npiv,
		integer, intent(in)	nass,
		integer, intent(in)	keep50,
		integer, intent(in)	level,
		double precision, intent(out)	cost )

Definition at line 72 of file estim_flops.F.

      IMPLICIT NONE
      INTEGER, intent(in) :: NFRONT,NPIV,KEEP50,LEVEL, NASS
      DOUBLE PRECISION, intent(out) :: COST
      IF (keep50.EQ.0) THEN
        IF (level.EQ.1 .OR. level.EQ.3) THEN
          cost = dble(2) * dble(nfront) * dble(npiv) *
     &      dble(nfront - npiv - 1) +
     &      dble(npiv) * dble(npiv + 1) * dble(2 * npiv + 1)
     &          / dble(3)
          cost = cost + dble(2 * nfront - npiv - 1)
     &      * dble(npiv) /dble(2)
        ELSEIF (level.EQ.2) THEN
          cost = dble(2*nass)*dble(nfront) -
     &          dble(nass+nfront)*dble(npiv+1)
          cost = dble(npiv)*cost +
     &     dble(2 * nass - npiv - 1) * dble(npiv) / dble(2) +
     &     dble(npiv) * dble(npiv + 1) *
     &     dble(2 * npiv + 1) /dble(3)
        ENDIF
      ELSE
        IF (level.EQ.1 .OR. (level.EQ.3 .AND. keep50.EQ.1)) THEN
          cost = dble(npiv) * (
     &          dble( nfront ) * dble( nfront ) +
     &          dble( nfront ) - (
     &          dble( nfront)*dble(npiv) + dble(npiv+1)
     &          )) +( dble(npiv)*dble(npiv+1)
     &          *dble(2*npiv+1))/ dble(6)
        ELSE IF (level.EQ.3.AND.keep50.EQ.2) THEN
          cost = dble(2) * dble(nfront) * dble(npiv) *
     &      dble(nfront - npiv - 1) +
     &      dble(npiv) * dble(npiv + 1) *
     &      dble(2 * npiv + 1) / dble(3)
          cost = cost + dble(2 * nfront - npiv - 1)
     &         * dble(npiv) / dble(2)
        ELSE
          cost = dble(npiv) * (
     &          dble( nass ) * dble( nass ) + dble( nass )
     &        - ( dble( nass) * dble(npiv) + dble( npiv + 1 ) ) )
     &        + ( dble(npiv)*dble(npiv+1)*dble(2*npiv+1) )
     &        / dble( 6 )
        ENDIF
      ENDIF
      RETURN

mumps_print_still_active()

subroutine mumps_print_still_active	(	integer, intent(in)	myid,
		integer, dimension (500), intent(in)	keep,
		double precision	dkeep17,
		double precision, intent(in)	opeliw,
		double precision, intent(inout)	oplast_printed,
		integer, intent(in)	mpa )

Definition at line 118 of file estim_flops.F.

      IMPLICIT NONE
      INTEGER, intent(in)             :: KEEP (500), MYID, MPA
      DOUBLE PRECISION                :: DKEEP17
      DOUBLE PRECISION, intent(in)    :: OPELIW
      DOUBLE PRECISION, intent(inout) :: OPLAST_PRINTED
      IF (mpa.GT.0) THEN
        IF ( (opeliw-oplast_printed).GT. dkeep17) THEN
         WRITE(mpa,'(A,I6,A,A,1PD10.3)')
     &   ' ... MPI process', myid, 
     &   ': theoretical number of flops locally performed',
     &   ' so far        = ',
     &    opeliw
         oplast_printed = opeliw
        ENDIF
      ENDIF
      RETURN

mumps_update_flops_root()

subroutine mumps_update_flops_root	(	double precision, intent(inout)	opeliw,
		integer, intent(in)	keep50,
		integer, intent(in)	nfront,
		integer, intent(in)	npiv,
		integer, intent(in)	nprow,
		integer, intent(in)	npcol,
		integer, intent(in)	myid )

Definition at line 59 of file estim_flops.F.

      DOUBLE PRECISION, intent(inout) :: OPELIW
      INTEGER, intent(in) :: KEEP50, NFRONT, NPIV,
     &           NPROW, NPCOL, MYID
      DOUBLE PRECISION :: COST, COST_PER_PROC
      INTEGER, PARAMETER :: LEVEL3 = 3
      CALL mumps_get_flops_cost(nfront, npiv, nfront, keep50, level3, 
     &                          cost)
      cost_per_proc = dble(int( cost,8) / int(nprow * npcol,8))
      opeliw = opeliw + cost_per_proc
      RETURN