def METHOD_NAME(self):
    self.comm = Epetra.PyComm()
    self.dataname = "testPyMatrix%d.dat" % comm.MyPID()
    myPID = self.comm.MyPID()
    # Create Spaces
    n = 10 * self.comm.NumProc()
    globalRows = list(range(myPID*10, (myPID+1)*10))
    rowSpace = ML.Space(n, globalRows)
    globalCols = globalRows[:]
    if myPID > 0:
        globalCols.insert(0,myPID*10-1)
    if myPID < self.comm.NumProc()-1:
        globalCols.append((myPID+1)*10)
    colSpace = ML.Space(-1, globalCols)
    # Create and fill matrix
    self.matrix = ML.PyMatrix(rowSpace, colSpace)
    for i in globalRows:
        if i > 0:
            self.matrix[i,i-1] =  1
        self.matrix[i,i] = -2
        if i < self.comm.NumProc()*10-1:
            self.matrix[i,i+1] =  1
    self.matrix.FillComplete()
    self.comm.Barrier()