Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

Extended-variable probabilistic computing with p-dits

T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Pith's one-line read Probabilistic p-dits beat binary p-bits by up to 34x and software solvers by 64x

desk verdict Real hardware and a useful generalization, but the p-int update's stationary distribution is unproven and may be sampling at 2β, which undercuts the integer-programming claims until fixed. read the letter →

arxiv 2506.00269 v1 pith:ZFCHOBWW submitted 2025-05-30 physics.app-ph

classification physics.app-ph
keywords probabilisticcomputingp-ditp-bitp-intIsingmachinecombinatorialoptimizationcategoricalvariablesintegerlinearprogramming
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper claims that the binary p-bit, the standard unit of probabilistic Ising machines, can be usefully generalized to a p-dit that holds more than two states, and that this generalization translates directly into faster probabilistic optimization on real hardware. The authors define a probabilistic d-dimensional bit and then specialize it to two practical forms: isotropic p-dits, which encode categorical variables with no invalid states, and probabilistic integers, which encode numeric variables directly on a number line. On a 130 nm CMOS ASIC, the isotropic p-dit3 solves a 14-number 3-partition problem in about 34x fewer trials than the p-bit implementation, and the p-int solves a change-making integer linear program in about 5.3x fewer trials. An FPGA p-int solver on a non-convex integer quadratic program reaches the true solution about 64x faster than the best of several state-of-the-art software solvers. If the formulation scales, p-dits give Ising machines a way around the resource and energy-barrier costs of binary and one-hot encodings.

What carries the argument

The central object is the probabilistic d-dimensional bit (p-dit), a spin that stochastically oscillates among discrete states in $|D|$-dimensional space. Its update rule is a multi-input sigmoid over the energy differences to all candidate states, so lower-energy states are selected with higher probability. Two restricted forms carry the hardware results: an isotropic p-dit, whose allowed states are the orthogonal coordinate axes and whose couplings reduce to a single $N \times N$ matrix, and a p-int, whose allowed states are consecutive integers along one axis and whose updates are limited to stepping to a neighboring value. These restrictions make extended probabilistic variables compact enough to fabricate, and the isotropic p-dit is explicitly related to the planar/clock Potts model.

What would settle it

Measure the long-run state histogram of a single p-int on the ASIC with a fixed local field $I$, no self-coupling, and no other spins, and compare it with the Boltzmann weights $\exp(\beta I m)$ over its full range; a systematic mismatch would show that the restricted update does not sample the distribution implied by the temperature parameter.

Watch

Extended reading notes

Core claim

The central claim is that probabilistic computing does not have to be built from two-state p-bits: one p-dit can represent a d-dimensional spin, and two restricted classes of p-dit remove the main encoding bottlenecks that arise when Ising machines are applied to categorical and integer optimization. An isotropic p-dit restricts the spin to lie along exactly one of the orthogonal axes in its state space, so a categorical variable is always in a valid state; this eliminates the exponentially growing set of invalid one-hot assignments and the need to tune a constraint-to-objective constant. A p-int restricts the spin to integer levels along one axis and updates by increments or decrements of one, which avoids the energy barriers that binary representations create between adjacent integer values. The authors derive Glauber-like update probabilities for both cases, argue in supplementary notes that the resulting Markov chains converge, and demonstrate the performance gains experimentally on an ASIC and an FPGA.

Load-bearing premise

The load-bearing premise is that a p-int restricted to stepping only to its current value and its two immediate neighbors converges to a stationary distribution matching the intended energy landscape, with the proof deferred to Supplementary Note 5.

Editorial extensions

