### ----------------------------------------------------
### -- Math::Libgsl::BLAS
### -- Licenses: Artistic-2.0
### -- Authors: Fernando Santagata
### -- File: t/01-raw.rakutest
### ----------------------------------------------------
### -- Chunk 3 of 4
.78e0, -.637e0, $c);
  mgsl_matrix_complex_float_set($cC1, 0, 1, $c);
  ok mgsl_blas_csymm(CblasLeft, CblasUpper, $cα, $cA1, $cB1, $cβ, $cC1) == GSL_SUCCESS, 'αAB + βC: complex32';
  is-deeply (gather for ^2 { mgsl_matrix_complex_float_get($cC1, 0, $_, $cres); take $cres.dat[0]; take $cres.dat[1] }), (1.4609999656677246e0, -0.1496639847755432e0, -0.8356919884681702e0, 0.36994391679763794e0), 'result: ok';

  mgsl_complex_rect(0e0, 0e0, $zα);
  mgsl_complex_rect(0e0, 0e0, $zβ);
  my gsl_matrix_complex $zA1 = gsl_matrix_complex_calloc(1, 1);
  mgsl_complex_rect(.511e0, -.486e0, $z);
  mgsl_matrix_complex_set($zA1, 0, 0, $z);
  my gsl_matrix_complex $zB1 = gsl_matrix_complex_calloc(1, 2);
  mgsl_complex_rect(.985e0, -.923e0, $z);
  mgsl_matrix_complex_set($zB1, 0, 0, $z);
  mgsl_complex_rect(-.234e0, -.756e0, $z);
  mgsl_matrix_complex_set($zB1, 0, 1, $z);
  my gsl_matrix_complex $zC1 = gsl_matrix_complex_calloc(1, 2);
  mgsl_complex_rect(-.16e0, .049e0, $z);
  mgsl_matrix_complex_set($zC1, 0, 0, $z);
  mgsl_complex_rect(.618e0, -.349e0, $z);
  mgsl_matrix_complex_set($zC1, 0, 1, $z);
  ok mgsl_blas_zsymm(CblasLeft, CblasUpper, $zα, $zA1, $zB1, $zβ, $zC1) == GSL_SUCCESS, 'αAB + βC: complex64';
  is-deeply (gather for ^2 { mgsl_matrix_complex_get($zC1, 0, $_, $zres); take $zres.dat[0]; take $zres.dat[1] }), (0e0, 0e0, 0e0, 0e0), 'result: ok';

  mgsl_complex_float_rect(0e0, .1e0, $cα);
  mgsl_complex_float_rect(0e0, .1e0, $cβ);
  my gsl_matrix_complex_float $cA2 = gsl_matrix_complex_float_calloc(1, 1);
  mgsl_complex_float_rect(-.126e0, .079e0, $c);
  mgsl_matrix_complex_float_set($cA2, 0, 0, $c);
  my gsl_matrix_complex_float $cB2 = gsl_matrix_complex_float_calloc(1, 2);
  mgsl_complex_float_rect(-.954e0, -.059e0, $c);
  mgsl_matrix_complex_float_set($cB2, 0, 0, $c);
  mgsl_complex_float_rect(.296e0, -.988e0, $c);
  mgsl_matrix_complex_float_set($cB2, 0, 1, $c);
  my gsl_matrix_complex_float $cC2 = gsl_matrix_complex_float_calloc(1, 2);
  mgsl_complex_float_rect(-.859e0, -.731e0, $c);
  mgsl_matrix_complex_float_set($cC2, 0, 0, $c);
  mgsl_complex_float_rect(.737e0, .593e0, $c);
  mgsl_matrix_complex_float_set($cC2, 0, 1, $c);
  ok mgsl_blas_chemm(CblasLeft, CblasUpper, $cα, $cA2, $cB2, $cβ, $cC2) == GSL_SUCCESS, 'sum αAB + βC: complex32';
  is-deeply (gather for ^2 { mgsl_matrix_complex_float_get($cC2, 0, $_, $cres); take $cres.dat[0]; take $cres.dat[1] }), (0.07235660403966904e0, -0.07387959957122803e0, -0.07174880057573318e0, 0.06997040659189224e0), 'result: ok';

  mgsl_complex_rect(0e0, .1e0, $zα);
  mgsl_complex_rect(0e0, .1e0, $zβ);
  my gsl_matrix_complex $zA2 = gsl_matrix_complex_calloc(1, 1);
  mgsl_complex_rect(-.359e0, .089e0, $z);
  mgsl_matrix_complex_set($zA2, 0, 0, $z);
  my gsl_matrix_complex $zB2 = gsl_matrix_complex_calloc(1, 2);
  mgsl_complex_rect(-.451e0, -.337e0, $z);
  mgsl_matrix_complex_set($zB2, 0, 0, $z);
  mgsl_complex_rect(-.901e0, -.871e0, $z);
  mgsl_matrix_complex_set($zB2, 0, 1, $z);
  my gsl_matrix_complex $zC2 = gsl_matrix_complex_calloc(1, 2);
  mgsl_complex_rect(.729e0, .631e0, $z);
  mgsl_matrix_complex_set($zC2, 0, 0, $z);
  mgsl_complex_rect(.364e0, .246e0, $z);
  mgsl_matrix_complex_set($zC2, 0, 1, $z);
  ok mgsl_blas_zhemm(CblasLeft, CblasUpper, $zα, $zA2, $zB2, $zβ, $zC2) == GSL_SUCCESS, 'sum αAB + βC: complex64';
  is-deeply (gather for ^2 { mgsl_matrix_complex_get($zC2, 0, $_, $zres); take $zres.dat[0]; take $zres.dat[1] }), (-0.0751983e0, 0.0890909e0, -0.0558689e0, 0.0687459e0), 'result: ok';

  $α = -.3e0;
  my gsl_matrix_float $A3 = gsl_matrix_float_calloc(2, 2);
  gsl_matrix_float_set($A3, 0, 0, .18e0);
  gsl_matrix_float_set($A3, 0, 1, .199e0);
  gsl_matrix_float_set($A3, 1, 0, .122e0);
  gsl_matrix_float_set($A3, 1, 1, -.547e0);
  my gsl_matrix_float $B3 = gsl_matrix_float_calloc(2, 3);
  gsl_matrix_float_set($B3, 0, 0, -.874e0);
  gsl_matrix_float_set($B3, 0, 1, -.383e0);
  gsl_matrix_float_set($B3, 0, 2, .458e0);
  gsl_matrix_float_set($B3, 1, 0, .124e0);
  gsl_matrix_float_set($B3, 1, 1, -.221e0);
  gsl_matrix_float_set($B3, 1, 2, -.107e0);
  ok gsl_blas_strmm(CblasLeft, CblasUpper, CblasNoTrans, CblasNonUnit, $α, $A3, $B3) == GSL_SUCCESS, 'αop(A)B: num32';
  is-deeply (gather for ^2 -> $i { for ^3 -> $j { take gsl_matrix_float_get($B3, $i, $j) } }), (0.03979320451617241e0, 0.03387570381164551e0, -0.018344102427363396e0, 0.020348399877548218e0, -0.03626609966158867e0, -0.01755870133638382e0), 'result: ok';

  $α = -.3e0;
  my gsl_matrix $dA3 = gsl_matrix_calloc(2, 2);
  gsl_matrix_set($dA3, 0, 0, .174e0);
  gsl_matrix_set($dA3, 0, 1, -.308e0);
  gsl_matrix_set($dA3, 1, 0, .997e0);
  gsl_matrix_set($dA3, 1, 1, -.484e0);
  my gsl_matrix $dB3 = gsl_matrix_calloc(2, 3);
  gsl_matrix_set($dB3, 0, 0, -.256e0);
  gsl_matrix_set($dB3, 0, 1, -.178e0);
  gsl_matrix_set($dB3, 0, 2, .098e0);
  gsl_matrix_set($dB3, 1, 0, .004e0);
  gsl_matrix_set($dB3, 1, 1, .97e0);
  gsl_matrix_set($dB3, 1, 2, -.408e0);
  ok gsl_blas_dtrmm(CblasLeft, CblasUpper, CblasNoTrans, CblasNonUnit, $α, $dA3, $dB3) == GSL_SUCCESS, 'αop(A)B: num64';
  is-deeply (gather for ^2 -> $i { for ^3 -> $j { take gsl_matrix_get($dB3, $i, $j) } }), (0.013732799999999998e0, 0.09891959999999998e0, -0.0428148e0, 0.0005808e0, 0.14084399999999997e0, -0.05924159999999999e0), 'result: ok';

  mgsl_complex_float_rect(-1e0, 0e0, $cα);
  my gsl_matrix_complex_float $cA3 = gsl_matrix_complex_float_calloc(2, 2);
  mgsl_complex_float_rect(-.023e0, .762e0, $c);
  mgsl_matrix_complex_float_set($cA3, 0, 0, $c);
  mgsl_complex_float_rect(-.687e0, -.039e0, $c);
  mgsl_matrix_complex_float_set($cA3, 0, 1, $c);
  mgsl_complex_float_rect(-.459e0, .047e0, $c);
  mgsl_matrix_complex_float_set($cA3, 1, 0, $c);
  mgsl_complex_float_rect(.189e0, .33e0, $c);
  mgsl_matrix_complex_float_set($cA3, 1, 1, $c);
  my gsl_matrix_complex_float $cB3 = gsl_matrix_complex_float_calloc(2, 3);
  mgsl_complex_float_rect(.827e0, -.561e0, $c);
  mgsl_matrix_complex_float_set($cB3, 0, 0, $c);
  mgsl_complex_float_rect(.641e0, -.229e0, $c);
  mgsl_matrix_complex_float_set($cB3, 0, 1, $c);
  mgsl_complex_float_rect(-.884e0, -.533e0, $c);
  mgsl_matrix_complex_float_set($cB3, 0, 2, $c);
  mgsl_complex_float_rect(-.624e0, -.138e0, $c);
  mgsl_matrix_complex_float_set($cB3, 1, 0, $c);
  mgsl_complex_float_rect(.073e0, .924e0, $c);
  mgsl_matrix_complex_float_set($cB3, 1, 1, $c);
  mgsl_complex_float_rect(-.501e0, -.164e0, $c);
  mgsl_matrix_complex_float_set($cB3, 1, 2, $c);
  ok mgsl_blas_ctrmm(CblasLeft, CblasUpper, CblasNoTrans, CblasNonUnit, $cα, $cA3, $cB3) == GSL_SUCCESS, 'αop(A)B: complex32';
  is-deeply (gather for ^2 -> $i { for ^3 -> $j { mgsl_matrix_complex_float_get($cB3, $i, $j, $cres); take $cres.dat[0]; take $cres.dat[1] } }), (-0.8317669630050659e0, -0.7622190117835999e0, -0.14564001560211182e0, 0.1439259946346283e0, -0.7642689943313599e0, 0.5291420221328735e0, 0.07239599525928497e0, 0.23200200498104095e0, 0.29112303256988525e0, -0.19872601330280304e0, 0.040568992495536804e0, 0.19632600247859955e0), 'result: ok';

  mgsl_complex_rect(0e0, 0e0, $zα);
  my gsl_matrix_complex $zA3 = gsl_matrix_complex_calloc(2, 2);
  mgsl_complex_rect(.463e0, .033e0, $z);
  mgsl_matrix_complex_set($zA3, 0, 0, $z);
  mgsl_complex_rect(-.929e0, .949e0, $z);
  mgsl_matrix_complex_set($zA3, 0, 1, $z);
  mgsl_complex_rect(.864e0, .986e0, $z);
  mgsl_matrix_complex_set($zA3, 1, 0, $z);
  mgsl_complex_rect(.393e0, .885e0, $z);
  mgsl_matrix_complex_set($zA3, 1, 1, $z);
  my gsl_matrix_complex $zB3 = gsl_matrix_complex_calloc(2, 3);
  mgsl_complex_rect(-.321e0, -.852e0, $z);
  mgsl_matrix_complex_set($zB3, 0, 0, $z);
  mgsl_complex_rect(-.337e0, -.175e0, $z);
  mgsl_matrix_complex_set($zB3, 0, 1, $z);
  mgsl_complex_rect(.607e0, -.613e0, $z);
  mgsl_matrix_complex_set($zB3, 0, 2, $z);
  mgsl_complex_rect(.688e0, .973e0, $z);
  mgsl_matrix_complex_set($zB3, 1, 0, $z);
  mgsl_complex_rect(-.331e0, -.35e0, $z);
  mgsl_matrix_complex_set($zB3, 1, 1, $z);
  mgsl_complex_rect(.719e0, -.553e0, $z);
  mgsl_matrix_complex_set($zB3, 1, 2, $z);
  ok mgsl_blas_ztrmm(CblasLeft, CblasUpper, CblasNoTrans, CblasNonUnit, $zα, $zA3, $zB3) == GSL_SUCCESS, 'αop(A)B: complex64';
  is-deeply (gather for ^2 -> $i { for ^3 -> $j { mgsl_matrix_complex_get($zB3, $i, $j, $zres); take $zres.dat[0]; take $zres.dat[1] } } ), (0e0, -0e0, 0e0, 0e0, 0e0, 0e0, -0e0, 0e0, 0e0, 0e0, 0e0, 0e0), 'result: ok';

  $α = -.3e0;
  my gsl_matrix_float $A4 = gsl_matrix_float_calloc(2, 2);
  gsl_matrix_float_set($A4, 0, 0, -.279e0);
  gsl_matrix_float_set($A4, 0, 1,  .058e0);
  gsl_matrix_float_set($A4, 1, 0,  .437e0);
  gsl_matrix_float_set($A4, 1, 1,  .462e0);
  my gsl_matrix_float $B4 = gsl_matrix_float_calloc(2, 3);
  gsl_matrix_float_set($B4, 0, 0,  .578e0);
  gsl_matrix_float_set($B4, 0, 1,  .473e0);
  gsl_matrix_float_set($B4, 0, 2,  -.34e0);
  gsl_matrix_float_set($B4, 1, 0, -.128e0);
  gsl_matrix_float_set($B4, 1, 1,  .503e0);
  gsl_matrix_float_set($B4, 1, 2,    .2e0);
  ok gsl_blas_strsm(CblasLeft, CblasUpper, CblasNoTrans, CblasNonUnit, $α, $A4, $B4) == GSL_SUCCESS, 'αop(inv(A))B: num32';
  is-deeply (gather for ^2 -> $i { for ^3 -> $j { take gsl_matrix_float_get($B4, $i, $j) } }), (0.6387841701507568e0, 0.4407019317150116e0, -0.3925895094871521e0, 0.08311688154935837e0, -0.3266233801841736e0, -0.12987013161182404e0), 'result: ok';

  $α = -1e0;
  my gsl_matrix $dA4 = gsl_matrix_calloc(2, 2);
  gsl_matrix_set($dA4, 0, 0, -.876e0);
  gsl_matrix_set($dA4, 0, 1, -.503e0);
  gsl_matrix_set($dA4, 1, 0, -.062e0);
  gsl_matrix_set($dA4, 1, 1, -.987e0);
  my gsl_matrix $dB4 = gsl_matrix_calloc(2, 3);
  gsl_matrix_set($dB4, 0, 0,  .219e0);
  gsl_matrix_set($dB4, 0, 1, -.986e0);
  gsl_matrix_set($dB4, 0, 2,    -0e0);
  gsl_matrix_set($dB4, 1, 0, -.605e0);
  gsl_matrix_set($dB4, 1, 1,  .289e0);
  gsl_matrix_set($dB4, 1, 2,  .641e0);
  ok gsl_blas_dtrsm(CblasLeft, CblasUpper, CblasNoTrans, CblasNonUnit, $α, $dA4, $dB4) == GSL_SUCCESS, 'αop(inv(A))B: num64';
  is-deeply (gather for ^2 -> $i { for ^3 -> $j { take gsl_matrix_get($dB4, $i, $j) } }), (0.6019671251382123e0, -1.2937005269415647e0, -0.3729106234935439e0, -0.6129685916919959e0, 0.292806484295846e0, 0.6494427558257345e0), 'result: ok';

  mgsl_complex_float_rect(0e0, 0e0, $cα);
  my gsl_matrix_complex_float $cA4 = gsl_matrix_complex_float_calloc(2, 2);
  mgsl_complex_float_rect(.491e0, -.317e0, $c);
  mgsl_matrix_complex_float_set($cA4, 0, 0, $c);
  mgsl_complex_float_rect(-.14e0, -.739e0, $c);
  mgsl_matrix_complex_float_set($cA4, 0, 1, $c);
  mgsl_complex_float_rect(-.969e0, -.518e0, $c);
  mgsl_matrix_complex_float_set($cA4, 1, 0, $c);
  mgsl_complex_float_rect(.702e0, -.287e0, $c);
  mgsl_matrix_complex_float_set($cA4, 1, 1, $c);
  my gsl_matrix_complex_float $cB4 = gsl_matrix_complex_float_calloc(2, 3);
  mgsl_complex_float_rect(-.962e0, -.38e0, $c);
  mgsl_matrix_complex_float_set($cB4, 0, 0, $c);
  mgsl_complex_float_rect(.656e0, .587e0, $c);
  mgsl_matrix_complex_float_set($cB4, 0, 1, $c);
  mgsl_complex_float_rect(-.195e0, -.862e0, $c);
  mgsl_matrix_complex_float_set($cB4, 0, 2, $c);
  mgsl_complex_float_rect(-.679e0, .598e0, $c);
  mgsl_matrix_complex_float_set($cB4, 1, 0, $c);
  mgsl_complex_float_rect(.919e0, .714e0, $c);
  mgsl_matrix_complex_float_set($cB4, 1, 1, $c);
  mgsl_complex_float_rect(-.513e0, .726e0, $c);
  mgsl_matrix_complex_float_set($cB4, 1, 2, $c);
  ok mgsl_blas_ctrsm(CblasLeft, CblasUpper, CblasNoTrans, CblasNonUnit, $cα, $cA4, $cB4) == GSL_SUCCESS, 'αop(inv(A))B: complex32';
  is-deeply (gather for ^2 -> $i { for ^3 -> $j { mgsl_matrix_complex_float_get($cB4, $i, $j, $cres); take $cres.dat[0]; take $cres.dat[1] } }), (0e0, 0e0, 0e0, 0e0, 0e0, 0e0, -0e0, 0e0, 0e0, 0e0, -0e0, 0e0), 'result: ok';

  mgsl_complex_rect(0e0, 0e0, $zα);
  my gsl_matrix_complex $zA4 = gsl_matrix_complex_calloc(2, 2);
  mgsl_complex_rect(.189e0, .519e0, $z);
  mgsl_matrix_complex_set($zA4, 0, 0, $z);
  mgsl_complex_rect(-.455e0, -.444e0, $z);
  mgsl_matrix_complex_set($zA4, 0, 1, $z);
  mgsl_complex_rect(-.21e0, -.507e0, $z);
  mgsl_matrix_complex_set($zA4, 1, 0, $z);
  mgsl_complex_rect(-.591e0, .859e0, $z);
  mgsl_matrix_complex_set($zA4, 1, 1, $z);
  my gsl_matrix_complex $zB4 = gsl_matrix_complex_calloc(2, 3);
  mgsl_complex_rect(-.779e0, -.484e0, $z);
  mgsl_matrix_complex_set($zB4, 0, 0, $z);
  mgsl_complex_rect(.249e0, -.107e0, $z);
  mgsl_matrix_complex_set($zB4, 0, 1, $z);
  mgsl_complex_rect(-.755e0, -.047e0, $z);
  mgsl_matrix_complex_set($zB4, 0, 2, $z);
  mgsl_complex_rect(.941e0, .675e0, $z);
  mgsl_matrix_complex_set($zB4, 1, 0, $z);
  mgsl_complex_rect(-.757e0, .645e0, $z);
  mgsl_matrix_complex_set($zB4, 1, 1, $z);
  mgsl_complex_rect(-.649e0, .242e0, $z);
  mgsl_matrix_complex_set($zB4, 1, 2, $z);
  ok mgsl_blas_ztrsm(CblasLeft, CblasUpper, CblasNoTrans, CblasNonUnit, $zα, $zA4, $zB4) == GSL_SUCCESS, 'αop(inv(A))B: complex64';
  is-deeply (gather for ^2 -> $i { for ^3 -> $j { mgsl_matrix_complex_get($zB4, $i, $j, $zres); take $zres.dat[0]; take $zres.dat[1] } } ), (0e0, -0e0, 0e0, 0e0, 0e0, 0e0, 0e0, -0e0, 0e0, 0e0, 0e0, 0e0), 'result: ok';

  $α = -1e0; $β = .1e0;
  my gsl_matrix_float $A5 = gsl_matrix_float_calloc(2, 1);
  gsl_matrix_float_set($A5, 0, 0,  .412e0);
  gsl_matrix_float_set($A5, 1, 0, -.229e0);
  my gsl_matrix_float $B5 = gsl_matrix_float_calloc(2, 2);
  gsl_matrix_float_set($B5, 0, 0,  .628e0);
  gsl_matrix_float_set($B5, 0, 1, -.664e0);
  gsl_matrix_float_set($B5, 1, 0, -.268e0);
  gsl_matrix_float_set($B5, 1, 1,  .096e0);
  ok gsl_blas_ssyrk(CblasUpper, CblasNoTrans, $α, $A5, $β, $B5) == GSL_SUCCESS, 'αAT(A) + βB: num32';
  is-deeply (gather for ^2 -> $i { for ^2 -> $j { take gsl_matrix_float_get($B5, $i, $j) } }), (-0.1069439947605133e0, 0.02794799953699112e0, -0.2680000066757202e0, -0.042841002345085144e0), 'result: ok';

  $α = .1e0; $β = -.3e0;
  my gsl_matrix $dA5 = gsl_matrix_calloc(2, 1);
  gsl_matrix_set($dA5, 0, 0,  .169e0);
  gsl_matrix_set($dA5, 1, 0, -.875e0);
  my gsl_matrix $dB5 = gsl_matrix_calloc(2, 2);
  gsl_matrix_set($dB5, 0, 0, .159e0);
  gsl_matrix_set($dB5, 0, 1, .277e0);
  gsl_matrix_set($dB5, 1, 0, .865e0);
  gsl_matrix_set($dB5, 1, 1, .346e0);
  ok gsl_blas_dsyrk(CblasUpper, CblasNoTrans, $α, $dA5, $β, $dB5) == GSL_SUCCESS, 'αAT(A) + βB: num64';
  is-deeply (gather for ^2 -> $i { for ^2 -> $j { take gsl_matrix_get($dB5, $i, $j) } }), (-0.0448439e0, -0.09788750000000002e0, 0.865e0, -0.027237499999999984e0), 'result: ok';

  mgsl_complex_float_rect(0e0, 0e0, $cα);
  mgsl_complex_float_rect(-.3e0, .1e0, $cβ);
  my gsl_matrix_complex_float $cA5 = gsl_matrix_complex_float_calloc(2, 1);
  mgsl_complex_float_rect(.447e0, -.507e0, $c);
  mgsl_matrix_complex_float_set($cA5, 0, 0, $c);
  mgsl_complex_float_rect(-.425e0, .701e0, $c);
  mgsl_matrix_complex_float_set($cA5, 1, 0, $c);
  my gsl_matrix_complex_float $cB5 = gsl_matrix_complex_float_calloc(2, 2);
  mgsl_complex_float_rect(.16e0, -.245e0, $c);
  mgsl_matrix_complex_float_set($cB5, 0, 0, $c);
  mgsl_complex_float_rect(.922e0, -.437e0, $c);
  mgsl_matrix_complex_float_set($cB5, 0, 1, $c);
  mgsl_complex_float_rect(.24e0, .008e0, $c);
  mgsl_matrix_complex_float_set($cB5, 1, 0, $c);
  mgsl_complex_float_rect(-.095e0, .749e0, $c);
  mgsl_matrix_complex_float_set($cB5, 1, 1, $c);
  ok mgsl_blas_csyrk(CblasUpper, CblasNoTrans, $cα, $cA5, $cβ, $cB5) == GSL_SUCCESS, 'αAT(A) + βB: complex32';
  is-deeply (gather for ^2 -> $i { for ^2 -> $j { mgsl_matrix_complex_float_get($cB5, $i, $j, $cres); take $cres.dat[0]; take $cres.dat[1] } }), (-0.023499999195337296e0, 0.08950001001358032e0, -0.2328999936580658e0, 0.2233000099658966e0, 0.23999999463558197e0, 0.00800000037997961e0, -0.046400003135204315e0, -0.23420001566410065e0), 'result: ok';

  mgsl_complex_rect(1e0, 0e0, $zα);
  mgsl_complex_rect(1e0, 0e0, $zβ);
  my gsl_matrix_complex $zA5 = gsl_matrix_complex_calloc(2, 1);
  mgsl_complex_rect(-.049e0, -.687e0, $z);
  mgsl_matrix_complex_set($zA5, 0, 0, $z);
  mgsl_complex_rect(-.434e0, .294e0, $z);
  mgsl_matrix_complex_set($zA5, 1, 0, $z);
  my gsl_matrix_complex $zB5 = gsl_matrix_complex_calloc(2, 2);
  mgsl_complex_rect(.937e0, -.113e0,