mbw.f90

来自「这是一个用fortran编写的“第一类边界条件的有限元算法”」· F90 代码 · 共 25 行

!***************************************************************************!
!                       SUBROUTINE MBW(ne,i3,iw)                            !
!                                                                           !
!   功能    计算总体系数矩阵的半带宽。                                      !
!   参数说明                                                                !
!     ne---整型变量,输入参数,单元总数。                                     !
!     i3---3*ne的二维整数组,输入参数,存放单元节点编号。                     !
!     iw---整型变量,输出参数,半带宽。                                       !
!***************************************************************************!

SUBROUTINE MBW(ne,i3,iw)

   DIMENSION i3(3,iw)

   iw=0
   do i=1,ne
     m=max(iabs(i3(1,i)-i3(2,i)),iabs(i3(2,i)-i3(3,i)),iabs(i3(3,i)-i3(1,i)))
	        
	   if (m+1.gt.iw) then
	      ie=m+1
       end if
	end do
	return
end subroutine MBW