If this is right

  • For problems with many categorical variables, the fraction of invalid one-hot states grows exponentially, so the advantage of isotropic p-dits should widen as problem size grows.
  • One p-int replaces roughly $\lceil \log_2(B-A) \rceil$ p-bits for an integer variable with bounds $A$ and $B$, reducing hardware footprint and the number of energy barriers a solver must cross.
  • Inequality constraints encoded with asymmetric violation variables avoid the deep energy valleys of slack variables, giving an additional ~10x improvement in trials-to-solution on a fixed-charge ILP.
  • A p-int probabilistic solver can outperform conventional software MILP/MIQCP solvers in time-to-solution on a non-convex IQP problem, suggesting a practical role for this hardware approach.
  • A single reconfigurable p-element can act as a p-bit, a p-int, or an isotropic p-dit3, so one CMOS chip can host several problem types with a two-clock-cycle update.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • We infer that the reported 34x and 5.3x ratios are likely to be lower bounds for larger instances, since both sources of overhead that p-dits remove grow with problem size and word length.
  • The violation-variable idea is not tied to p-ints: the same asymmetric constraint representation could be adapted to isotropic p-dits or other extended variables, which the paper does not explore.
  • Because isotropic p-dits coincide with the clock Potts model, annealing schedules and theoretical results from Potts-model physics could be transferred directly to tune p-dit-based solvers.
  • Graph coloring and scheduling are natural next benchmarks for isotropic p-dits, since they share the labeling-symmetry condition that makes the isotropic coupling matrix valid.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper generalizes probabilistic Ising machines from binary p-bits to p-dits, probabilistic units with d-dimensional states, and analyzes two restrictions: isotropic p-dits for categorical variables and p-ints for integer variables. The authors derive a full-state Glauber-like update that converges to the Boltzmann distribution and a restricted update for p-ints that is claimed to converge. They fabricate a 130-nm ASIC with 64 configurable p-elements and implement an FPGA-based p-int solver. Measured results include ~34x fewer trials-to-solution for a 3-partition problem using isotropic p-dit3 versus p-bits, ~5.3x for a change-making ILP using p-ints versus p-bits, and ~64x faster time-to-solution for a non-convex IQP compared with several GAMS solvers. The paper also introduces violation variables and scaled sampling as practical heuristics for inequality constraints and variable ranges.

Significance. The empirical hardware results are substantial: a fabricated ASIC implementing all three variable types, an elegant isotropic p-dit formulation that avoids invalid one-hot states, and large measured speedups on representative optimization problems. The full-state and isotropic p-dit theory is sound as far as presented, and the performance numbers are measured rather than fitted. The main significance is the demonstration that extended variables can reduce hardware resource usage and improve time-to-solution for problems with categorical or integer structure. However, the theoretical status of the restricted p-int update needs clarification before the 'probabilistic computing' framing is fully justified; the p-int results currently stand as strong heuristics rather than exact Boltzmann sampling.

major comments (3)
  1. [Probabilistic Computing with P-ints, Eq. (30)] The restricted p-int update does not yield the Boltzmann stationary distribution of the stated Hamiltonian. For a single p-int with constant local field I and J_ii=0, the transition probabilities in Eq. (30) give P_up/P_down = exp(2βI), so the birth-death chain has π_m ∝ exp(2βI m), not exp(-βE) ∝ exp(βI m). Thus p-int annealing operates at an effective inverse temperature 2β for this case, and in general the stationary distribution is not specified. The main text only states that the update 'will result in a converging probability distribution' (deferred to Supplementary Note 5). Because the paper's framing presents p-ints as probabilistic computing elements, the authors should either provide the actual stationary distribution and its relation to the Hamiltonian, or explicitly state that p-ints are a non-Boltzmann heuristic and adjust the theoretical narrative accordingly.
  2. [Violation Variables] The violation-variable construction produces an asymmetric J matrix. The Boltzmann result in Eq. (10) and the energy-difference formulas are derived under the symmetric-coupling assumption stated before Eq. (8). The paper does not specify what energy function, if any, the asymmetric update kernel is sampling, nor does it prove convergence for this case. Since the fixed-charge and IQP benchmarks rely on violation variables, the theoretical basis for these results should be stated explicitly, or the scheme should be presented as a heuristic with appropriate caveats.
  3. [Results/Integer Programming, Fig. 6] The reported 5.3x improvement in trials-to-solution is an ASIC p-int result compared with a simulated p-bit baseline ('simulated results for p-bits, and both simulated and experimental results for p-ints'). Because the ASIC supports p-bits, the comparison should either use the ASIC p-bit baseline or be explicitly labeled as hardware-versus-simulation. Similarly, the IQP comparison in Fig. 8 is based on a single problem instance, 25 hardware trials, and one trial per software solver; the generality of the ~64x claim would be strengthened by multiple instances and seeds, and these limitations should be stated.
