169    const vl_csc_matrix *cstd;  in vl_csc_get_matrix()  local
185          cstd = &bt_601;  in vl_csc_get_matrix()
188          cstd = &bt_709;  in vl_csc_get_matrix()
191          cstd = &smpte240m;  in vl_csc_get_matrix()
203    (*matrix)[0][0] = c * (*cstd)[0][0];  in vl_csc_get_matrix()
204    (*matrix)[0][1] = (*cstd)[0][1] * x - (*cstd)[0][2] * y;  in vl_csc_get_matrix()
205    (*matrix)[0][2] = (*cstd)[0][2] * x + (*cstd)[0][1] * y;  in vl_csc_get_matrix()
206    (*matrix)[0][3] = (*cstd)[0][3] + (*cstd)[0][0] * b +  in vl_csc_get_matrix()
207                      (*cstd)[0][1] * (x * cbbias + y * crbias) +  in vl_csc_get_matrix()
208                      (*cstd)[0][2] * (x * crbias - y * cbbias);  in vl_csc_get_matrix()
210    (*matrix)[1][0] = c * (*cstd)[1][0];  in vl_csc_get_matrix()
211    (*matrix)[1][1] = (*cstd)[1][1] * x - (*cstd)[1][2] * y;  in vl_csc_get_matrix()
212    (*matrix)[1][2] = (*cstd)[1][2] * x + (*cstd)[1][1] * y;  in vl_csc_get_matrix()
213    (*matrix)[1][3] = (*cstd)[1][3] + (*cstd)[1][0] * b +  in vl_csc_get_matrix()
214                      (*cstd)[1][1] * (x * cbbias + y * crbias) +  in vl_csc_get_matrix()
215                      (*cstd)[1][2] * (x * crbias - y * cbbias);  in vl_csc_get_matrix()
217    (*matrix)[2][0] = c * (*cstd)[2][0];  in vl_csc_get_matrix()
218    (*matrix)[2][1] = (*cstd)[2][1] * x - (*cstd)[2][2] * y;  in vl_csc_get_matrix()
219    (*matrix)[2][2] = (*cstd)[2][2] * x + (*cstd)[2][1] * y;  in vl_csc_get_matrix()
220    (*matrix)[2][3] = (*cstd)[2][3] + (*cstd)[2][0] * b +  in vl_csc_get_matrix()
221                      (*cstd)[2][1] * (x * cbbias + y * crbias) +  in vl_csc_get_matrix()
222                      (*cstd)[2][2] * (x * crbias - y * cbbias);  in vl_csc_get_matrix()