I have written a code for solving Laplace’s equation which works well on some grids but not so well when there is a short distance from an element to another node point. I would like to try increasing the order of the Gaussian quadrature scheme. How do I do this?

This prints the quadrature order for regular and singular integrals. Bempp-cl does not make a distinction between near, medium, and far interactions anymore.
Notice that these are the order of quadrature, not the number of points. This information can be accessed by

By default, all operators take the global parameters when assembling the matrix. So set this before defining the operators.

If op is an operator, you can check the quadrature order by

print(op.parameters.quadrature.regular)

If you need to define different operators with different quadrature order in the same script, you can pass an operator object to the parameters argument of an operator.