A sample-based method for computing the radiosity inverse matrix