def get_A(self):
A = NamedMatrix(square=True, rows=self.get_sectors())
if self.harmonized_rows is not None:
A.set_harmonized_rows(self.harmonized_rows)
stmt = self.generate_A_stmt()
result = stmt(*self.condition_args)
for row in result:
A.set_element(row[0], row[1], row[2]) # from_sec, to_sec, value
return A
