sne2.vb

一个lbm的程序

Imports System
Public Class Sne2

    Friend Shared Function getState(ByVal i As Int32, ByVal j As Int32) As Byte
        Dim o As Byte
        Dim fluid, empty As Boolean
        For o = 1 To q
            Select Case flag(i + ei(o, 0), j + ei(o, 1))
                Case CONST_FLUID
                    fluid = True
                Case CONST_EMPTY
                    empty = True
            End Select
        Next
        If Not fluid Then
            Return 0
        ElseIf Not empty Then
            Return 2
        Else
            Return 1
        End If
    End Function

    Friend Shared Function getCellNum(ByVal i As Int32, ByVal j As Int32, ByVal celltype As Byte) As Byte
        Dim c, o As Byte
        For o = 1 To q
            If flag(i + ei(o, 0), j + ei(o, 1)) = celltype Then c += 1
        Next
        Return c
    End Function

    Friend Shared Function getNormalCellNum(ByVal i As Int32, ByVal j As Int32) As Byte()
        Dim c(1), o As Byte
        For o = 1 To q
            Select Case flag(i + ei(o, 0), j + ei(o, 1))
                Case CONST_EMPTY
                    c(0) += 1
                Case CONST_FLUID
                    c(1) += 1
            End Select
        Next
        Return c
    End Function

    Friend Shared Function isIn(ByVal i As Int32, ByVal j As Int32) As Boolean
        If i >= 0 And i <= m And j >= 0 And j <= n Then
            Return True
        Else
            Return False
        End If
    End Function

    Friend Shared Function hasFluid(ByVal i As Int32, ByVal j As Int32) As Boolean
        Dim o As Byte
        For o = 1 To q
            If flag(i + ei(o, 0), j + ei(o, 1)) = CONST_FLUID Then Return True
        Next
        Return False
    End Function

    Friend Shared Function hasCellType(ByVal i As Int32, ByVal j As Int32, ByVal celltype As Byte) As Boolean
        Dim o As Byte
        For o = 1 To q
            If flag(i + ei(o, 0), j + ei(o, 1)) = celltype Then Return True
        Next
        Return False
    End Function

    Friend Shared Function getNormalDirection(ByVal i As Int32, ByVal j As Int32) As Double()
        Dim k As Byte
        Dim nd(1) As Double
        For k = 1 To q
            nd(0) += s(k) * ei(k, 0) * eps(i + ei(k, 0), j + ei(k, 1))
            nd(1) += s(k) * ei(k, 1) * eps(i + ei(k, 0), j + ei(k, 1))
        Next
        Return nd
    End Function
End Class