minor comments (4)
  1. [Throughout] Typos and leftover fragments: 'dimentional' in Fig. 1, 'quadradic' in the Discussion, '3-partion' in Table 1, and the orphaned 'or53.' in §Isotropic P-dits should be corrected.
  2. [References] References [32] and [47] cite the same paper (Camsari et al., Stochastic p-bits for invertible logic) and should be consolidated.
  3. [Scaled Sampling] The scaled-sampling procedure is described only qualitatively; the formula for the sampling ratio as a function of p-int ranges should be given in the main text or Methods.
  4. [Probabilistic Computing with P-ints, Eq. (28)] The statement in Eq. (28) is consistent with Eq. (24) when I_i is evaluated at the current state; the self-coupling terms cancel in the two-step difference. The text should clarify this to avoid apparent inconsistency.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the p-dit theory and benchmarks are self-contained and measured; the main weakness is an unsupported stationary-distribution claim for the restricted p-int update, which is a correctness issue rather than a circular reduction.

full rationale

The paper's derivation chain is largely self-contained. The generalized spin Hamiltonian (Eq. 5) and the Glauber-like update rules (Eqs. 9 and 11) define the stochastic process; the claim that the full-update chain converges to the Boltzmann distribution (Eq. 10) is a standard detailed-balance result, not an input fitted from data. The isotropic p-dit and p-int restrictions are specializations of this framework, and the reported ~34x, ~5.3x, ~10x, and ~64x improvements are measured trials-to-solution or time-to-solution numbers from an ASIC and an FPGA, not quantities derived from the model, so there is no fitted-input-called-prediction. Hyperparameters (beta, C, O, cutoffs) are tuned per benchmark, which is standard practice and not circular. There are some self-citations (e.g., refs. 3, 9, 39), but they are background or roadmap references and are not load-bearing for any central claim; no uniqueness theorem is imported from the authors' prior work, and no ansatz is smuggled in via citation. The one substantive weakness is that the restricted p-int update of Eq. (30) is asserted to converge, with proof deferred to Supplementary Note 5, and the two-step energy difference in Eq. (28) appears to omit the -2J_ii term required by Eq. (24). This is a missing or invalid mathematical justification, i.e., a correctness risk, not a circular reduction: the update rule is not defined in terms of the claimed stationary distribution, and the benchmark speedups do not rest on that distribution being exactly Boltzmann. Hence the circularity score is 1.

Assumptions & free parameters 5 free parameters · 5 assumptions · 4 invented entities

The generalized spin model is a straightforward vector-spin Hamiltonian with standard Boltzmann or Glauber dynamics for the full-state update. No parameters are fitted in the theory itself. The empirical claims rest on hand-tuned hyperparameters (C, O, β, cutoffs) and on two heuristics: the restricted one-step update for p-ints and the violation-variable encoding for inequality constraints. The convergence of the restricted update to a Boltzmann distribution is asserted rather than derived and appears false in the single-spin case.

free parameters (5)
  • Inverse temperature β = e.g., 1/64 (p-int change-making), 1/128 (p-bit change-making), β0=1/256 and 1/4 for fixed-charge, β0=3.3 for IQP
    Tuned per problem and implementation to maximize success rate; not derived.
  • Constraint scaling C = e.g., 94 (3-partition), 1 (change-making, fixed-charge)
    Chosen to balance constraint enforcement and objective optimization; swept in experiments.
  • Objective scaling O = e.g., 1 (3-partition), 96 (change-making), 2 (slack fixed-charge), 1/4 (violation fixed-charge), 0.0008 (IQP)
    Tuned per problem and solver configuration.
  • IQP reset cutoff = 2^21 and 2^28 iterations
    Selected to optimize time-to-solution based on the distribution of solution times.
  • Sampling ratio for scaled sampling = proportional to p-int range
    Ad hoc method to equalize exploration across variables with different ranges.
