For any given field, you can compute the residue of the zeta function to arbitrary precision: just take the product of sufficiently many Euler factors. In practice, this can be quite slow for high degree fields, since you might need lots of factors to get the desired precision.
Of course, computations can never prove the class number 1 problem, since a computation can only check finitely many fields.
The main interest of the formula is that it provides a prototypical example of a "local-global" principle: the left hand side is something that you compute locally at each prime, while the right hand side is something manifestly global. As you probably know, the class number formula has served as a blue print for very deep and far reaching generalisations, such as the Birch and Swinnerton-Dyer conjecture for example. Most of these generalisations are wide open.
Since you also asked for references, e.g. about extracting information about fields from relations between zeta functions, here is a paper of mine in which concrete information about the Galois module structure of the units is extracted from the relationship you mention.
Edit: I failed to mention perhaps the most widely used practical application of the class number formula: when you ask MAGMA or most other algebraic number theory computer packages to compute the class group of a number field, the routines will make use of the class number formula. The way this works is as follows: it is easy to find generators of the class group (e.g. using the Minkowski bound). You then want to look for all possible relations, but how do you know when you are done? Of course, you can prove by extremely brute force that you have found all the relations, but the way the computer will do this, which is much faster, is to also find a finite index subgroup of the unit group, and as you are finding more relations for the ideal classes, you also saturate your finite index subgroup of the units, i.e. keep checking whether certain elements of your subgroup are divisible by an integer inside the unit group, thus reducing the index of your subgroup inside the full unit group. More relations lead to a smaller candidate for the class group, and further saturation leads to a smaller regulator. In the meantime, you compute the residue of the zeta function to some precision, thus knowing in advance what the product of class number and regulator should be. When you reach the expected product, you stop. Thus, most routines for computing the class group secretly also compute the unit group of the number field at the same time, and the analytic class number formula provides a stopping condition for this calculation.