Do not return CSL matrix

aimsgb/grain.py

@@ -220,9 +220,7 @@ def build_grains(self, csl, direction, uc_a=1, uc_b=1):
             uc_b (int): Number of unit cell of grain B. Default to 1.
 
         Returns:
-            Grain objects for grain A and B, and CSL matrix. The CSL matrix 
-            may be different from the input CSL matrix if the lattice system 
-            of the grain is not orthogonal.
+            Grain objects for grain A and B.
         """
         csl_t = csl.transpose()
         # rotate along a longer axis between a and b
@@ -247,6 +245,7 @@ def build_grains(self, csl, direction, uc_a=1, uc_b=1):
             csl = np.array([reduce_vector(i) for i in _csl])
             grain_a = self.copy()
             grain_a.make_supercell(csl)
+            print(csl)
 
             # grain_a.make_supercell(get_sc_fromstruct(grain_a, min_length=min(grain_a.lattice.abc),
                                                     #  force_diagonal=True))
@@ -274,7 +273,7 @@ def build_grains(self, csl, direction, uc_a=1, uc_b=1):
 
         # grain_b.to(filename='POSCAR_b_1')
         # exit()
-        return grain_a, grain_b, csl
+        return grain_a, grain_b
 
     @classmethod
     def stack_grains(cls, grain_a, grain_b, vacuum=0.0, gap=0.0, direction=2,