def rank(self, grid, gp, alphas, params, *args, **kws): gs = grid.getStorage() # estimate the convergence of ancestors of gp ratio = estimateConvergence(grid, gp, alphas) # get grid point associated to ix p = [gs.getCoordinates(gp, j) for j in range(gs.getDimension())] # get joint distribution ap = params.activeParams() U = ap.getIndependentJointDistribution() T = ap.getJointTransformation() q = T.unitToProbabilistic(p) # get area of basis function A = 1.0 for d in range(gs.getDimension()): A *= getIntegral(grid, gp.getLevel(d), gp.getIndex(d)) return abs(ratio) * U.pdf(q) * A
def rank(self, grid, gp, alphas, params, *args, **kws): gs = grid.getStorage() # estimate the convergence of ancestors of gp ratio = estimateConvergence(grid, gp, alphas) # get grid point associated to ix p = [gp.getCoord(j) for j in xrange(gs.dim())] # get joint distribution ap = params.activeParams() U = ap.getIndependentJointDistribution() T = ap.getJointTransformation() q = T.unitToProbabilistic(p) # get area of basis function A = 1.0 for d in xrange(gs.dim()): A *= getIntegral(grid, gp.getLevel(d), gp.getIndex(d)) return abs(ratio) * U.pdf(q) * A
def rank(self, grid, gp, alphas, *args, **kws): return np.abs(estimateConvergence(grid, gp, alphas))
def rank(self, grid, gp, alphas, *args, **kws): return abs(estimateConvergence(grid, gp, alphas))
for i in range(gs.getSize()): gs.getCoordinates(gs.getPoint(i), p) nodalValues[i] = f(p.array()) v = hierarchize(grid, nodalValues) ipar = 1 assert ipar < len(v) gp = gs.getPoint(ipar) print("Grand father:", gs.getCoordinate(gs.getPoint(0), 0), v[0]) print("Father:", gs.getCoordinate(gs.getPoint(0), 0), v[ipar]) gpl = HashGridPoint(gp) gpl.getLeftChild(0) gpr = HashGridPoint(gp) gpr.getRightChild(0) assert gsc.isContaining(gpl) assert gsc.isContaining(gpr) print(gs.getCoordinate(gpl, 0), \ estimateSurplus(grid, gpl, v), \ estimateConvergence(grid, gpl, v), \ w[gsc.getSequenceNumber(gpl)]) print(gs.getCoordinate(gpr, 0), \ estimateSurplus(grid, gpr, v), \ estimateConvergence(grid, gpr, v), \ w[gsc.getSequenceNumber(gpr)])