assert (abs(r_center.direction*Vector([0,0,1])-r_center.direction.norm())<0.01,"c.ray_at(c.image_nrows / 2, c.image_ncols / 2) should have the direction [0,0,1].")
assertabs(r_center.direction*Vector([0,0,1])-r_center.direction.norm())<0.01,"c.ray_at(c.image_nrows / 2, c.image_ncols / 2) should have the direction [0,0,1]."