test script now with scipy solver

    #switch to fuse simulation to make it easier to read the resulting stiffness matrix

    dict_of_simulation_parameters = {"niter":10_000_000, "externalV":1.0, "accuracy":12}

    builder = mn.sim.FuseDirichletBuilder()

    solver = mn.sim.PARDISOBindingsSolver()

    solver = mn.sim.ScipySolver()

    intrescal = mn.sim.StartInterimResultsCalculator() #InterimResultsCalculator stack pre-computes data needed for output and simulation step

    intrescal = mn.sim.FuseCurrentsVoltagesCalculator(intrescal)

    outgen = mn.sim.StartOutputGenerator()

    outgen = mn.sim.IVCurveGenerator(outgen)

    applier = mn.sim.HottestFuseBreakApplier()

    checker = mn.sim.NonInitialChecker(mn.sim.ConnectedPathChecker())

    simulation = mn.sim.Simulation(builder, solver, intrescal, outgen, applier, checker)

    simulation = mn.sim.RunTillStopconditionSimulation(builder, solver, intrescal, outgen, applier, checker)

    dict_of_crs_edge_data = bounded_network.get_network_data() 

    simdata = simulation.start_simulation(dict_of_crs_edge_data, dict_of_simulation_parameters, maximum_runtime=1)

    simdata = simulation.start_simulation(dict_of_crs_edge_data, dict_of_simulation_parameters, maximum_runtime=3600)

    simdata_dict = simdata.get_simdata_dict()

    network_descriptor_dict = bounded_network.get_construction_data()

    mn.datman.save_simulation_output(network_descriptor_dict, simdata_dict, simulation.simulationname, simulation.simulationdoc, targetdirectory=destination, suffix=str(seed), overwritemode = False)

    mn.datman.save_simulation_output(network_descriptor_dict, simdata_dict, simulation.simulationname, simulation.simulationdoc,

                                     targetdirectory=destination, suffix=str(seed), overwritemode = False)

    print("done with sim of seed", seed)

    return (network_descriptor_dict, simdata_dict, simulation.simulationname, simulation.simulationdoc)

    offset = 2

    for seed in range(1000, 1001):

            list_of_seeds.append(seed) 

    crack_lengths = [1, 2]

    crack_lengths = [1]

    for crack_length in crack_lengths: 

        destination = f"/FASTTEMP/p7/lpyka/hierarchical_interface/dummy_test/test/a_{crack_length}/"

        run_a_simulation_fuse(crack_length=crack_length, offset=offset, destination=destination, seed=1000)