How to specify different linear solver (MUMPS) in C++

Question

What's the C++ equivalent of this:

solve(a == L, u, bc, solver_parameters={'linear_solver':'mumps'})

When I write it as above, I get errors.

Answer 1

problem = LinearVariationalProblem(a, L, u, bc)
solver = LinearVariationalSolver(problem)
solver.parameters["linear_solver"] = "mumps"
solver.solve

Comments

You're probably meaning

solver.solve()

on the last line :)