Skip to content

Kblock

from rgrow import State
from rgrow.kblock import KBlock, KBlockParams, KBlockTile

params = KBlockParams(
    [
        KBlockTile("S", 40e-9, ["null", "rb", "bb", "null"], "purple"),
        KBlockTile("RB", 40e-9, ["null", "rb", "e1", "rb*"], "red"),
        KBlockTile("BB", 100e-9, ["bb*", "e1", "bb", "null"], "blue"),
        KBlockTile("00", 100e-9, ["e0*", "e0", "e0", "e0*"], "white"),
        KBlockTile("10", 100e-9, ["e1*", "e1", "e1", "e0*"], "green"),
        KBlockTile("01", 100e-9, ["e0*", "e1", "e1", "e1*"], "yellow"),
        KBlockTile("11", 100e-9, ["e1*", "e0", "e0", "e1*"], "orange"),
    ],
    {"e1": 1000e-9, "e0": 1000e-9},
    seed={
        (4,4): "S"
    },
    binding_strength={
        "rb": "ACTGCTAGTACGTAGCGTAA",
        "bb": "ACTGCTAGTACGTAGCGTAA",
        "e0": "ACTGCTAGTA",
        "e1": "ACTGCTAGTA",
    },
    temp=47
)

sys = KBlock(params)
state = State((32, 32), "square", "none")
sys.setup_state(state)
sys.update_state(state)
sys.evolve(state, show_window=True)

EvolveOutcome.ReachedWallTimeMax

state.canvas_view[5,8]

np.uint32(0)

sys.plot_canvas(state)

<Axes: >

png

sys.temperature = 80

sys.evolve(state, show_window=True)

---------------------------------------------------------------------------

PanicException                            Traceback (most recent call last)

Cell In[16], line 1
----> 1 sys.evolve(state, show_window=True)


