docstr: Find the BD and FD finite differences along each component vector.
code: def surface_topo_ref(self):  # Validated
      """
      Find the BD and FD finite differences along each component vector.
      """
      # Replace numpy inf, -inf and nan objects with floating point numbers
      # nan --> float
      self.F[np.isnan(self.F)] = np.inf
      # inf, -inf  --> floats
      self.F = np.nan_to_num(self.F)
  
      self.Ft = self.F[self.Ind_sorted]
      self.Ftp = np.diff(self.Ft, axis=0)  # FD
      self.Ftm = np.diff(self.Ft[::-1], axis=0)[::-1]
  
imports: import numpy
  import time
  import logging
  import warnings
  from scipy import spatial
  from scipy.optimize import OptimizeResult, minimize
  from scipy.optimize._shgo_lib.triangulation import Complex
label_desc: Use image node_id 33 for calling the SHGO obj's underlying member method code with example usage: obj.surface_topo_ref() without return types