212 #define CN(I,J,K) cn[(I)*uinc+(J)*dim+(K)]  in _math_de_casteljau_surf()  macro
218 	    du[k] = vs * (CN(1, 0, k) - CN(0, 0, k)) +  in _math_de_casteljau_surf()
219 	       v * (CN(1, 1, k) - CN(0, 1, k));  in _math_de_casteljau_surf()
222 	    dv[k] = us * (CN(0, 1, k) - CN(0, 0, k)) +  in _math_de_casteljau_surf()
223 	       u * (CN(1, 1, k) - CN(1, 0, k));  in _math_de_casteljau_surf()
226 	    out[k] = us * (vs * CN(0, 0, k) + v * CN(0, 1, k)) +  in _math_de_casteljau_surf()
227 	       u * (vs * CN(1, 0, k) + v * CN(1, 1, k));  in _math_de_casteljau_surf()
233 	    DCN(1, 0) = CN(1, 0, k) - CN(0, 0, k);  in _math_de_casteljau_surf()
234 	    DCN(0, 0) = us * CN(0, 0, k) + u * CN(1, 0, k);  in _math_de_casteljau_surf()
238 	       DCN(1, j + 1) = CN(1, j + 1, k) - CN(0, j + 1, k);  in _math_de_casteljau_surf()
242 	       DCN(0, j + 1) = us * CN(0, j + 1, k) + u * CN(1, j + 1, k);  in _math_de_casteljau_surf()
270 	    DCN(0, 1) = CN(0, 1, k) - CN(0, 0, k);  in _math_de_casteljau_surf()
271 	    DCN(0, 0) = vs * CN(0, 0, k) + v * CN(0, 1, k);  in _math_de_casteljau_surf()
274 	       DCN(i + 1, 1) = CN(i + 1, 1, k) - CN(i + 1, 0, k);  in _math_de_casteljau_surf()
278 	       DCN(i + 1, 0) = vs * CN(i + 1, 0, k) + v * CN(i + 1, 1, k);  in _math_de_casteljau_surf()
307 	    DCN(i, 0) = us * CN(i, 0, k) + u * CN(i + 1, 0, k);  in _math_de_casteljau_surf()
309 	       DCN(i, j + 1) = us * CN(i, j + 1, k) + u * CN(i + 1, j + 1, k);  in _math_de_casteljau_surf()
339 	    DCN(i, 0) = us * CN(i, 0, k) + u * CN(i + 1, 0, k);  in _math_de_casteljau_surf()
341 	       DCN(i, j + 1) = us * CN(i, j + 1, k) + u * CN(i + 1, j + 1, k);  in _math_de_casteljau_surf()
394 	    DCN(i, 0) = us * CN(i, 0, k) + u * CN(i + 1, 0, k);  in _math_de_casteljau_surf()
396 	       DCN(i, j + 1) = us * CN(i, j + 1, k) + u * CN(i + 1, j + 1, k);  in _math_de_casteljau_surf()
445 #undef CN  in _math_de_casteljau_surf()