KaneAI - World’s First E2E Software Testing Agent.Book a Demo Join TestMu Conference

How to use Block method in root

Best JavaScript code snippet using root

lobpcg.py

Source:lobpcg.py

1"""2Locally Optimal Block Preconditioned Conjugate Gradient Method (LOBPCG).3References4----------5.. [1] A. V. Knyazev (2001),6       Toward the Optimal Preconditioned Eigensolver: Locally Optimal7       Block Preconditioned Conjugate Gradient Method.8       SIAM Journal on Scientific Computing 23, no. 2,9       pp. 517-541. http://dx.doi.org/10.1137/S106482750036612410.. [2] A. V. Knyazev, I. Lashuk, M. E. Argentati, and E. Ovchinnikov (2007),11       Block Locally Optimal Preconditioned Eigenvalue Xolvers (BLOPEX)12       in hypre and PETSc.  https://arxiv.org/abs/0705.262613.. [3] A. V. Knyazev's C and MATLAB implementations:14       https://github.com/lobpcg/blopex15"""16from __future__ import division, print_function, absolute_import17import numpy as np18from scipy.linalg import (inv, eigh, cho_factor, cho_solve, cholesky, orth,19                          LinAlgError)20from scipy.sparse.linalg import aslinearoperator21from scipy.sparse.sputils import bmat22__all__ = ['lobpcg']23def _report_nonhermitian(M, name):24    """25    Report if `M` is not a hermitian matrix given its type.26    """27    from scipy.linalg import norm28    md = M - M.T.conj()29    nmd = norm(md, 1)30    tol = 10 * np.finfo(M.dtype).eps31    tol = max(tol, tol * norm(M, 1))32    if nmd > tol:33        print('matrix %s of the type %s is not sufficiently Hermitian:'34              % (name, M.dtype))35        print('condition: %.e < %e' % (nmd, tol))36def _as2d(ar):37    """38    If the input array is 2D return it, if it is 1D, append a dimension,39    making it a column vector.40    """41    if ar.ndim == 2:42        return ar43    else:  # Assume 1!44        aux = np.array(ar, copy=False)45        aux.shape = (ar.shape[0], 1)46        return aux47def _makeOperator(operatorInput, expectedShape):48    """Takes a dense numpy array or a sparse matrix or49    a function and makes an operator performing matrix * blockvector50    products."""51    if operatorInput is None:52        return None53    else:54        operator = aslinearoperator(operatorInput)55    if operator.shape != expectedShape:56        raise ValueError('operator has invalid shape')57    return operator58def _applyConstraints(blockVectorV, factYBY, blockVectorBY, blockVectorY):59    """Changes blockVectorV in place."""60    YBV = np.dot(blockVectorBY.T.conj(), blockVectorV)61    tmp = cho_solve(factYBY, YBV)62    blockVectorV -= np.dot(blockVectorY, tmp)63def _b_orthonormalize(B, blockVectorV, blockVectorBV=None, retInvR=False):64    """B-orthonormalize the given block vector using Cholesky."""65    normalization = blockVectorV.max(axis=0)+np.finfo(blockVectorV.dtype).eps66    blockVectorV = blockVectorV / normalization67    if blockVectorBV is None:68        if B is not None:69            blockVectorBV = B(blockVectorV)70        else:71            blockVectorBV = blockVectorV  # Shared data!!!72    else:73        blockVectorBV = blockVectorBV / normalization74    VBV = np.matmul(blockVectorV.T.conj(), blockVectorBV)75    try:76        # VBV is a Cholesky factor from now on...77        VBV = cholesky(VBV, overwrite_a=True)78        VBV = inv(VBV, overwrite_a=True)79        blockVectorV = np.matmul(blockVectorV, VBV)80        # blockVectorV = (cho_solve((VBV.T, True), blockVectorV.T)).T81        if B is not None:82            blockVectorBV = np.matmul(blockVectorBV, VBV)83            # blockVectorBV = (cho_solve((VBV.T, True), blockVectorBV.T)).T84        else:85            blockVectorBV = None86    except LinAlgError:87        #raise ValueError('Cholesky has failed')88        blockVectorV = None89        blockVectorBV = None90        VBV = None91    if retInvR:92        return blockVectorV, blockVectorBV, VBV, normalization93    else:94        return blockVectorV, blockVectorBV95def _get_indx(_lambda, num, largest):96    """Get `num` indices into `_lambda` depending on `largest` option."""97    ii = np.argsort(_lambda)98    if largest:99        ii = ii[:-num-1:-1]100    else:101        ii = ii[:num]102    return ii103def lobpcg(A, X,104           B=None, M=None, Y=None,105           tol=None, maxiter=None,106           largest=True, verbosityLevel=0,107           retLambdaHistory=False, retResidualNormsHistory=False):108    """Locally Optimal Block Preconditioned Conjugate Gradient Method (LOBPCG)109    LOBPCG is a preconditioned eigensolver for large symmetric positive110    definite (SPD) generalized eigenproblems.111    Parameters112    ----------113    A : {sparse matrix, dense matrix, LinearOperator}114        The symmetric linear operator of the problem, usually a115        sparse matrix.  Often called the "stiffness matrix".116    X : ndarray, float32 or float64117        Initial approximation to the ``k`` eigenvectors (non-sparse). If `A`118        has ``shape=(n,n)`` then `X` should have shape ``shape=(n,k)``.119    B : {dense matrix, sparse matrix, LinearOperator}, optional120        The right hand side operator in a generalized eigenproblem.121        By default, ``B = Identity``.  Often called the "mass matrix".122    M : {dense matrix, sparse matrix, LinearOperator}, optional123        Preconditioner to `A`; by default ``M = Identity``.124        `M` should approximate the inverse of `A`.125    Y : ndarray, float32 or float64, optional126        n-by-sizeY matrix of constraints (non-sparse), sizeY < n127        The iterations will be performed in the B-orthogonal complement128        of the column-space of Y. Y must be full rank.129    tol : scalar, optional130        Solver tolerance (stopping criterion).131        The default is ``tol=n*sqrt(eps)``.132    maxiter : int, optional133        Maximum number of iterations.  The default is ``maxiter = 20``.134    largest : bool, optional135        When True, solve for the largest eigenvalues, otherwise the smallest.136    verbosityLevel : int, optional137        Controls solver output.  The default is ``verbosityLevel=0``.138    retLambdaHistory : bool, optional139        Whether to return eigenvalue history.  Default is False.140    retResidualNormsHistory : bool, optional141        Whether to return history of residual norms.  Default is False.142    Returns143    -------144    w : ndarray145        Array of ``k`` eigenvalues146    v : ndarray147        An array of ``k`` eigenvectors.  `v` has the same shape as `X`.148    lambdas : list of ndarray, optional149        The eigenvalue history, if `retLambdaHistory` is True.150    rnorms : list of ndarray, optional151        The history of residual norms, if `retResidualNormsHistory` is True.152    Notes153    -----154    If both ``retLambdaHistory`` and ``retResidualNormsHistory`` are True,155    the return tuple has the following format156    ``(lambda, V, lambda history, residual norms history)``.157    In the following ``n`` denotes the matrix size and ``m`` the number158    of required eigenvalues (smallest or largest).159    The LOBPCG code internally solves eigenproblems of the size ``3m`` on every160    iteration by calling the "standard" dense eigensolver, so if ``m`` is not161    small enough compared to ``n``, it does not make sense to call the LOBPCG162    code, but rather one should use the "standard" eigensolver, e.g. numpy or163    scipy function in this case.164    If one calls the LOBPCG algorithm for ``5m > n``, it will most likely break165    internally, so the code tries to call the standard function instead.166    It is not that ``n`` should be large for the LOBPCG to work, but rather the167    ratio ``n / m`` should be large. It you call LOBPCG with ``m=1``168    and ``n=10``, it works though ``n`` is small. The method is intended169    for extremely large ``n / m``, see e.g., reference [28] in170    https://arxiv.org/abs/0705.2626171    The convergence speed depends basically on two factors:172    1. How well relatively separated the seeking eigenvalues are from the rest173       of the eigenvalues. One can try to vary ``m`` to make this better.174    2. How well conditioned the problem is. This can be changed by using proper175       preconditioning. For example, a rod vibration test problem (under tests176       directory) is ill-conditioned for large ``n``, so convergence will be177       slow, unless efficient preconditioning is used. For this specific178       problem, a good simple preconditioner function would be a linear solve179       for `A`, which is easy to code since A is tridiagonal.180    References181    ----------182    .. [1] A. V. Knyazev (2001),183           Toward the Optimal Preconditioned Eigensolver: Locally Optimal184           Block Preconditioned Conjugate Gradient Method.185           SIAM Journal on Scientific Computing 23, no. 2,186           pp. 517-541. http://dx.doi.org/10.1137/S1064827500366124187    .. [2] A. V. Knyazev, I. Lashuk, M. E. Argentati, and E. Ovchinnikov188           (2007), Block Locally Optimal Preconditioned Eigenvalue Xolvers189           (BLOPEX) in hypre and PETSc. https://arxiv.org/abs/0705.2626190    .. [3] A. V. Knyazev's C and MATLAB implementations:191           https://bitbucket.org/joseroman/blopex192    Examples193    --------194    Solve ``A x = lambda x`` with constraints and preconditioning.195    >>> import numpy as np196    >>> from scipy.sparse import spdiags, issparse197    >>> from scipy.sparse.linalg import lobpcg, LinearOperator198    >>> n = 100199    >>> vals = np.arange(1, n + 1)200    >>> A = spdiags(vals, 0, n, n)201    >>> A.toarray()202    array([[  1.,   0.,   0., ...,   0.,   0.,   0.],203           [  0.,   2.,   0., ...,   0.,   0.,   0.],204           [  0.,   0.,   3., ...,   0.,   0.,   0.],205           ...,206           [  0.,   0.,   0., ...,  98.,   0.,   0.],207           [  0.,   0.,   0., ...,   0.,  99.,   0.],208           [  0.,   0.,   0., ...,   0.,   0., 100.]])209    Constraints:210    >>> Y = np.eye(n, 3)211    Initial guess for eigenvectors, should have linearly independent212    columns. Column dimension = number of requested eigenvalues.213    >>> X = np.random.rand(n, 3)214    Preconditioner in the inverse of A in this example:215    >>> invA = spdiags([1./vals], 0, n, n)216    The preconditiner must be defined by a function:217    >>> def precond( x ):218    ...     return invA @ x219    The argument x of the preconditioner function is a matrix inside `lobpcg`,220    thus the use of matrix-matrix product ``@``.221    The preconditioner function is passed to lobpcg as a `LinearOperator`:222    >>> M = LinearOperator(matvec=precond, matmat=precond,223    ...                    shape=(n, n), dtype=float)224    Let us now solve the eigenvalue problem for the matrix A:225    >>> eigenvalues, _ = lobpcg(A, X, Y=Y, M=M, largest=False)226    >>> eigenvalues227    array([4., 5., 6.])228    Note that the vectors passed in Y are the eigenvectors of the 3 smallest229    eigenvalues. The results returned are orthogonal to those.230    """231    blockVectorX = X232    blockVectorY = Y233    residualTolerance = tol234    if maxiter is None:235        maxiter = 20236    if blockVectorY is not None:237        sizeY = blockVectorY.shape[1]238    else:239        sizeY = 0240    # Block size.241    if len(blockVectorX.shape) != 2:242        raise ValueError('expected rank-2 array for argument X')243    n, sizeX = blockVectorX.shape244    if verbosityLevel:245        aux = "Solving "246        if B is None:247            aux += "standard"248        else:249            aux += "generalized"250        aux += " eigenvalue problem with"251        if M is None:252            aux += "out"253        aux += " preconditioning\n\n"254        aux += "matrix size %d\n" % n255        aux += "block size %d\n\n" % sizeX256        if blockVectorY is None:257            aux += "No constraints\n\n"258        else:259            if sizeY > 1:260                aux += "%d constraints\n\n" % sizeY261            else:262                aux += "%d constraint\n\n" % sizeY263        print(aux)264    A = _makeOperator(A, (n, n))265    B = _makeOperator(B, (n, n))266    M = _makeOperator(M, (n, n))267    if (n - sizeY) < (5 * sizeX):268        # warn('The problem size is small compared to the block size.' \269        #        ' Using dense eigensolver instead of LOBPCG.')270        sizeX = min(sizeX, n)271        if blockVectorY is not None:272            raise NotImplementedError('The dense eigensolver '273                                      'does not support constraints.')274        # Define the closed range of indices of eigenvalues to return.275        if largest:276            eigvals = (n - sizeX, n-1)277        else:278            eigvals = (0, sizeX-1)279        A_dense = A(np.eye(n, dtype=A.dtype))280        B_dense = None if B is None else B(np.eye(n, dtype=B.dtype))281        vals, vecs = eigh(A_dense, B_dense, eigvals=eigvals,282                          check_finite=False)283        if largest:284            # Reverse order to be compatible with eigs() in 'LM' mode.285            vals = vals[::-1]286            vecs = vecs[:, ::-1]287        return vals, vecs288    if (residualTolerance is None) or (residualTolerance <= 0.0):289        residualTolerance = np.sqrt(1e-15) * n290    # Apply constraints to X.291    if blockVectorY is not None:292        if B is not None:293            blockVectorBY = B(blockVectorY)294        else:295            blockVectorBY = blockVectorY296        # gramYBY is a dense array.297        gramYBY = np.dot(blockVectorY.T.conj(), blockVectorBY)298        try:299            # gramYBY is a Cholesky factor from now on...300            gramYBY = cho_factor(gramYBY)301        except LinAlgError:302            raise ValueError('cannot handle linearly dependent constraints')303        _applyConstraints(blockVectorX, gramYBY, blockVectorBY, blockVectorY)304    ##305    # B-orthonormalize X.306    blockVectorX, blockVectorBX = _b_orthonormalize(B, blockVectorX)307    ##308    # Compute the initial Ritz vectors: solve the eigenproblem.309    blockVectorAX = A(blockVectorX)310    gramXAX = np.dot(blockVectorX.T.conj(), blockVectorAX)311    _lambda, eigBlockVector = eigh(gramXAX, check_finite=False)312    ii = _get_indx(_lambda, sizeX, largest)313    _lambda = _lambda[ii]314    eigBlockVector = np.asarray(eigBlockVector[:, ii])315    blockVectorX = np.dot(blockVectorX, eigBlockVector)316    blockVectorAX = np.dot(blockVectorAX, eigBlockVector)317    if B is not None:318        blockVectorBX = np.dot(blockVectorBX, eigBlockVector)319    ##320    # Active index set.321    activeMask = np.ones((sizeX,), dtype=bool)322    lambdaHistory = [_lambda]323    residualNormsHistory = []324    previousBlockSize = sizeX325    ident = np.eye(sizeX, dtype=A.dtype)326    ident0 = np.eye(sizeX, dtype=A.dtype)327    ##328    # Main iteration loop.329    blockVectorP = None  # set during iteration330    blockVectorAP = None331    blockVectorBP = None332    iterationNumber = -1333    restart = True334    explicitGramFlag = False335    while iterationNumber < maxiter:336        iterationNumber += 1337        if verbosityLevel > 0:338            print('iteration %d' % iterationNumber)339        if B is not None:340            aux = blockVectorBX * _lambda[np.newaxis, :]341        else:342            aux = blockVectorX * _lambda[np.newaxis, :]343        blockVectorR = blockVectorAX - aux344        aux = np.sum(blockVectorR.conj() * blockVectorR, 0)345        residualNorms = np.sqrt(aux)346        residualNormsHistory.append(residualNorms)347        ii = np.where(residualNorms > residualTolerance, True, False)348        activeMask = activeMask & ii349        if verbosityLevel > 2:350            print(activeMask)351        currentBlockSize = activeMask.sum()352        if currentBlockSize != previousBlockSize:353            previousBlockSize = currentBlockSize354            ident = np.eye(currentBlockSize, dtype=A.dtype)355        if currentBlockSize == 0:356            break357        if verbosityLevel > 0:358            print('current block size:', currentBlockSize)359            print('eigenvalue:', _lambda)360            print('residual norms:', residualNorms)361        if verbosityLevel > 10:362            print(eigBlockVector)363        activeBlockVectorR = _as2d(blockVectorR[:, activeMask])364        if iterationNumber > 0:365            activeBlockVectorP = _as2d(blockVectorP[:, activeMask])366            activeBlockVectorAP = _as2d(blockVectorAP[:, activeMask])367            if B is not None:368                activeBlockVectorBP = _as2d(blockVectorBP[:, activeMask])369        if M is not None:370            # Apply preconditioner T to the active residuals.371            activeBlockVectorR = M(activeBlockVectorR)372        ##373        # Apply constraints to the preconditioned residuals.374        if blockVectorY is not None:375            _applyConstraints(activeBlockVectorR,376                              gramYBY, blockVectorBY, blockVectorY)377        ##378        # B-orthogonalize the preconditioned residuals to X.379        if B is not None:380            activeBlockVectorR = activeBlockVectorR - np.matmul(blockVectorX,381                                 np.matmul(blockVectorBX.T.conj(),382                                 activeBlockVectorR))383        else:384            activeBlockVectorR = activeBlockVectorR - np.matmul(blockVectorX,385                                 np.matmul(blockVectorX.T.conj(),386                                 activeBlockVectorR))387        ##388        # B-orthonormalize the preconditioned residuals.389        aux = _b_orthonormalize(B, activeBlockVectorR)390        activeBlockVectorR, activeBlockVectorBR = aux391        activeBlockVectorAR = A(activeBlockVectorR)392        if iterationNumber > 0:393            if B is not None:394                aux = _b_orthonormalize(B, activeBlockVectorP,395                                        activeBlockVectorBP, retInvR=True)396                activeBlockVectorP, activeBlockVectorBP, invR, normal = aux397            else:398                aux = _b_orthonormalize(B, activeBlockVectorP, retInvR=True)399                activeBlockVectorP, _, invR, normal = aux400            # Function _b_orthonormalize returns None if Cholesky fails401            if activeBlockVectorP is not None:402                activeBlockVectorAP = activeBlockVectorAP / normal403                activeBlockVectorAP = np.dot(activeBlockVectorAP, invR)404                restart = False405            else:406                restart = True407        ##408        # Perform the Rayleigh Ritz Procedure:409        # Compute symmetric Gram matrices:410        if activeBlockVectorAR.dtype == 'float32':411            myeps = 1412        elif activeBlockVectorR.dtype == 'float32':413            myeps = 1e-4414        else:415            myeps = 1e-8416        if residualNorms.max() > myeps and not explicitGramFlag:417            explicitGramFlag = False418        else:419            # Once explicitGramFlag, forever explicitGramFlag.420            explicitGramFlag = True421        # Shared memory assingments to simplify the code422        if B is None:423            blockVectorBX = blockVectorX424            activeBlockVectorBR = activeBlockVectorR425            if not restart:426                activeBlockVectorBP = activeBlockVectorP427        # Common submatrices:428        gramXAR = np.dot(blockVectorX.T.conj(), activeBlockVectorAR)429        gramRAR = np.dot(activeBlockVectorR.T.conj(), activeBlockVectorAR)430        if explicitGramFlag:431            gramRAR = (gramRAR + gramRAR.T.conj())/2432            gramXAX = np.dot(blockVectorX.T.conj(), blockVectorAX)433            gramXAX = (gramXAX + gramXAX.T.conj())/2434            gramXBX = np.dot(blockVectorX.T.conj(), blockVectorBX)435            gramRBR = np.dot(activeBlockVectorR.T.conj(), activeBlockVectorBR)436            gramXBR = np.dot(blockVectorX.T.conj(), activeBlockVectorBR)437        else:438            gramXAX = np.diag(_lambda)439            gramXBX = ident0440            gramRBR = ident441            gramXBR = np.zeros((sizeX, currentBlockSize), dtype=A.dtype)442        def _handle_gramA_gramB_verbosity(gramA, gramB):443            if verbosityLevel > 0:444                _report_nonhermitian(gramA, 'gramA')445                _report_nonhermitian(gramB, 'gramB')446            if verbosityLevel > 10:447                # Note: not documented, but leave it in here for now448                np.savetxt('gramA.txt', gramA)449                np.savetxt('gramB.txt', gramB)450        if not restart:451            gramXAP = np.dot(blockVectorX.T.conj(), activeBlockVectorAP)452            gramRAP = np.dot(activeBlockVectorR.T.conj(), activeBlockVectorAP)453            gramPAP = np.dot(activeBlockVectorP.T.conj(), activeBlockVectorAP)454            gramXBP = np.dot(blockVectorX.T.conj(), activeBlockVectorBP)455            gramRBP = np.dot(activeBlockVectorR.T.conj(), activeBlockVectorBP)456            if explicitGramFlag:457                gramPAP = (gramPAP + gramPAP.T.conj())/2458                gramPBP = np.dot(activeBlockVectorP.T.conj(),459                                 activeBlockVectorBP)460            else:461                gramPBP = ident462            gramA = bmat([[gramXAX, gramXAR, gramXAP],463                          [gramXAR.T.conj(), gramRAR, gramRAP],464                          [gramXAP.T.conj(), gramRAP.T.conj(), gramPAP]])465            gramB = bmat([[gramXBX, gramXBR, gramXBP],466                          [gramXBR.T.conj(), gramRBR, gramRBP],467                          [gramXBP.T.conj(), gramRBP.T.conj(), gramPBP]])468            _handle_gramA_gramB_verbosity(gramA, gramB)469            try:470                _lambda, eigBlockVector = eigh(gramA, gramB,471                                               check_finite=False)472            except LinAlgError:473                # try again after dropping the direction vectors P from RR474                restart = True475        if restart:476            gramA = bmat([[gramXAX, gramXAR],477                          [gramXAR.T.conj(), gramRAR]])478            gramB = bmat([[gramXBX, gramXBR],479                          [gramXBR.T.conj(), gramRBR]])480            _handle_gramA_gramB_verbosity(gramA, gramB)481            try:482                _lambda, eigBlockVector = eigh(gramA, gramB,483                                               check_finite=False)484            except LinAlgError:485                raise ValueError('eigh has failed in lobpcg iterations')486        ii = _get_indx(_lambda, sizeX, largest)487        if verbosityLevel > 10:488            print(ii)489            print(_lambda)490        _lambda = _lambda[ii]491        eigBlockVector = eigBlockVector[:, ii]492        lambdaHistory.append(_lambda)493        if verbosityLevel > 10:494            print('lambda:', _lambda)495#         # Normalize eigenvectors!496#         aux = np.sum( eigBlockVector.conj() * eigBlockVector, 0 )497#         eigVecNorms = np.sqrt( aux )498#         eigBlockVector = eigBlockVector / eigVecNorms[np.newaxis, :]499#         eigBlockVector, aux = _b_orthonormalize( B, eigBlockVector )500        if verbosityLevel > 10:501            print(eigBlockVector)502        # Compute Ritz vectors.503        if B is not None:504            if not restart:505                eigBlockVectorX = eigBlockVector[:sizeX]506                eigBlockVectorR = eigBlockVector[sizeX:sizeX+currentBlockSize]507                eigBlockVectorP = eigBlockVector[sizeX+currentBlockSize:]508                pp = np.dot(activeBlockVectorR, eigBlockVectorR)509                pp += np.dot(activeBlockVectorP, eigBlockVectorP)510                app = np.dot(activeBlockVectorAR, eigBlockVectorR)511                app += np.dot(activeBlockVectorAP, eigBlockVectorP)512                bpp = np.dot(activeBlockVectorBR, eigBlockVectorR)513                bpp += np.dot(activeBlockVectorBP, eigBlockVectorP)514            else:515                eigBlockVectorX = eigBlockVector[:sizeX]516                eigBlockVectorR = eigBlockVector[sizeX:]517                pp = np.dot(activeBlockVectorR, eigBlockVectorR)518                app = np.dot(activeBlockVectorAR, eigBlockVectorR)519                bpp = np.dot(activeBlockVectorBR, eigBlockVectorR)520            if verbosityLevel > 10:521                print(pp)522                print(app)523                print(bpp)524            blockVectorX = np.dot(blockVectorX, eigBlockVectorX) + pp525            blockVectorAX = np.dot(blockVectorAX, eigBlockVectorX) + app526            blockVectorBX = np.dot(blockVectorBX, eigBlockVectorX) + bpp527            blockVectorP, blockVectorAP, blockVectorBP = pp, app, bpp528        else:529            if not restart:530                eigBlockVectorX = eigBlockVector[:sizeX]531                eigBlockVectorR = eigBlockVector[sizeX:sizeX+currentBlockSize]532                eigBlockVectorP = eigBlockVector[sizeX+currentBlockSize:]533                pp = np.dot(activeBlockVectorR, eigBlockVectorR)534                pp += np.dot(activeBlockVectorP, eigBlockVectorP)535                app = np.dot(activeBlockVectorAR, eigBlockVectorR)536                app += np.dot(activeBlockVectorAP, eigBlockVectorP)537            else:538                eigBlockVectorX = eigBlockVector[:sizeX]539                eigBlockVectorR = eigBlockVector[sizeX:]540                pp = np.dot(activeBlockVectorR, eigBlockVectorR)541                app = np.dot(activeBlockVectorAR, eigBlockVectorR)542            if verbosityLevel > 10:543                print(pp)544                print(app)545            blockVectorX = np.dot(blockVectorX, eigBlockVectorX) + pp546            blockVectorAX = np.dot(blockVectorAX, eigBlockVectorX) + app547            blockVectorP, blockVectorAP = pp, app548    if B is not None:549        aux = blockVectorBX * _lambda[np.newaxis, :]550    else:551        aux = blockVectorX * _lambda[np.newaxis, :]552    blockVectorR = blockVectorAX - aux553    aux = np.sum(blockVectorR.conj() * blockVectorR, 0)554    residualNorms = np.sqrt(aux)555    # Future work: Need to add Postprocessing here:556    # Making sure eigenvectors "exactly" satisfy the blockVectorY constrains?557    # Making sure eigenvecotrs are "exactly" othonormalized by final "exact" RR558    # Computing the actual true residuals559    if verbosityLevel > 0:560        print('final eigenvalue:', _lambda)561        print('final residual norms:', residualNorms)562    if retLambdaHistory:563        if retResidualNormsHistory:564            return _lambda, blockVectorX, lambdaHistory, residualNormsHistory565        else:566            return _lambda, blockVectorX, lambdaHistory567    else:568        if retResidualNormsHistory:569            return _lambda, blockVectorX, residualNormsHistory570        else:...