PanicException: Rate: PerSecond(1.3420742256956508), PointSafe2((4, 13)), SquareCanvasNullTracker(QuadTreeState { rates: QuadTreeSquareArray([[[PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), ..., PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), ..., PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), ..., PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), ..., PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), ..., PerSecond(0.9784088192093325), PerSecond(0.9784088192093325), PerSecond(0.9784101700403529), PerSecond(0.0), PerSecond(0.0)],
 ...,
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.9784088192093325), ..., PerSecond(0.19999999999999998), PerSecond(0.19999999999999998), PerSecond(0.19999999999999998), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(23.491802768354166), ..., PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.9784101700403529), ..., PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), ..., PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), ..., PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)]], shape=[32, 32], strides=[32, 1], layout=Cc (0x5), const ndim=2, [[PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.00016479046899986097), PerSecond(3.9077301606964268), PerSecond(24.670211587563497), PerSecond(2.3568176384186654), PerSecond(24.670211587563497), PerSecond(2.3568176384186654), PerSecond(2.3568176384186654), PerSecond(2.3568176384186654), PerSecond(2.3568176384186654), PerSecond(2.3568176384186654), PerSecond(2.34280530776177), PerSecond(2.349811473090218), PerSecond(2.356818989249686), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(19.989014970435466), PerSecond(0.0), PerSecond(0.9784088192093348), PerSecond(1.956817638418665), PerSecond(1.956817638418665), PerSecond(1.956817638418665), PerSecond(1.956817638418665), PerSecond(1.956817638418665), PerSecond(0.9784088192093348), PerSecond(0.9784088192093348), PerSecond(1.956817638418665), PerSecond(1.956817638418665), PerSecond(1.9568189892496854), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(2.3487992879183133), PerSecond(0.0), PerSecond(2.44280530776177), PerSecond(0.39999999999999997), PerSecond(0.39999999999999997), PerSecond(0.39999999999999997), PerSecond(0.39999999999999997), PerSecond(0.39198164949964864), PerSecond(2.4367839214358953), PerSecond(2.456818989249686), PerSecond(0.39999999999999997), PerSecond(0.39999999999999997), PerSecond(0.39999999999999997), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(2.2568189892496875), PerSecond(2.348800638749334), PerSecond(2.3498114730902175), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(24.662193237063146), PerSecond(2.356817638418665), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(2.356817638418665), PerSecond(0.0), PerSecond(24.670211587563497), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(2.340780937417962), PerSecond(2.3498114730902175), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(2.3498114730902175), PerSecond(0.0), PerSecond(2.3498114730902175), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(2.3381701882180175), PerSecond(2.34280530776177), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(18.5428710800907), PerSecond(0.0), PerSecond(2.34280530776177), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(2.340780937417962), PerSecond(24.670211587563497), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(24.6702115875635), PerSecond(0.0), PerSecond(2.3498114730902175), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(2.338170188218018), PerSecond(18.5428710800907), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(2.356817638418665), PerSecond(0.0), PerSecond(2.356817638418665), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(2.340780937417962), PerSecond(20.6221614947114), PerSecond(0.19999999999999998), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(2.356817638418665), PerSecond(0.0), PerSecond(2.3487992879183133), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(2.3330854098088096), PerSecond(2.448800638749334), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(2.356817638418665), PerSecond(0.0), PerSecond(2.2568189892496875), PerSecond(2.3487992879183137), PerSecond(0.000551165239588117), PerSecond(3.9077301606964268), PerSecond(2.3568176384186654), PerSecond(2.3568176384186654), PerSecond(2.3568176384186654), PerSecond(2.340780937417962), PerSecond(2.340780937417962), PerSecond(2.3337747720895146), PerSecond(18.53485335260785), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(2.356817638418665), PerSecond(0.0), PerSecond(18.5428710800907), PerSecond(0.0), PerSecond(3.9077301606964268), PerSecond(0.3), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(2.256817638418667), PerSecond(2.3568176384186654), PerSecond(2.3568176384186654), PerSecond(2.356818989249686), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(24.670212938394517), PerSecond(0.0), PerSecond(2.349812823921238), PerSecond(0.0), PerSecond(2.3568189892496854), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(24.670212938394517), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)]], shape=[16, 16], strides=[16, 1], layout=Cc (0x5), const ndim=2, [[PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(23.896909921600894), PerSecond(29.962255683610163), PerSecond(30.940664502819494), PerSecond(8.627270553674661), PerSecond(6.670452915256001), PerSecond(8.606252057689318), PerSecond(4.313637978499371)],
 [PerSecond(0.0), PerSecond(6.954418915917334), PerSecond(5.192616780851988), PerSecond(0.7999999999999999), PerSecond(0.7919816494996486), PerSecond(31.912613786167395), PerSecond(0.7999999999999999), PerSecond(0.39999999999999997)],
 [PerSecond(0.0), PerSecond(4.7066291115088825), PerSecond(27.020023060653713), PerSecond(0.0), PerSecond(0.0), PerSecond(9.371567906487966), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(43.2130826676542), PerSecond(4.692616780851988), PerSecond(0.0), PerSecond(0.0), PerSecond(47.892033793290175), PerSecond(0.0), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(4.71363527683733), PerSecond(4.705616926336978), PerSecond(0.0), PerSecond(0.0), PerSecond(27.744828480687502), PerSecond(0.19999999999999998), PerSecond(0.0)],
 [PerSecond(0.0), PerSecond(4.71363527683733), PerSecond(23.1484893572587), PerSecond(8.116011486632441), PerSecond(4.713635276837331), PerSecond(6.9544162142552945), PerSecond(9.388190986344807), PerSecond(20.891672341857536)],
 [PerSecond(0.0), PerSecond(24.670212938394517), PerSecond(2.349812823921238), PerSecond(2.3568189892496854), PerSecond(0.0), PerSecond(24.670212938394517), PerSecond(0.0), PerSecond(0.0)]], shape=[8, 8], strides=[8, 1], layout=Cc (0x5), const ndim=2, [[PerSecond(23.896909921600894), PerSecond(60.902920186429654), PerSecond(15.297723468930663), PerSecond(12.91989003618869)],
 [PerSecond(11.661048027426217), PerSecond(33.0126398415057), PerSecond(42.07616334215501), PerSecond(1.2)],
 [PerSecond(47.92671794449153), PerSecond(9.398233707188965), PerSecond(75.63686227397767), PerSecond(0.19999999999999998)],
 [PerSecond(29.383848215231847), PerSecond(35.97113265706207), PerSecond(36.33826442948714), PerSecond(30.279863328202342)]], shape=[4, 4], strides=[4, 1], layout=Cc (0x5), const ndim=2, [[PerSecond(129.47351797696246), PerSecond(71.49377684727436)],
 [PerSecond(122.67993252397441), PerSecond(142.45499003166717)]], shape=[2, 2], strides=[2, 1], layout=Cc (0x5), const ndim=2], PerSecond(466.1022173798784)), canvas: CanvasSquare([[0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 16, ..., 32, 32, 32, 0, 0],
 ...,
 [0, 0, 0, 0, 48, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 50, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 48, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0]], shape=[32, 32], strides=[32, 1], layout=Cc (0x5), const ndim=2), n_tiles: 174, total_events: 2749887, energy: 0.0, time: Second(15642.908290077261), tracker: NullStateTracker, tile_counts: [0], shape=[1], strides=[1], layout=CFcf (0xf), const ndim=1 }), [[0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 16, ..., 32, 32, 32, 0, 0],
 ...,
 [0, 0, 0, 0, 48, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 50, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 48, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0],
 [0, 0, 0, 0, 0, ..., 0, 0, 0, 0, 0]], shape=[32, 32], strides=[32, 1], layout=Cc (0x5), const ndim=2