assumptions (5)
  • standard math Boltzmann distribution and Glauber dynamics for multi-state spins
    Used in Eq. (10) and the derivation of update probabilities.
  • domain assumption Quadratic Hamiltonian captures the relevant energy landscape
    The generalized model (Eq. 5) considers only up to second-order interactions; higher-order terms are mentioned but not used.
  • ad hoc to paper Restricted state-dependent updates converge to the Boltzmann stationary distribution
    Assumed for p-ints at Eq. (30) and deferred to Supplementary Note 5; the main text equations suggest the stationary distribution is not Boltzmann.
  • ad hoc to paper Asymmetric J matrix in violation variables preserves the optimization dynamics
    Violation variables break symmetry of J to improve exploration, with no proof that the resulting dynamics converge to any desired distribution.
  • ad hoc to paper Scaled sampling improves exploration
    Proportional sampling of p-ints by range is introduced heuristically.
invented entities (4)
  • p-dit (probabilistic d-dimensional bit) independent evidence
    purpose: Base unit of the generalized spin model; oscillates among discrete states in d-dimensional space.
    Implemented as 64 multi-purpose p-elements on a fabricated ASIC and measured on partition and ILP problems.
  • Isotropic p-dit independent evidence
    purpose: Represent categorical variables with orthogonal states, avoiding invalid one-hot encodings.
    Implemented on ASIC; 3-partition and 6-partition results show solution quality improvements.
  • p-int (probabilistic integer) independent evidence
    purpose: Represent integer variables with regular steps along one dimension.
    Implemented on ASIC and FPGA; shows speedups on ILP and IQP problems.
  • Violation variables independent evidence
    purpose: Replace slack variables for inequality constraints in ILP and IQP Hamiltonians.
    Tested on fixed-charge ILP; shows about 10x improvement over slack variables with even sampling.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Extended-variable probabilistic computing with p-dits." pith.science (2026). https://pith.science/paper/ZFCHOBWW

@misc{pith2026250600269,
  author       = {Pith},
  title        = {Pith review of: Extended-variable probabilistic computing with p-dits},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZFCHOBWW}},
  note         = {Machine review of arXiv:2506.00269}
}
read the original abstract

Ising machines can solve combinatorial optimization problems by representing them as energy minimization problems. A common implementation is the probabilistic Ising machine (PIM), which uses probabilistic (p-) bits to represent coupled binary spins. However, many real-world problems have complex data representations that do not map naturally into a binary encoding, leading to a significant increase in hardware resources and time-to-solution. Here, we describe a generalized spin model that supports an arbitrary number of spin dimensions, each with an arbitrary real component. We define the probabilistic d-dimensional bit (p-dit) as the base unit of a p-computing implementation of this model. We further describe two restricted forms of p-dits for specific classes of common problems and implement them experimentally on an application-specific integrated circuit (ASIC): (A) isotropic p-dits, which simplify the implementation of categorical variables resulting in ~34x performance improvement compared to a p-bit implementation on an example 3-partition problem. (B) Probabilistic integers (p-ints), which simplify the representation of numeric values and provide ~5x improvement compared to a p-bit implementation of an example integer linear programming (ILP) problem. Additionally, we report a field-programmable gate array (FPGA) p-int-based integer quadratic programming (IQP) solver which shows ~64x faster time-to-solution compared to the best of a series of state-of-the-art software solvers. The generalized formulation of probabilistic variables presented here provides a path to solving large-scale optimization problems on various hardware platforms including digital CMOS.

Discussion (0). Sign in to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Pathology-Informed Latent Diffusion Model for Anomaly Detection in Lymph Node Metastasis

    eess.IV 2025-08 unverdicted novelty 5.0 of 10

    A pathology prompt-guided latent diffusion model reconstructs normal tissue and flags abnormal patches by reconstruction error, targeting unsupervised lymph node metastasis detection.

Reference graph

Works this paper leans on

