189    const vl_csc_matrix *cstd;  in vl_csc_get_matrix()  local
195          cstd = full_range ? &bt_601_full : &bt_601;  in vl_csc_get_matrix()
198          cstd = full_range ? &bt_709_full : &bt_709;  in vl_csc_get_matrix()
201          cstd = full_range ? &smpte240m_full : &smpte240m;  in vl_csc_get_matrix()
210    (*matrix)[0][0] = c * (*cstd)[0][0];  in vl_csc_get_matrix()
211    (*matrix)[0][1] = c * (*cstd)[0][1] * s * cosf(h) - c * (*cstd)[0][2] * s * sinf(h);  in vl_csc_get_matrix()
212    (*matrix)[0][2] = c * (*cstd)[0][2] * s * cosf(h) + c * (*cstd)[0][1] * s * sinf(h);  in vl_csc_get_matrix()
213    (*matrix)[0][3] = (*cstd)[0][3] + (*cstd)[0][0] * (b + c * ybias) +  in vl_csc_get_matrix()
214                      (*cstd)[0][1] * (c * cbbias * s * cosf(h) + c * crbias * s * sinf(h)) +  in vl_csc_get_matrix()
215                      (*cstd)[0][2] * (c * crbias * s * cosf(h) - c * cbbias * s * sinf(h));  in vl_csc_get_matrix()
217    (*matrix)[1][0] = c * (*cstd)[1][0];  in vl_csc_get_matrix()
218    (*matrix)[1][1] = c * (*cstd)[1][1] * s * cosf(h) - c * (*cstd)[1][2] * s * sinf(h);  in vl_csc_get_matrix()
219    (*matrix)[1][2] = c * (*cstd)[1][2] * s * cosf(h) + c * (*cstd)[1][1] * s * sinf(h);  in vl_csc_get_matrix()
220    (*matrix)[1][3] = (*cstd)[1][3] + (*cstd)[1][0] * (b + c * ybias) +  in vl_csc_get_matrix()
221                      (*cstd)[1][1] * (c * cbbias * s * cosf(h) + c * crbias * s * sinf(h)) +  in vl_csc_get_matrix()
222                      (*cstd)[1][2] * (c * crbias * s * cosf(h) - c * cbbias * s * sinf(h));  in vl_csc_get_matrix()
224    (*matrix)[2][0] = c * (*cstd)[2][0];  in vl_csc_get_matrix()
225    (*matrix)[2][1] = c * (*cstd)[2][1] * s * cosf(h) - c * (*cstd)[2][2] * s * sinf(h);  in vl_csc_get_matrix()
226    (*matrix)[2][2] = c * (*cstd)[2][2] * s * cosf(h) + c * (*cstd)[2][1] * s * sinf(h);  in vl_csc_get_matrix()
227    (*matrix)[2][3] = (*cstd)[2][3] + (*cstd)[2][0] * (b + c * ybias) +  in vl_csc_get_matrix()
228                      (*cstd)[2][1] * (c * cbbias * s * cosf(h) + c * crbias * s * sinf(h)) +  in vl_csc_get_matrix()
229                      (*cstd)[2][2] * (c * crbias * s * cosf(h) - c * cbbias * s * sinf(h));  in vl_csc_get_matrix()