js.js

Source:js.js

...45}46function robot(i,j){47    if(i<=startD[1]&&j<=endD[1]){48        if(block[i][j]==2){49            drawBlock(i+1,j,"red");50            drawBlock(i,j,"white");51            block[i+1][j]=3;52            block[i][j]=0;53        }else{54            drawBlock(i,j+1,"red");55            drawBlock(i,j,"white");56            block[i][j+1]=3;57            block[i][j]=0;58        }59    }else if(i>startD[1]&&j<=endD[1]){60        if(block[i][j]==2){61            drawBlock(i,j+1,"red");62            drawBlock(i,j,"white");63            block[i][j+1]=3;64            block[i][j]=0;65        }else{66            drawBlock(i-1,j,"red");67            drawBlock(i,j,"white");68            block[i-1][j]=3;69            block[i][j]=0;70        }71    }else if(i>startD[1]&&j>endD[1]){72        if(block[i][j]==2){73            drawBlock(i-1,j,"red");74            drawBlock(i,j,"white");75            block[i-1][j]=3;76            block[i][j]=0;77        }else{78            drawBlock(i,j-1,"red");79            drawBlock(i,j,"white");80            block[i][j-1]=3;81            block[i][j]=0;82        }83    }else if(i<=startD[1]&&j>endD[1]){84        if(block[i][j]==2){85            drawBlock(i,j-1,"red");86            drawBlock(i,j,"white");87            block[i][j-1]=3;88            block[i][j]=0;89        }else{90            drawBlock(i+1,j,"red");91            drawBlock(i,j,"white");92            block[i+1][j]=3;93            block[i][j]=0;94        }95    }96}97//ÒÆ¶¯»úÆ÷ÈË98//function mvRobot(){99//    for(var i=1;i<10;i++){100//        for(var j=1;j<10;j++){101//            if(block[i][j]==3){102//                if(i<=startD[1]){103//                    if(j<=endD[1]){104//                        if(block[i][j+1]==2){105//                            drawBlock(i+1,j,"red");106//                            drawBlock(i,j,"white");107//                            block[i+1][j]=3;108//                            block[i][j]=0;109//                        }else {110//                            if(j-endD[1]>=i-startD[1]){111//                                drawBlock(i,j+1,"red");112//                                drawBlock(i,j,"white");113//                                block[i][j+1]=3;114//                                block[i][j]=0;115//                            }else if(j-endD[1]<i-startD[1]){116//                                drawBlock(i+1,j,"red");117//                                drawBlock(i,j,"white");118//                                block[i+1][j]=3;119//                                block[i][j]=0;120//                            }121//                        }122//                    }else{123//                        if(block[i][j+1]==2){124//                            drawBlock(i-1,j,"red");125//                            drawBlock(i,j,"white");126//                            block[i-1][j]=3;127//                            block[i][j]=0;128//                        }else{129//                            if(j-endD[1]>=i-startD[1]){130//                                drawBlock(i,j-1,"red");131//                                drawBlock(i,j,"white");132//                                block[i][j-1]=3;133//                                block[i][j]=0;134//                            }else if(j-endD[1]<i-startD[1]){135//                                drawBlock(i-1,j,"red");136//                                drawBlock(i,j,"white");137//                                block[i-1][j]=3;138//                                block[i][j]=0;139//                            }140//                        }141//                    }142//                }else{143//                    if(j<=endD[1]){144//                        if(block[i][j+1]==2){145//                            drawBlock(i-1,j,"red");146//                            drawBlock(i,j,"white");147//                            block[i-1][j]=3;148//                            block[i][j]=0;149//                        }else{150//                            if(j-endD[1]>=i-startD[1]){151//                                drawBlock(i,j+1,"red");152//                                drawBlock(i,j,"white");153//                                block[i][j+1]=3;154//                                block[i][j]=0;155//                            }else if(j-endD[1]<i-startD[1]){156//                                drawBlock(i+1,j,"red");157//                                drawBlock(i,j,"white");158//                                block[i+1][j]=3;159//                                block[i][j]=0;160//                            }161//                        }162//                    }else{163//                        if(block[i][j+1]==2){164//                            drawBlock(i+1,j,"red");165//                            drawBlock(i,j,"white");166//                            block[i+1][j]=3;167//                            block[i][j]=0;168//                        }else{169//                            if(j-endD[1]>=i-startD[1]){170//                                drawBlock(i,j-1,"red");171//                                drawBlock(i,j,"white");172//                                block[i][j-1]=3;173//                                block[i][j]=0;174//                            }else if(j-endD[1]<i-startD[1]){175//                                drawBlock(i-1,j,"red");176//                                drawBlock(i,j,"white");177//                                block[i-1][j]=3;178//                                block[i][j]=0;179//                            }180//                        }181//                    }182//                }183//            }184//        }185//    }186//}187//ÒÆ¶¯hero188function mvSuper(pX,pY){189    if(block[startD[1]-pX][endD[1]-pY]==0){//-1 0190        startD[1] -= pX;191        endD[1]-=pY;192        if(superBlock==-1){193            drawBlock(startD[1],endD[1],"green");194        }else if(superBlock<3){195            drawBlock(startD[1],endD[1],"#6BF65E");196            superBlock+=1;197            if(superBlock==2){198                superBlock=-1;199            }200        }201        drawBlock(startD[1]+pX,endD[1]+pY,"white");//+1 0202        block[startD[1]][endD[1]]=-1;203        block[startD[1]+pX][endD[1]+pY]=0;//+1 0204    }else if(block[startD[1]-pX][endD[1]-pY]==2){//-1 0205        alert("no");206    }else if(block[startD[1]-pX][endD[1]-pY]==1){207        startD[1] -= pX;208        endD[1]-=pY;209        drawBlock(startD[1],endD[1],"#6BF65E");210        drawBlock(startD[1]+pX,endD[1]+pY,"white");211        block[startD[1]][endD[1]]=-1;212        block[startD[1]+pX][endD[1]+pY]=0;213        superBlock=0;214    }else if(block[startD[1]-pX][endD[1]-pY]==3){215        if(superBlock==-1) {216            alert("failed");217        } else if(superBlock!=-1){218            startD[1] -= pX;219            endD[1]-=pY;220            drawBlock(startD[1],endD[1],"#6BF65E");221            drawBlock(startD[1]+pX,endD[1]+pY,"white");222            block[startD[1]][endD[1]]=-1;223            block[startD[1]+1][endD[1]]=0;224            superBlock++;225        }226    }227    var robot=0;228    var suptool=0;229    for(var i=1;i<=10;i++){230        for(var j=1;j<=10;j++){231            if(block[i][j]==3){232                robot++;233            }234            if(block[i][j]==1){235                suptool++;236            }237        }238    }239    if(robot==0){240        alert("win");241    }else if(robot!=0 || (robot!=0&&suptool==0)){242        if(superBlock==-1) alert("failed");243        else{244            return;245        }246    }247}248//²úÉúhero robot super as....249function clickXY(){250    i=1;251    var x,y;252    chess.onclick=function(e){253        x= e.offsetX;254        y= e.offsetY;255        startD[i]= Math.floor(x / 30);256        endD [i]= Math.floor(y / 30);257        if(startD[i]<1 || startD[i]>10 || endD[i]<1 || endD[i]>10) return;258            if(i==1){259                drawBlock(startD[i],endD[i],"green");260                block[startD[i]][endD[i]]=-1;261                //console.log(startD[i],endD[i]);262                i++;263            }else if(i>1){264                if(block[startD[i]][endD[i]]==0){265                    var str=window.prompt("1,0,2","");266                    if(str==0){267                        drawBlock(startD[i],endD[i],"red");268                        block[startD[i]][endD[i]]=3;269                        i++;270                    }else if(str==1){271                        drawBlock(startD[i],endD[i],"blue");272                        block[startD[i]][endD[i]]=1;273                        i++;274                    }else if(str==2){275                        drawBlock(startD[i],endD[i],"#F6C462");276                        block[startD[i]][endD[i]]=2;277                        i++;278                    }279                }else alert("repeat");280            }281        else return;282    }283}284//»±í¸ñ285var drawchessbox=function(){286    for(var i=0;i<11;i++)287    {288        context.moveTo(30+i*30,30);289        context.lineTo(30+i*30,360-30);290        context.stroke();291        context.moveTo(30,30+i*30);292        context.lineTo(360-30,30+i*30);293        context.stroke();294    }295}296//»·½¿é297function drawBlock(startD,endD,color) {298    if(startD>0&&startD<11&&endD>0&&endD<11){299        context.strokeStyle = "block";300        context.fillStyle = color;301        context.fillRect(startD*30,endD*30,30,30);302        context.strokeRect(startD*30,endD*30,30,30);303    }304}305//if(block[startD[1]-1][endD[1]]==0){306//    startD[1] -= 1;307//    if(superBlock==-1){308//        drawBlock(startD[1],endD[1],"green");309//    }else if(superBlock<3){310//        drawBlock(startD[1],endD[1],"#6BF65E");311//        superBlock+=1;312//        if(superBlock==2){313//            superBlock=-1;314//        }315//    }316//    drawBlock(startD[1]+1,endD[1],"white");317//    block[startD[1]][endD[1]]=-1;318//    block[startD[1]+1][endD[1]]=0;319//}else if(block[startD[1]-1][endD[1]]==2){320//    alert("no");321//}else if(block[startD[1]-1][endD[1]]==1){322//    startD[1] -= 1;323//    drawBlock(startD[1],endD[1],"#6BF65E");324//    drawBlock(startD[1]+1,endD[1],"white");325//    block[startD[1]][endD[1]]=-1;326//    block[startD[1]+1][endD[1]]=0;327//    superBlock=0;328//}else if(block[startD[1]-1][endD[1]]==3){329//    if(superBlock==-1) {330//        alert("failed");331//    } else if(superBlock!=-1){332//        startD[1] -= 1;333//        drawBlock(startD[1],endD[1],"#6BF65E");334//        drawBlock(startD[1]+1,endD[1],"white");335//        block[startD[1]][endD[1]]=-1;336//        block[startD[1]+1][endD[1]]=0;337//        superBlock++;338//    }339//}340//startD[1] += 1;341//drawBlock(startD[1],endD[1],"green");342//drawBlock(startD[1]-1,endD[1],"white");343//endD[1] += 1;344//drawBlock(startD[1],endD[1],"green");345//drawBlock(startD[1],endD[1]-1,"white");346//endD[1] -= 1;347//drawBlock(startD[1],endD[1],"green");348//drawBlock(startD[1],endD[1]+1,"white");349//if(i<=startD[1]){350//    if(j<=endD[1]){351//        if(block[i][j+1]!=2){352//            drawBlock(i,j+1,"red");353//            drawBlock(i,j-1,"white");354//            block[i][j+1]=3;355//            block[i][j]=0;356//        }else{357//            drawBlock(i+1,j,"red");358//            drawBlock(i-1,j,"white");359//            block[i+1][j]=3;360//            block[i][j]=0;361//        }362//    }else{363//        if(block[i][j-1]!=2){364//            drawBlock(i,j-1,"red");365//            drawBlock(i,j+1,"white");366//            block[i][j-1]=3;367//            block[i][j]=0;368//        }else{369//            drawBlock(i-1,j,"red");370//            drawBlock(i+1,j,"white");371//            block[i-1][j]=3;372//            block[i][j]=0;373//        }374//    }375//}376//if(i>startD[1]){377//    if(j<=endD[1]){378//        if(block[i][j-1]!=2){379//            drawBlock(i,j-1,"red");380//            drawBlock(i,j+1,"white");381//            block[i][j-1]=3;382//            block[i][j]=0;383//        }else{384//            drawBlock(i-1,j,"red");385//            drawBlock(i+1,j,"white");386//            block[i-1][j]=3;387//            block[i][j]=0;388//        }389//    }else{390//        if(block[i][j+1]!=2){391//            drawBlock(i,j+1,"red");392//            drawBlock(i,j-1,"white");393//            block[i][j+1]=3;394//            block[i][j]=0;395//        }else{396//            drawBlock(i+1,j,"red");397//            drawBlock(i-1,j,"white");398//            block[i+1][j]=3;399//            block[i][j]=0;400//        }401//    }...

