def METHOD_NAME(R,s,kup,kdn):
  M=np.zeros((len(R),len(kup)),dtype=complex);
  for iat in range(0,len(R)):
    for norb in range(0,len(kup)):
      M[iat][norb]=spinor_val(R[iat],s[iat],kup[norb],kdn[norb])
  return M