60 extracted references · 59 canonical work pages · cited by 1 Pith paper

  1. [1]

    Aramon, M. et al. Physics-inspired optimization for quadratic unconstrained problems using a digital annealer. Frontiers in Physics 7, 48 (2019)

  2. [2]

    & Lin, M

    Zou, Y. & Lin, M. Massively simulating adiabatic bifurcations with fpga to solve combinatorial optimization. In Proceedings of the 2020 ACM/SIGDA International Symposium on Field Programmable Gate Arrays, 65–75 (2020)

  3. [3]

    Finocchio, G. et al. Roadmap for unconventional computing with nanotechnology. Nano Futures (2023)

  4. [4]

    Simulating physics with computers

    Feynman, R. Simulating physics with computers. International Journal of Theoretical Physics 21 (1982)

  5. [5]

    & Huang, M

    Sharma, A., Burns, M., Hahn, A. & Huang, M. Augmenting an electronic ising machine to effectively solve boolean satisfiability. Scientific Reports 13, 22858 (2023)

  6. [6]

    & Han, J

    Zhang, T. & Han, J. Efficient traveling salesman problem solvers using the ising model with simulated bifurcation. In 2022 Design, Automation & Test in Europe Conference & Exhibition (DATE), 548–551 (IEEE, 2022)

  7. [7]

    Si, J. et al. Energy-efficient superparamagnetic ising machine and its application to traveling salesman problems. Nature Communications 15, 3457 (2024)

  8. [8]

    Lu, A. et al. Scalable in-memory clustered annealer with temporal noise of charge trap transistor for large scale travelling salesman problems. IEEE Journal on Emerging and Selected Topics in Circuits and Systems 13, 422–435 (2023)

