 smat_option = 'opt'
 decomp_option = '1d'
