code: def spatialDomainNoGrid(self):
      
      self.w = np.zeros(self.xw.shape)
  
      if self.Debug:
        print("w = ")
        print(self.w.shape)
      
      for i in range(len(self.q)):
        # More efficient if we have created some 0-load points
        # (e.g., for where we want output)
        if self.q[i] != 0:
          dist = np.abs(self.xw - self.x[i])
          self.w -= self.q[i] * self.coeff * np.exp(-dist/self.alpha) * \
            ( np.cos(dist/self.alpha) + np.sin(dist/self.alpha) )
docs: Superposition of analytical solutions without a gridded domain