# in-plane orthotropy
ipo_now = ipo_param_1_12(result.lampam, mat_prop, constraints.sym)
table_result.loc[i, 'In-plane orthotropy parameter 1'] = ipo_now[0]
table_result.loc[i, 'In-plane orthotropy parameter 2'] = ipo_now[1]
table_result.loc[i, 'In-plane orthotropy parameter 3'] = ipo_now[2]
table_result.loc[i, 'In-plane orthotropy parameter 4'] = ipo_now[3]
table_result.loc[i, 'In-plane orthotropy parameter 5'] = ipo_now[4]
table_result.loc[i, 'In-plane orthotropy parameter 6'] = ipo_now[5]
table_result.loc[i, 'In-plane orthotropy parameter 7'] = ipo_now[6]
table_result.loc[i, 'In-plane orthotropy parameter 8'] = ipo_now[7]
table_result.loc[i, 'In-plane orthotropy parameter 9'] = ipo_now[8]
table_result.loc[i, 'In-plane orthotropy parameter 10'] = ipo_now[9]
table_result.loc[i, 'In-plane orthotropy parameter 11'] = ipo_now[10]
table_result.loc[i, 'In-plane orthotropy parameter 12'] = ipo_now[11]
A = A_from_lampam(result.lampam, mat_prop)
A11 = A[0, 0]
A22 = A[1, 1]
A12 = A[0, 1]
A66 = A[2, 2]
A16 = A[0, 2]
A26 = A[1, 2]
B = B_from_lampam(result.lampam, mat_prop)
B11 = B[0, 0]
B22 = B[1, 1]
B12 = B[0, 1]
B66 = B[2, 2]
B16 = B[0, 2]
B26 = B[1, 2]
Pop.loc[ipop, 'ply_counts'] = n_plies_in_panels
Pop.loc[ipop, 'lampam[1]'] = lampam[0]
Pop.loc[ipop, 'lampam[2]'] = lampam[1]
Pop.loc[ipop, 'lampam[3]'] = lampam[2]
Pop.loc[ipop, 'lampam[4]'] = lampam[3]
Pop.loc[ipop, 'lampam[5]'] = lampam[4]
Pop.loc[ipop, 'lampam[6]'] = lampam[5]
Pop.loc[ipop, 'lampam[7]'] = lampam[6]
Pop.loc[ipop, 'lampam[8]'] = lampam[7]
Pop.loc[ipop, 'lampam[9]'] = lampam[8]
Pop.loc[ipop, 'lampam[10]'] = lampam[9]
Pop.loc[ipop, 'lampam[11]'] = lampam[10]
Pop.loc[ipop, 'lampam[12]'] = lampam[11]
A = A_from_lampam(lampam, mat)
Pop.loc[ipop, 'A11'] = A[0, 0]
Pop.loc[ipop, 'A22'] = A[1, 1]
Pop.loc[ipop, 'A12'] = A[0, 1]
Pop.loc[ipop, 'A66'] = A[2, 2]
Pop.loc[ipop, 'A16'] = A[0, 2]
Pop.loc[ipop, 'A26'] = A[1, 2]
B = B_from_lampam(lampam, mat, constraints.sym)
Pop.loc[ipop, 'B11'] = B[0, 0]
Pop.loc[ipop, 'B22'] = B[1, 1]
Pop.loc[ipop, 'B12'] = B[0, 1]
Pop.loc[ipop, 'B66'] = B[2, 2]
Pop.loc[ipop, 'B16'] = B[0, 2]
lampam = calc_lampam(stack, constraints)
popu.loc[ind, 'lampam1'] = lampam[0]
popu.loc[ind, 'lampam2'] = lampam[1]
popu.loc[ind, 'lampam3'] = lampam[2]
popu.loc[ind, 'lampam4'] = lampam[3]
popu.loc[ind, 'lampam5'] = lampam[4]
popu.loc[ind, 'lampam6'] = lampam[5]
popu.loc[ind, 'lampam7'] = lampam[6]
popu.loc[ind, 'lampam8'] = lampam[7]
popu.loc[ind, 'lampam9'] = lampam[8]
popu.loc[ind, 'lampam10'] = lampam[9]
popu.loc[ind, 'lampam11'] = lampam[10]
popu.loc[ind, 'lampam12'] = lampam[11]
A = A_from_lampam(lampam, mat_prop)
B = B_from_lampam(lampam, mat_prop, constraints.sym)
D = D_from_lampam(lampam, mat_prop)
# filter_ABD(A=A, B=B, D=D, sym=constraints.sym, ply_t=mat_prop.ply_t)
popu.loc[ind, 'A11'] = A[0, 0]
popu.loc[ind, 'A22'] = A[1, 1]
popu.loc[ind, 'A12'] = A[0, 1]
popu.loc[ind, 'A66'] = A[2, 2]
popu.loc[ind, 'A16'] = A[0, 2]
popu.loc[ind, 'A26'] = A[1, 2]
popu.loc[ind, 'B11'] = B[0, 0]
popu.loc[ind, 'B22'] = B[1, 1]
popu.loc[ind, 'B12'] = B[0, 1]