What is the logic behind deviding arrays of of double precision matrix to a large number before sending that matrix to a solver? : I am given a code which divides the members of [A] in [A]{X}=[B] to 10^4 when assembling matrix [A]. 
Then it also divides the results array to this number to calculate correct value.
I cannot understand why this should be done? Does double precision has limitation on numbers of integer it can hold; so the maker of this code wanted to increase number of float digits? Or maybe he misunderstood the conception of double precision?