smooth.f90

边界元程序,供力学工作者参考,希望对大家有所帮助

        SUBROUTINE Smooth
        USE CtrlData
        IMPLICIT DOUBLE PRECISION( a-h, o-z )

        iswth = 13
        CALL PreSmooth
        DO ielem = 1, nelem
          CALL ElmOpt( ielem, iswth )
        END DO
        CALL PstSmooth
        RETURN
        END

        SUBROUTINE PreSmooth
        USE CtrlData
        USE GlobData
        IMPLICIT DOUBLE PRECISION( a-h, o-z )
        CALL InitFloat( wsmts, npoin, 0.0D0 )
        CALL InitFloat( vsmts, ndofs, 0.0D0 )
        END

        SUBROUTINE PstSmooth
        USE CtrlData
        USE SolvData
        USE GlobData
        IMPLICIT DOUBLE PRECISION( a-h, o-z )
        DO ipoin = 1, npoin
          weigh = wsmts( ipoin )
          DO idofn = 1, mdofn
            IF( lsmts( idofn ) .NE. 0 ) THEN
              idofs = ( ipoin - 1 ) * mdofn + idofn
              IF( weigh .GT. 1.0D-8 )                                          &
              vslvs( idofs ) = vsmts( idofs ) / weigh
            ENDIF
          END DO
        END DO
        END