  In the framework of various statistical models as well as of mechanisms for
color confinement, disorder parameters can be developed which are generally
expressed as ratios of partition functions and whose numerical determination is
usually challenging. We develop an efficient method for their computation and
apply it to the study of dual superconductivity in 4d compact U(1) gauge
theory.
