6. Solving the finite element problem in Florence

The final step is to solve our problem. To do this we need to create a FEMSolver object

solver = FEMSolver()

and then call its Solve method by specifying all our previously defined data

solution = solver.Solve(material=material, boundary_condition=boundary_condition, 
    formulation=formulation, mesh=mesh)

This will trigger the solver in Florence. Once the analysis is complete we can write the results (visualisation of low and high order elements is supported) to VTK file(s) to visualise them using paraview

solution.WriteVTK("my_results")

or alternatively, using Florence's built-in plotter (visualisation of low and high order elements is supported)

solution.Plot(configuration="deformed", quantity=0)