def setUp(self):

        cs = 10.
        ncx, ncy, ncz = 10, 10, 10
        npad = 4
        freq = 1e2

        hx = [(cs,npad,-1.3), (cs,ncx), (cs,npad,1.3)]
        hy = [(cs,npad,-1.3), (cs,ncy), (cs,npad,1.3)]
        hz = [(cs,npad,-1.3), (cs,ncz), (cs,npad,1.3)]
        mesh = Mesh.TensorMesh([hx,hy,hz], 'CCC')

        mapping = Maps.ExpMap(mesh)

        x = np.linspace(-10,10,5)
        XYZ = Utils.ndgrid(x,np.r_[0],np.r_[0])
        rxList = EM.FDEM.Rx(XYZ, 'exi')
        Src0 = EM.FDEM.Src.MagDipole([rxList],loc=np.r_[0.,0.,0.], freq=freq)

        survey = EM.FDEM.Survey([Src0])

        prb = EM.FDEM.Problem_b(mesh, mapping=mapping)
        prb.pair(survey)

        try:
            from pymatsolver import MumpsSolver
            prb.Solver = MumpsSolver
        except __HOLE__, e:
            prb.Solver = SolverLU

        sig = 1e-1
        sigma = np.ones(mesh.nC)*sig
        sigma[mesh.gridCC[:,2] > 0] = 1e-8
        m = np.log(sigma)

        self.prb = prb
        self.mesh = mesh
        self.m = m
        self.Src0 = Src0
        self.sig = sig
dataset/ETHPy150Open simpeg/simpeg/tests/em/fdem/forward/test_FDEM_analytics.py/FDEM_analyticTests.setUp