block-prop-container.ts

Source:block-prop-container.ts

1import {AfterViewInit, ChangeDetectorRef, Component, Inject, ViewChild} from "@angular/core";2import {Compbaser} from "ng-mslib";3import {YellowPepperService} from "../../services/yellowpepper.service";4import {BlockService, IBlockData, ISceneData} from "./block-service";5import {CampaignTimelineChanelPlayersModel} from "../../store/imsdb.interfaces_auto";6import {ColorPickerService} from "ngx-color-picker";7import {Tab} from "../../comps/tabs/tab";8import {BlockLabels, PLACEMENT_CHANNEL, PLACEMENT_SCENE} from "../../interfaces/Consts";9@Component({10    selector: 'block-prop-container',11    template: `12        <small class="debug">{{me}}</small>13        <tabs>14            <tab [tabtitle]="'style'">15                <block-prop-common [setBlockData]="m_blockData"></block-prop-common>16            </tab>17            <tab [tabtitle]="m_tabTitle">18                <div [ngSwitch]="m_blockTypeSelected">19                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_IMAGE">20                        <block-prop-image [external]="false" [setBlockData]="m_blockData"></block-prop-image>21                    </div>22                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_SVG">23                        <block-prop-image [external]="false" [setBlockData]="m_blockData"></block-prop-image>24                    </div>25                    <div *ngSwitchCase="m_blockLabels.IMAGE">26                        <block-prop-image [external]="true" [setBlockData]="m_blockData"></block-prop-image>27                    </div>28                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_VIDEO">29                        <block-prop-video [external]="false" [setBlockData]="m_blockData"></block-prop-video>30                    </div>31                    <div *ngSwitchCase="m_blockLabels.EXTERNAL_VIDEO">32                        <block-prop-video [external]="true" [setBlockData]="m_blockData"></block-prop-video>33                    </div>34                    <div *ngSwitchCase="m_blockLabels.MRSS">35                        <block-prop-mrss [setBlockData]="m_blockData"></block-prop-mrss>36                    </div>37                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_COLLECTION">38                        <block-prop-collection [setBlockData]="m_blockData"></block-prop-collection>39                    </div>40                    <div *ngSwitchCase="m_blockLabels.RSS">41                        <block-prop-rss [setBlockData]="m_blockData"></block-prop-rss>42                    </div>43                    <div *ngSwitchCase="m_blockLabels.LABEL">44                        <block-prop-label [setBlockData]="m_blockData"></block-prop-label>45                    </div>46                    <div *ngSwitchCase="m_blockLabels.QR">47                        <block-prop-qr [setBlockData]="m_blockData"></block-prop-qr>48                    </div>49                    <div *ngSwitchCase="m_blockLabels.CLOCK">50                        <block-prop-clock [setBlockData]="m_blockData"></block-prop-clock>51                    </div>52                    <div *ngSwitchCase="m_blockLabels.HTML">53                        <block-prop-html [setBlockData]="m_blockData"></block-prop-html>54                    </div>55                    <div *ngSwitchCase="m_blockLabels.FASTERQ">56                        <block-prop-fasterq [setBlockData]="m_blockData"></block-prop-fasterq>57                    </div>58                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_SCENE">59                        <block-prop-scene [setBlockData]="m_blockData"></block-prop-scene>60                    </div>61                    <div *ngSwitchCase="m_blockLabels.LOCATION">62                        <block-prop-location [setBlockData]="m_blockData"></block-prop-location>63                    </div>64                    <div *ngSwitchCase="m_blockLabels.YOUTUBE">65                        <block-prop-youtube [setBlockData]="m_blockData"></block-prop-youtube>66                    </div>67                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_JSON">68                        <block-prop-json-player [standAlone]="true" [setBlockData]="m_blockData"></block-prop-json-player>69                    </div>70                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_JSON_ITEM">71                        <block-prop-json-item [setBlockData]="m_blockData"></block-prop-json-item>72                    </div>73                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_WORLD_WEATHER">74                        <block-prop-weather [setBlockData]="m_blockData"></block-prop-weather>75                    </div>76                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_INSTAGRAM">77                        <block-prop-instagram [setBlockData]="m_blockData"></block-prop-instagram>78                    </div>79                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_DIGG">80                        <block-prop-digg [setBlockData]="m_blockData"></block-prop-digg>81                    </div>82                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_CALENDAR">83                        <block-prop-calendar [setBlockData]="m_blockData"></block-prop-calendar>84                    </div>85                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_TWITTERV3">86                        <block-prop-twitter [setBlockData]="m_blockData"></block-prop-twitter>87                    </div>88                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_GOOGLE_SHEETS">89                        <block-prop-sheets [setBlockData]="m_blockData"></block-prop-sheets>90                    </div>91                    <div *ngSwitchDefault>92                        <h3>no block prop found, new?</h3>93                        {{m_blockTypeSelected}}94                    </div>95                </div>96            </tab>97            <!-- below are JSON based blocked which bind to scenes -->98            <tab #settings [tabtitle]="'settings'">99                <div [ngSwitch]="m_blockTypeSelected">100                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_WORLD_WEATHER">101                        <block-prop-weather [jsonMode]="true" [setBlockData]="m_blockData"></block-prop-weather>102                    </div>103                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_INSTAGRAM">104                        <block-prop-instagram [jsonMode]="true" [setBlockData]="m_blockData"></block-prop-instagram>105                    </div>106                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_DIGG">107                        <block-prop-digg [jsonMode]="true" [setBlockData]="m_blockData"></block-prop-digg>108                    </div>109                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_CALENDAR">110                        <block-prop-calendar [jsonMode]="true" [setBlockData]="m_blockData"></block-prop-calendar>111                    </div>112                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_TWITTERV3">113                        <block-prop-twitter [jsonMode]="true" [setBlockData]="m_blockData"></block-prop-twitter>114                    </div>115                    <div *ngSwitchCase="m_blockLabels.BLOCKCODE_GOOGLE_SHEETS">116                        <block-prop-sheets [jsonMode]="true" [setBlockData]="m_blockData"></block-prop-sheets>117                    </div>118                </div>119            </tab>120        </tabs>121    `,122})123export class BlockPropContainer extends Compbaser implements AfterViewInit {124    m_blockTypeSelected: string = 'none';125    m_blockLabels = BlockLabels;126    m_blockData: IBlockData;127    m_tabTitle: string = 'none';128    m_showSettingsTab = true;129    constructor(@Inject('BLOCK_PLACEMENT') private blockPlacement: string, private yp: YellowPepperService, private bs: BlockService, private cpService: ColorPickerService, private cd: ChangeDetectorRef) {130        super();131        // console.log(blockPlacement);132        if (this.blockPlacement == PLACEMENT_CHANNEL)133            this._listenOnChannels();134        if (this.blockPlacement == PLACEMENT_SCENE)135            this._listenOnScenes();136    }137    @ViewChild('settings')138    settings: Tab;139    ngAfterViewInit() {140        this.toggleSettingsTab();141    }142    private _listenOnChannels() {143        this.cancelOnDestroy(144            //145            this.yp.listenBlockChannelSelectedOrChanged()146                .mergeMap((i_campaignTimelineChanelPlayersModel: CampaignTimelineChanelPlayersModel) => {147                    return this.bs.getBlockData(i_campaignTimelineChanelPlayersModel.getCampaignTimelineChanelPlayerId())148                })149                .subscribe((blockData: IBlockData) => {150                    this.m_blockTypeSelected = blockData.blockCode;151                    this.m_tabTitle = blockData.blockAcronym;152                    this.m_blockData = blockData;153                    // for json based scenes show the settings, unless its the actual Json Player which we don't154                    if (blockData.playerMimeScene && this.m_blockTypeSelected != '4300') {155                        this.m_showSettingsTab = true;156                        this.toggleSettingsTab();157                    } else {158                        this.m_showSettingsTab = false;159                        this.toggleSettingsTab();160                    }161                }, (e) => console.error(e))162        )163    }164    private _listenOnScenes() {165        this.cancelOnDestroy(166            //167            this.yp.listenSceneOrBlockSelectedChanged()168                .mergeMap((i_sceneData: ISceneData) => {169                    return this.bs.getBlockDataInScene(i_sceneData)170                })171                .subscribe((blockData: IBlockData) => {172                    this.m_blockTypeSelected = blockData.blockCode;173                    this.m_tabTitle = blockData.blockAcronym;174                    this.m_blockData = blockData;175                    this.m_showSettingsTab = false;176                    this.toggleSettingsTab();177                    this.cd.markForCheck();178                }, (e) => console.error(e))179        )180    }181    // private getTabTitle(i_blockData: IBlockData): string {182    //     // this.bs.getCommonBorderXML()if (i_blockData.blockAcronym.indexOf('JSON')) > -1 ? blockData.playerMimeScene : blockData.blockAcronym183    // }184    private toggleSettingsTab() {185        if (!this.settings) return;186        this.settings.show = this.m_showSettingsTab;187    }188    ngOnInit() {189    }190    destroy() {191    }...

zerocoin_publicSpend_reorg.py

Source:zerocoin_publicSpend_reorg.py

...37        self.unspent = self.node.listunspent()38        block_count = self.node.getblockcount()39        pastBlockHash = self.node.getblockhash(block_count)40        self.log.info("Block count: %d - Current best: %s..." % (self.node.getblockcount(), self.node.getbestblockhash()[:5]))41        pastBlock = CBlock()42        pastBlock.deserialize(BytesIO(hex_str_to_bytes(self.node.getblock(pastBlockHash, False))))43        checkpoint = pastBlock.nAccumulatorCheckpoint44        # 5) get the raw zerocoin spend txes45        self.log.info("Getting the raw zerocoin public spends...")46        public_spend_A = self.node.createrawzerocoinpublicspend(mints[0].get("serial hash"))47        tx_A = CTransaction()48        tx_A.deserialize(BytesIO(hex_str_to_bytes(public_spend_A)))49        tx_A.rehash()50        public_spend_B = self.node.createrawzerocoinpublicspend(mints[1].get("serial hash"))51        tx_B = CTransaction()52        tx_B.deserialize(BytesIO(hex_str_to_bytes(public_spend_B)))53        tx_B.rehash()54        # Spending same coins to different recipients to get different txids55        my_addy = "yAVWM5urwaTyhiuFQHP2aP47rdZsLUG5PH"...

block.js

Source:block.js

...16            return null17        }18        const web3 = await Web3Util.getWeb3()19        const start = new Date()20        const _block = await web3.eth.getBlock(blockNumber)21        const end = new Date() - start22        logger.info(`Execution time : %dms web3.eth.getBlock(${blockNumber})`, end)23        if (!_block) {24            return null25        }26        const { m1, m2 } = await utils.getM1M2(_block)27        _block.m2 = m228        _block.signer = m129        const timestamp = _block.timestamp * 100030        // Update end tx count.31        _block.timestamp = timestamp32        _block.e_tx = _block.transactions.length33        try {34            const data = {35                jsonrpc: '2.0',36                method: 'eth_getBlockFinalityByHash',37                params: [_block.hash],38                id: 8839            }40            const response = await axios.post(config.get('WEB3_URI'), data)41            const result = response.data42            let finalityNumber = parseInt(result.result)43            // blockNumber = 0 is genesis block44            if (parseInt(blockNumber) === 0) {45                finalityNumber = 10046            }47            _block.finality = finalityNumber48        } catch (e) {49            logger.warn('Cannot get block finality %s', blockNumber)50            logger.warn(e)51        }52        const txs = _block.transactions53        delete _block.transactions54        _block.status = true55        _block.updateFinalityTime = 056        delete _block._id57        delete _block.signers58        await db.Block.updateOne({ number: _block.number }, _block,59            { upsert: true, new: true })60        await elastic.index(_block.hash, 'blocks', _block)61        if (_block.number % config.get('BLOCK_PER_EPOCH') === 0) {62            const slashedNode = []63            const blk = await BlockHelper.getBlock(_block.number)64            if (blk.penalties && blk.penalties !== '0x') {65                const sbuff = Buffer.from((blk.penalties || '').substring(2), 'hex')66                if (sbuff.length > 0) {67                    for (let i = 1; i <= sbuff.length / 20; i++) {68                        const address = sbuff.slice((i - 1) * 20, i * 20)69                        const add = '0x' + address.toString('hex')70                        slashedNode.push(add.toLowerCase())71                    }72                }73            }74            const buff = Buffer.from(blk.extraData.substring(2), 'hex')75            const sbuff = buff.slice(32, buff.length - 65)76            const signers = []77            if (sbuff.length > 0) {78                for (let i = 1; i <= sbuff.length / 20; i++) {79                    const address = sbuff.slice((i - 1) * 20, i * 20)80                    signers.push('0x' + address.toString('hex'))81                }82            }83            const epoch = _block.number / config.get('BLOCK_PER_EPOCH')84            // TODO: update slash for next 5 epochs85            if (slashedNode.length > 0) {86                for (let i = 0; i < 5; i++) {87                    const nextEpoch = epoch + 1 + i88                    const e = await db.Epoch.findOne({ epoch: nextEpoch })89                    if (e) {90                        const sn = e.slashedNode91                        const newArr = sn.concat(slashedNode)92                        e.slashedNode = Array.from(new Set(newArr))93                        await e.save()94                    } else {95                        const ne = new db.Epoch({96                            epoch: nextEpoch,97                            startBlock: (nextEpoch - 1) * 900 + 1,98                            endBlock: nextEpoch * 900,99                            slashedNode: slashedNode,100                            isActive: false101                        })102                        await ne.save()103                    }104                }105            }106        }107        return { txs, timestamp }108    },109    getBlockDetail: async (hashOrNumber) => {110        try {111            let block = await db.Block.findOne({ number: hashOrNumber })112            if (!block) {113                block = await db.Block.findOne({ hash: String(hashOrNumber).toLowerCase() })114            }115            if (block && block.finality === 100) {116                return block117            }118            const web3 = await Web3Util.getWeb3()119            const _block = await web3.eth.getBlock(hashOrNumber)120            if (!_block) {121                return null122            }123            try {124                const { m1, m2 } = await utils.getM1M2(_block)125                _block.m2 = m2126                _block.signer = m1127            } catch (e) {128                logger.warn('Cannot get M1, M2 of block %s', hashOrNumber)129                logger.warn(e)130            }131            // Update end tx count.132            _block.timestamp = _block.timestamp * 1000133            _block.e_tx = _block.transactions.length134            try {135                const data = {136                    jsonrpc: '2.0',137                    method: 'eth_getBlockFinalityByHash',138                    params: [_block.hash],139                    id: 88140                }141                const response = await axios.post(config.get('WEB3_URI'), data)142                const result = response.data143                _block.finality = parseInt(result.result)144            } catch (e) {145                logger.warn('Cannot get block finality %s', hashOrNumber)146                logger.warn(e)147            }148            _block.status = true149            if (block) {150                if (!Object.prototype.hasOwnProperty.call(block, 'updateFinalityTime')) {151                    _block.updateFinalityTime = 0152                }153            } else {154                _block.updateFinalityTime = 0155            }156            delete _block._id157            delete _block.signers158            block = await db.Block.findOneAndUpdate({ number: _block.number }, _block,159                { upsert: true, new: true })160            return block161        } catch (e) {162            logger.warn('cannot get block %s with error %s', hashOrNumber, e)163            return {}164        }165    },166    getBlockOnChain: async (number) => {167        try {168            const _block = await BlockHelper.getBlock(number)169            if (!_block) {170                return null171            }172            // Get signer.173            const signer = await utils.toAddress(await utils.getSigner(_block), 100)174            _block.signer = signer.toLowerCase()175            // Update end tx count.176            _block.timestamp = _block.timestamp * 1000177            _block.e_tx = _block.transactions.length178            try {179                const data = {180                    jsonrpc: '2.0',181                    method: 'eth_getBlockFinalityByHash',182                    params: [_block.hash],183                    id: 88184                }185                const response = await axios.post(config.get('WEB3_URI'), data)186                const result = response.data187                let finalityNumber = parseInt(result.result)188                // blockNumber = 0 is genesis block189                if (parseInt(_block.number) === 0) {190                    finalityNumber = 100191                }192                _block.finality = finalityNumber193            } catch (e) {194                logger.warn('Cannot get block finality %s', number)195                logger.warn(e)196            }197            _block.status = true198            return _block199        } catch (e) {200            logger.warn('cannot get block on chain with error %s', e)201            return {}202        }203    },204    getBlock: async (number) => {205        const web3 = await Web3Util.getWeb3()206        return web3.eth.getBlock(number).catch(e => {207            logger.warn('Cannot get block %s detail. Sleep 2 seconds and try more. Error %s', number, e)208            return sleep(2000).then(() => {209                return BlockHelper.getBlock(number)210            })211        })212    },213    getLastBlockNumber: async () => {214        const web3 = await Web3Util.getWeb3()215        return web3.eth.getBlockNumber().catch(e => {216            logger.warn('Cannot get last block number. Sleep 2 seconds and try more. Error %s', e)217            return sleep(2000).then(() => {218                return BlockHelper.getLastBlockNumber()219            })220        })221    }222}223module.exports = BlockHelper

block.py

Source:block.py

...7        return hash(self) - hash(rhs)8    def __hash__(self):9        return (self.id << 8) + self.data10    def withData(self, data):11        return Block(self.id, data)12    def __iter__(self):13        """Allows a Block to be sent whenever id [and data] is needed"""14        return iter((self.id, self.data))15        16    def __repr__(self):17        return "Block(%d, %d)"%(self.id, self.data)18AIR                 = Block(0)19STONE               = Block(1)20GRASS               = Block(2)21DIRT                = Block(3)22COBBLESTONE         = Block(4)23WOOD_PLANKS         = Block(5)24SAPLING             = Block(6)25BEDROCK             = Block(7)26WATER_FLOWING       = Block(8)27WATER               = WATER_FLOWING28WATER_STATIONARY    = Block(9)29LAVA_FLOWING        = Block(10)30LAVA                = LAVA_FLOWING31LAVA_STATIONARY     = Block(11)32SAND                = Block(12)33GRAVEL              = Block(13)34GOLD_ORE            = Block(14)35IRON_ORE            = Block(15)36COAL_ORE            = Block(16)37WOOD                = Block(17)38LEAVES              = Block(18)39GLASS               = Block(20)40LAPIS_LAZULI_ORE    = Block(21)41LAPIS_LAZULI_BLOCK  = Block(22)42SANDSTONE           = Block(24)43BED                 = Block(26)44COBWEB              = Block(30)45GRASS_TALL          = Block(31)46WOOL                = Block(35)47FLOWER_YELLOW       = Block(37)48FLOWER_CYAN         = Block(38)49MUSHROOM_BROWN      = Block(39)50MUSHROOM_RED        = Block(40)51GOLD_BLOCK          = Block(41)52IRON_BLOCK          = Block(42)53STONE_SLAB_DOUBLE   = Block(43)54STONE_SLAB          = Block(44)55BRICK_BLOCK         = Block(45)56TNT                 = Block(46)57BOOKSHELF           = Block(47)58MOSS_STONE          = Block(48)59OBSIDIAN            = Block(49)60TORCH               = Block(50)61FIRE                = Block(51)62STAIRS_WOOD         = Block(53)63CHEST               = Block(54)64DIAMOND_ORE         = Block(56)65DIAMOND_BLOCK       = Block(57)66CRAFTING_TABLE      = Block(58)67FARMLAND            = Block(60)68FURNACE_INACTIVE    = Block(61)69FURNACE_ACTIVE      = Block(62)70DOOR_WOOD           = Block(64)71LADDER              = Block(65)72STAIRS_COBBLESTONE  = Block(67)73DOOR_IRON           = Block(71)74REDSTONE_ORE        = Block(73)75SNOW                = Block(78)76ICE                 = Block(79)77SNOW_BLOCK          = Block(80)78CACTUS              = Block(81)79CLAY                = Block(82)80SUGAR_CANE          = Block(83)81FENCE               = Block(85)82GLOWSTONE_BLOCK     = Block(89)83BEDROCK_INVISIBLE   = Block(95)84STONE_BRICK         = Block(98)85GLASS_PANE          = Block(102)86MELON               = Block(103)87FENCE_GATE          = Block(107)88GLOWING_OBSIDIAN    = Block(246)...

test_block.py

Source:test_block.py

...26    time.sleep(MINE_RATE / SECONDS)27    mined_block = Block.mine_block(last_block, 'bar')28    assert mined_block.difficulty == last_block.difficulty - 129def test_mined_block_difficulty_limits_at_1():30    last_block = Block(time.time_ns(), 'test_last_hash', 'test_hash',31                       'test_data', 1, 0)32    time.sleep(MINE_RATE / SECONDS)33    mined_block = Block.mine_block(last_block, 'bar')34    assert mined_block.difficulty == 135@pytest.fixture36def last_block():37    return Block.genesis_block()38@pytest.fixture39def block(last_block):40    return Block.mine_block(last_block, 'test_data')41def test_is_valid_block(last_block, block):42    Block.is_valid_block(last_block, block)43def test_is_valid_block_bad_last_hash(last_block, block):44    block.last_hash = 'evil_last_hash'...

index.ts

Source:index.ts

...15        this.data = data;16        this.timestamp = timestamp;17    }18}19const genesisBlock:Block = new Block(0, "2021", "", "Hello", 123456);20let blockchain: Block[] = [genesisBlock];21const getBlockchain = () : Block[] => blockchain;22const getLatestBlock = () : Block => blockchain[blockchain.length -1];23const getNewTimeStamp = () : number => Math.round(new Date().getTime() / 1000);24const createNewBlock = (data: string) : Block => {25    const previousBlock : Block = getLatestBlock();26    const newIndex : number = previousBlock.index + 1;27    const newTimestamp : number = getNewTimeStamp();28    const newHash : string = Block.calculateBlockHash(newIndex, previousBlock.hash, newTimestamp, data);29    const newBlock : Block = new Block(newIndex, newHash, previousBlock.hash, data, newTimestamp);30    addBlock(newBlock);31    return newBlock;32};33const getHashForBlock = (aBlock: Block) : string => Block.calculateBlockHash(aBlock.index, aBlock.previousHash, aBlock.timestamp, aBlock.data);34const isBlockValid = (candidateBlock : Block, previousBlock : Block): boolean => {35    if(!Block.validateStructure(candidateBlock)) {36        return false;37    } else if(previousBlock.index + 1 !== candidateBlock.index){38        return false;39    } else if(previousBlock.hash !== candidateBlock.previousHash){40        return false;41    } else if(getHashForBlock(candidateBlock) !== candidateBlock.hash){42        return false;43    } else {44        return true;45    }46};47const addBlock = (candidateBlock: Block) : void => {48    if(isBlockValid(candidateBlock, getLatestBlock())) {49        blockchain.push(candidateBlock);50    }51};52createNewBlock("second block");53createNewBlock("third block");54createNewBlock("fourth block");55console.log(blockchain);...

