code: def create(atoms_list,N, L, cutoff = 0, all_atomtypes=[]):
      
      myAlphas, myBetas = genBasis.getBasisFunc(cutoff, N)
      # get information about feature length
      n_datapoints = len(atoms_list)
      atoms = atoms_list[0]
      x = get_lastatom_soap(atoms, cutoff, myAlphas, myBetas,N,L, all_atomtypes=all_atomtypes)
      n_features = x.shape[1]
      print("soap first", x.shape)
      print(n_datapoints, n_features)
      soapmatrix = np.zeros((n_datapoints, n_features))
  
      i = -1
      for atoms in atoms_list:
          i +=1
          #atoms
          print("Processing " + str(atoms.info)," Run time: " + str(time.time()-t0_total), end="\r")
          soapmatrix[i,:] = get_lastatom_soap(atoms, cutoff, myAlphas, myBetas, N, L, all_atomtypes=all_atomtypes)
      print("")
  
      # infos
      print("shape", soapmatrix.shape)
      return soapmatrix
docs: Takes a trajectory xyz file and writes soap features