Esempio n. 1
0
    def _explicitBuildMatrix(self, SparseMatrix, oldArray, id1, id2, b, weight, mesh, boundaryConditions, interiorFaces, dt):

        coeffMatrix = self._getCoeffMatrix(mesh, weight)

        inline._optionalInline(self._explicitBuildMatrixIn, self._explicitBuildMatrixPy, oldArray, id1, id2, b, coeffMatrix, mesh, interiorFaces, dt, weight)

        N = mesh.getNumberOfCells()
        M = mesh._getMaxFacesPerCell()

        for boundaryCondition in boundaryConditions:

            LL,bb = boundaryCondition._buildMatrix(SparseMatrix, N, M, coeffMatrix)
            if LL != 0:
##              b -= LL.takeDiagonal() * numerix.array(oldArray)
                b -= LL * numerix.array(oldArray)
            b += bb
Esempio n. 2
0
    def _buildMatrix(self, var, SparseMatrix, boundaryConditions=(), dt=1., equation=None):
        N = len(var)
        b = numerix.zeros((N),'d')
        L = SparseMatrix(size=N)
        
        # The sign of the matrix diagonal doesn't seem likely to change
        # after initialization, but who knows?
        if equation is not None:
            from fipy.tools.numerix import sign, add
            self._diagonalSign.setValue(sign(add.reduce(equation.matrix.takeDiagonal())))
        else:
            self._diagonalSign.setValue(1)
            
        coeffVectors = self._getCoeffVectors(var=var)

        inline._optionalInline(self._buildMatrixIn, self._buildMatrixPy, L, var.getOld(), b, dt, coeffVectors)
        
        return (L, b)
Esempio n. 3
0
    def _calcValue(self):
        N = self.mesh.getNumberOfCells()
        M = self.mesh._getMaxFacesPerCell()

        ids = self.mesh._getCellFaceIDs()

        orientations = self.mesh._getCellFaceOrientations()
        volumes = self.mesh.getCellVolumes()

        return inline._optionalInline(self._calcValueIn, self._calcValuePy, N, M, ids, orientations, volumes)
Esempio n. 4
0
File: grid2D.py Progetto: regmi/fipy
 def _createCells(self):
     """
     cells = (f1, f2, f3, f4) going anticlock wise.
     f1 etc. refer to the faces
     """
     return inline._optionalInline(self._createCellsIn, self._createCellsPy)
        def _calcValue(self):
            P  = self.P.getNumericValue()

            return inline._optionalInline(self._calcValueIn, self._calcValuePy, P)
Esempio n. 6
0
 def _calcValue(self):
     alpha = self.mesh._getFaceToCellDistanceRatio()
     id1, id2 = self.mesh._getAdjacentCellIDs()
     return inline._optionalInline(self._calcValueIn, self._calcValuePy, alpha, id1, id2)
    def _calcValue(self):
        from fipy.tools import inline

        return inline._optionalInline(self._calcValueInline, self._calcValuePy)
Esempio n. 8
0
 def _calcValue(self):        
     return inline._optionalInline(self._calcValueInline, self._calcValuePy)
Esempio n. 9
0
def putAdd(vector, ids, additionVector):
    """ This is a temporary replacement for Numeric.put as it was not doing
    what we thought it was doing.
    """
    from fipy.tools import inline
    inline._optionalInline(_putAddIn, _putAddPy, vector, ids, additionVector)
Esempio n. 10
0
 def _calcValue(self):
     if not self.canInline:
         return self._calcValuePy()
     else:
         from fipy.tools import inline
         return inline._optionalInline(self._calcValueIn, self._calcValuePy)
Esempio n. 11
0
 def _getAreaProjections(self):
     return inline._optionalInline(self._getAreaProjectionsIn, self._getAreaProjectionsPy)
Esempio n. 12
0
 def getFaceCellIDs(self):
     return inline._optionalInline(self._getFaceCellIDsIn, self._getFaceCellIDsPy)
Esempio n. 13
0
 def _getAdjacentCellIDs(self):
     return inline._optionalInline(self._getAdjacentCellIDsIn, self._getAdjacentCellIDsPy)