A numerical method for separation of stresses in photo-orthotropic elasticity