def METHOD_NAME(casefile, init_feelpp):
    e = init_feelpp
    feelpp.Environment.setConfigFile(casefile)
    model_path = "$cfgdir/"+os.path.splitext(os.path.basename(casefile))[0] + ".json"
    model_properties = CRBModelProperties(worldComm=feelpp.Environment.worldCommPtr())
    model_properties.setup(model_path)
    modelParameters = model_properties.parameters()
    Dmu = feelpp.mor._mor.ParameterSpace.New(modelParameters, feelpp.Environment.worldCommPtr())
    print(feelpp.mor.__file__)
    Nsamples = 10
    s = Dmu.sampling()
    s.sample(Nsamples, samplingMode="random")
    assert len(s) == Nsamples , "wrong number of samples"
    v = s.getVector()
    assert len(v) == Nsamples , "wrong number of samples"
    s.writeOnFile("samplingTest.sample")
    s_read  = Dmu.sampling()
    N_read = s_read.readFromFile("samplingTest.sample")
    v_read  = s_read.getVector()
    assert N_read == Nsamples, "wrong number of samples read"
    for i in range(Nsamples):
        for p in Dmu.parameterNames():
            diff = abs(v[i].parameterNamed(p) - v_read [i].parameterNamed(p))
            assert diff < 1e-10 , f"wrong sample read at index {i} : {v[i].parameterNamed(p)} != {v_read[i].parameterNamed(p)}"