Using AI Code Generation

1var Block = require('./block');2var root = new Block('root');3root.add(new Block('block1'));4root.add(new Block('block2'));5root.add(new Block('block3'));6root.add(new Block('block4'));7root.add(new Block('block5'));8console.log(root);9var Block = require('./block');10var root = new Block('root');11root.add(new Block('block1'));12root.add(new Block('block2'));13root.add(new Block('block3'));14root.add(new Block('block4'));15root.add(new Block('block5'));16console.log(root);17console.log(root.childs[0]);18var Block = require('./block');19var root = new Block('root');20root.add(new Block('block1'));21root.add(new Block('block2'));22root.add(new Block('block3'));23root.add(new Block('block4'));24root.add(new Block('block5'));25console.log(root);26console.log(root.childs[0]);27console.log(root.childs[0].childs[0]);28var Block = require('./block');29var root = new Block('root');30root.add(new Block('block1'));31root.add(new Block('block2'));32root.add(new Block('block3'));33root.add(new Block('block4'));34root.add(new Block('block5'));35console.log(root);36console.log(root.childs[0]);37console.log(root.childs[0].childs[0]);38console.log(root.childs[0].childs[0].childs[0]);39var Block = require('./block');40var root = new Block('root');41root.add(new Block('block1'));42root.add(new Block('block2'));43root.add(new Block('block3'));44root.add(new Block('block4'));45root.add(new Block('block5'));46console.log(root);47console.log(root.childs[0]);48console.log(root.childs[0].childs[0]);49console.log(root.childs[0].childs[0].childs[0]);50console.log(root.childs[0].childs[0].childs[0].childs[0]);51var Block = require('./block');52var root = new Block('root');53root.add(new Block('block1'));54root.add(new Block('block2'));55root.add(new Block('block3

Using AI Code Generation

1root.Block(0, 0, 100, 100, 0.5, 0.5, 0.5, 0.5);2root.Block(0, 0, 100, 100, 0.5, 0.5, 0.5, 0.5);3root.Block(0, 0, 100, 100, 0.5, 0.5, 0.5, 0.5);4root.Block(0, 0, 100, 100, 0.5, 0.5, 0.5, 0.5);5root.Block(0, 0, 100, 100, 0.5, 0.5, 0.5, 0.5);6root.Block(0, 0, 100, 100, 0.5, 0.5, 0.5, 0.5);7root.Block(0, 0, 100, 100, 0.5, 0.5, 0.5, 0.5);8root.Block(0, 0, 100, 100, 0.5, 0.5, 0.5, 0.5);9root.Block(0, 0, 100, 100, 0.5, 0.5, 0.5, 0.5);10root.Block(0, 0, 100, 100, 0.5, 0.5, 0.5, 0.5);11root.Block(0, 0, 100, 100, 0.5, 0.5, 0.5, 0.5);

Using AI Code Generation

1var root = require('root');2var block = root.Block(function() {3  var x = 5;4  var y = 10;5  console.log(x + y);6});7block();8var root = require('root');9var block = root.Block(function() {10  var x = 5;11  var y = 10;12  console.log(x + y);13});14block();15var root = require('root');16var block = root.Block(function() {17  var x = 5;18  var y = 10;19  console.log(x + y);20});21block();22var root = require('root');23var block = root.Block(function() {24  var x = 5;25  var y = 10;26  console.log(x + y);27});28block();29var root = require('root');30var block = root.Block(function() {31  var x = 5;32  var y = 10;33  console.log(x + y);34});35block();36var root = require('root');37var block = root.Block(function() {38  var x = 5;39  var y = 10;40  console.log(x + y);41});42block();43var root = require('root');44var block = root.Block(function() {45  var x = 5;46  var y = 10;47  console.log(x + y);48});49block();50var root = require('root');51var block = root.Block(function() {52  var x = 5;53  var y = 10;54  console.log(x + y);55});56block();57var root = require('root');58var block = root.Block(function() {59  var x = 5;60  var y = 10;61  console.log(x + y);62});63block();64var root = require('root');

Automation Testing Tutorials

Learn to execute automation testing from scratch with LambdaTest Learning Hub. Right from setting up the prerequisites to run your first automation test, to following best practices and diving deeper into advanced test scenarios. LambdaTest Learning Hubs compile a list of step-by-step guides to help you be proficient with different test automation frameworks i.e. Selenium, Cypress, TestNG etc.