Array calculations for all functions

[scotthavens]

Most of the equations are capable of handling array multiplication. For example, the cosine of the illumination angle in SRMF can be calculated over the whole domain. Passing that to the shade function will probably not work as it is now.