def METHOD_NAME(self):
    G1 = Gf(mesh=self.iw_mesh, target_shape=(1,1))
    G2 = G1.copy()
    G3 = Gf(mesh=self.iw_mesh, target_shape=(1,2))
    G4 = Gf(mesh=self.iw_mesh, target_shape=(1,3))
    G1 << inverse(2 + iOmega_n)
    G2 << inverse(2 - iOmega_n)
    for i in range(2):
        G3[0,i] << inverse(iOmega_n + 2 + i)
    for i in range(3):
        G4[0,i] << inverse(iOmega_n + 2 + i)
    B1 = BlockGf(name_list=['0', '1'], block_list=[G1, G2])
    B2 = BlockGf(name_list=['0', '1'], block_list=[G3, G4])
    B3 = B1 * B2
    B4 = BlockGf(name_list=['0', '1'], block_list=[G1*G3, G2*G4])
    assert_block_gfs_are_close(B3, B4)