Source code for sfepy.solvers.ls_mumps_parallel

from __future__ import absolute_import

import numpy as nm
from mpi4py import MPI
import sfepy.solvers.ls_mumps as mumps
from tempfile import gettempdir
import os.path as op




def tmpfile(fname):
    return op.join(gettempdir(), fname)





def mumps_parallel_solve():
    comm = MPI.COMM_WORLD
    flags = nm.memmap(tmpfile('vals_flags.array'), dtype='int32', mode='r')

    mumps_ls = mumps.MumpsSolver(system={0: 'real', 1: 'complex'}[flags[1]],
                                 is_sym=flags[2])

    if flags[3]:
        mumps_ls.set_verbose()

    if comm.rank == 0:
        dtype = {0: 'float64', 1: 'complex128'}[flags[1]]

        vals_mtx = nm.memmap(tmpfile('vals_mtx.array'), dtype=dtype, mode='r')
        vals_rhs = nm.memmap(tmpfile('vals_rhs.array'), dtype=dtype, mode='r')
        idxs = nm.memmap(tmpfile('idxs.array'), dtype='int32',
                         mode='r').reshape((2, vals_mtx.shape[0]))

        mumps_ls.set_rcd_centralized(idxs[0, :], idxs[1, :], vals_mtx,
                                     flags[0])
        x = vals_rhs.copy()
        mumps_ls.set_rhs(x)

    mumps_ls(6)  # analyse, factorize, solve

    if comm.rank == 0:
        vals_x = nm.memmap(tmpfile('vals_x.array'), dtype=dtype, mode='w+',
                           shape=x.shape)
        vals_x[:] = x

    del(mumps_ls)



if __name__ == '__main__':
    comm = MPI.Comm.Get_parent()
    mumps_parallel_solve()
    comm.Disconnect()