A Boundary Element-Response Matrix method for the solution of 3D criticality problems