Show all 60 references
  1. [9]

    Finocchio, G. et al. The promise of spintronics for unconventional computing. Journal of Magnetism and Magnetic Materials 521, 167506 (2021)

  2. [10]

    Mohseni, N., McMahon, P. L. & Byrnes, T. Ising machines as hardware solvers of combinatorial optimization problems. Nature Reviews Physics 4, 363–379 (2022)

  3. [11]

    & Tanaka, S

    Tanahashi, K., Takayanagi, S., Motohashi, T. & Tanaka, S. Application of ising machines and a software development for ising machines. Journal of the Physical Society of Japan 88, 061010 (2019)

  4. [13]

    W., Rosenbluth, M

    Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H., & Teller, E. Equation of state calculations by fast computing machines. The Journal of Chemical Physics 21(6), 1087-1092 (1953)

  5. [14]

    Litvinenko, A. et al. A spinwave ising machine. Communications Physics 6, 227 (2023). Page 35 of 38

  6. [15]

    & Åkerman, J

    González, V., Litvinenko, A., Khymyn, R. & Åkerman, J. Global biasing using a hardwarebased artificial zeeman term in spinwave ising machines. In 2023 IEEE International Magnetic Conference-Short Papers (INTERMAG Short Papers), 1–2 (IEEE, 2023)

  7. [16]

    & Åkerman, J

    Litvinenko, A., Khymyn, R., Ovcharov, R. & Åkerman, J. A 50-spin surface acoustic wave ising machine. arXiv preprint arXiv:2311.06830 (2023)

  8. [17]

    & Nakamura, Y

    Goto, H., Lin, Z. & Nakamura, Y. Boltzmann sampling from the ising model using quantum heating of coupled nonlinear oscillators. Scientific reports 8, 7154 (2018)

  9. [18]

    & Goto, H

    Kanao, T. & Goto, H. High -accuracy ising machine using kerr -nonlinear parametric oscillators with local four -body interactions. npj Quantum Information 7, 18 (2021)

  10. [19]

    Marandi, A., Wang, Z., Takata, K., Byer, R. L. & Yamamoto, Y. Network of time -multiplexed optical parametric oscillators as a coherent ising machine. Nature Photonics 8,937–942(2014)

  11. [20]

    Okawachi, Y. et al. Demonstration of chip -based coupled degenerate optical parametric oscillators for realizing a nanophotonic spin-glass. Nature Communications 11, 4119 (2020)

  12. [21]

    Integrated coherent ising machines based on self-phase modulation in micro ring resonators

    Tezak,N.etal. Integrated coherent ising machines based on self-phase modulation in micro ring resonators. IEEE Journal of Selected Topics in Quantum Electronics 26, 1–15 (2019)

  13. [22]

    & Kim, C

    Lo, H., Moy, W., Yu, H., Sapatnekar, S. & Kim, C. H. An ising solver chip based on coupled ring oscillators with a 48 - node all-to-all connected array architecture. Nature Electronics 6, 771–778 (2023)

  14. [23]

    Karpuzcu, U. et al. Cobi: A coupled oscillator based ising chip for combinatorial optimization. Nature Portfolio: Manuscript submitted for publication (2024)

  15. [24]

    K., Mallick, A

    Bashar, M. K., Mallick, A. & Shukla, N. Experimental investigation of the dynamics of coupled oscillators as ising machines. IEEE Access 9, 148184–148190 (2021)

  16. [25]

    Mallick, A. et al. Using synchronized oscillators to compute the maximum independent set. Nature communications 11, 4689 (2020)

  17. [26]

    K., Lin, Z

    Bashar, M. K., Lin, Z. & Shukla, N. Stability of oscillator ising machines: Not all solutions are created equal. Journal of Applied Physics 134 (2023)

  18. [27]

    K., Li, Z., Narayanan, V

    Bashar, M. K., Li, Z., Narayanan, V. & Shukla, N. An fpga -based max -k-cut accelerator exploiting oscillator synchronization model. In 2024 25th International Symposium on Quality Electronic Design (ISQED), 1–8 (IEEE, 2024)

  19. [28]

    Kalinin, K. P. & Berloff, N. G. Simulating ising and n-state planar potts models and external fields with nonequilibrium condensates. Physical review letters 121, 235302 (2018). Page 36 of 38

  20. [29]

    Polaritonic xy-ising machine

    Kalinin,K.P.,Amo,A.,Bloch,J.&Berloff,N.G. Polaritonic xy-ising machine. Nanophotonics 9, 4127–4138 (2020)

  21. [30]

    Luo, S. et al. Classical spin chains mimicked by room -temperature polariton condensates. Physical Review Applied 13, 044052 (2020)

  22. [31]

    & Datta, S

    Camsari, K. & Datta, S. Waiting for quantum computing? try probabilistic computing. IEEE Spectrum (2021)

  23. [32]

    Y., Faria, R., Sutton, B

    Camsari, K. Y., Faria, R., Sutton, B. M. & Datta, S. Stochastic p -bits for invertible logic. Physical Review X 7, 031014 (2017)

  24. [33]

    Liu, Y. et al. Time-division multiplexing ising computer using single stochastic magnetic tunneling junction. IEEE Transactions on Electron Devices 69, 4700–4707 (2022)

  25. [34]

    Yin, J. et al. Scalable ising computer based on ultra-fast field-free spin orbit torque stochastic device with extreme 1 -bit quantization. In 2022 International Electron Devices Meeting (IEDM), 36–1 (IEEE, 2022)

  26. [35]

    Lu, A. et al. Scalable in-memory clustered annealer with temporal noise of finfet for the travelling salesman problem. In 2022 International Electron Devices Meeting (IEDM), 22–5 (IEEE, 2022)

  27. [36]

    Kaiser, J. et al. Hardware-aware in situ learning based on stochastic magnetic tunnel junctions. Physical Review Applied 17, 014016 (2022)

  28. [37]

    Y., Sutton, B

    Camsari, K. Y., Sutton, B. M. & Datta, S. P-bits for probabilistic spin logic. Applied Physics Reviews 6 (2019)

  29. [38]

    & Datta, S

    Kaiser, J. & Datta, S. Probabilistic computing with p-bits. Applied Physics Letters 119 (2021)

  30. [39]

    Probabilistic and oscillatory Ising machines for combinatorial optimization: from software to hardware - acceleration with spintronics

    Grimaldi, A. Probabilistic and oscillatory Ising machines for combinatorial optimization: from software to hardware - acceleration with spintronics. Ph.D. thesis, Università degli Studi di Messina (2023)

  31. [40]

    Singh, N. S. et al. Cmos plus stochastic nanomagnets enabling heterogeneous computers for probabilistic inference and learning. Nature Communications 15, 2685 (2024)

  32. [41]

    & Yamaoka, M

    Okuyama, T., Hayashi, M. & Yamaoka, M. An ising computer based on simulated quantum annealing by path integral monte carlo method. In 2017 IEEE international conference on rebooting computing (ICRC), 1–6 (IEEE, 2017)

  33. [42]

    Zhang, T. et al. A review of ising machines implemented in conventional and emerging technologies. IEEE Transactions on Nanotechnology (2024)

  34. [43]

    & Han, J

    Zhang, T., Tao, Q., Liu, B. & Han, J. A review of simulation algorithms of classical ising machines for combinatorial optimization. In 2022 IEEE International Symposium on Circuits and Systems (ISCAS), 1877–1881 (IEEE, 2022)

  35. [44]

    Enhancing the solution quality of hardware ising -model solver via parallel tempering

    Gyoten, H., Hiromoto, M., & Sato, T. Enhancing the solution quality of hardware ising -model solver via parallel tempering. In 2018 IEEE/ACM International Conference on Computer-Aided Design (ICCAD) (pp. 1-8). (IEEE, 2018). Page 37 of 38

  36. [45]

    Ising formulations of many np problems

    Lucas, A. Ising formulations of many np problems. Frontiers in physics 2, 5 (2014)

  37. [46]

    & Virasoro, M

    Mézard, M., Parisi, G. & Virasoro, M. A. Spin glass theory and beyond: An Introduction to the Replica Method and Its Applications, vol. 9 (World Scientific Publishing Company, 1987)

  38. [47]

    Y., Faria, R., Sutton, B

    Camsari, K. Y., Faria, R., Sutton, B. M. & Datta, S. Stochastic p-bits for invertible logic. Physical Review X 7, 031014 (2017)

  39. [48]

    & Bermel, P

    Zhu, J., Xie, Z. & Bermel, P. Numerical simulation of probabilistic computing to np -complete number theory problems. Journal of Photonics for Energy 13, 028501–028501 (2023)

  40. [49]

    Aadit, N. A. et al. Massively parallel probabilistic computing with sparse ising machines. Nature Electronics 5, 460–468 (2022)

  41. [50]

    H., Hinton, G

    Ackley, D. H., Hinton, G. E. & Sejnowski, T. J. A learning algorithm for boltzmann machines. Cognitive science 9, 147– 169 (1985)

  42. [51]

    Whitehead, W., Nelson, Z., Camsari, K. Y. & Theogarajan, L. Cmos -compatible ising and potts annealing using single - photon avalanche diodes. Nature Electronics 6, 1009–1019 (2023)

  43. [52]

    Inaba, K. et al. Potts model solver based on hybrid physical and digital architecture. Communications Physics 5, 137 (2022)

  44. [53]

    K., Hasan, A

    Bashar, M. K., Hasan, A. & Shukla, N. Designing a k-state p-bit engine. arXiv preprint arXiv:2403.06436 (2024)

  45. [54]

    Tamura, K. et al. Performance comparison of typical binary-integer encodings in an ising machine. IEEE Access 9, 81032- 81039 (2021)

  46. [55]

    Winston, W. L. Operations research: applications and algorithms (Cengage Learning, 2022)

  47. [56]

    & Klein, J

    Wisniewski, M. & Klein, J. H. Linear programming. Critical path analysis (Palgrave, 2001)

  48. [57]

    Vella, D. C. Invitation to linear programming and game theory (Cambridge University Press, 2021)

  49. [58]

    & Yang, C

    Yang, Z. & Yang, C. How far can a biased random walker go? Journal of Applied Mathematics and Physics 3, 1159– 1167 (2015)

  50. [59]

    https://caravel-mgmt-soc-litex.readthedocs.io/en/latest/ (2022)

    Caravel management soc - litex. https://caravel-mgmt-soc-litex.readthedocs.io/en/latest/ (2022). Accessed: 2024-06-25

  51. [60]

    General Algebraic Modeling System (GAMS) Release 49.3.0 (2025)

    GAMS Development Corporation. General Algebraic Modeling System (GAMS) Release 49.3.0 (2025)

  52. [61]

    The Italian factory of micromagnetic modelling and spintronics

    O’Neill, M. E. Pcg: A family of simple fast space -efficient statistically good algorithms for random number generation. Tech. Rep. HMC-CS-2014-0905 (2014). Page 38 of 38 Acknowledgments: This work was supported by the U.S. National Science Foundation (NSF) under award numbers...

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.