SE2250979A1 - Spinwave coherent ising machine - Google Patents

Abstract

A spinwave Coherent Ising computational machine comprising: a spinwave delay line that propagates a plurality of pseudo Ising spin pulses in the form of parametric spinwave oscillators, a spinwave-microwave transducer, an electronic phase-sensitive amplifier configured to cause phase degeneracy of a plurality of propagating spinwave pulses, a measuring unit configured to temporarily measure pseudo spins of the plurality of pseudo spin pulses, an interaction unit configured to temporarily implement magnitudes and phases of all the interactions relating to the plurality of pseudo Ising spin pulses; and, a pseudo Ising spin measuring unit which measures the pseudo Ising spins of the plurality of pseudo Ising spin pulses.

Description

The disclosed invention generally relates to non-von-Neuman computing architectures and more particularly to a combinatorial spinwave Ising-model solver using the physical annealing method.

Background of the invention The Ising Model [Reference: l] is an efficient computational tool that can be used to solve a variety of difficult computational problems time- and hardware-efficiently by using its physical implementation - Ising machine. Ising model manifests that if there is a magnetic structure that is made of an array of magnetic domains and each domain is either up or down Ci, cj: +1 01" - 1 magnetic spin and they are coupled to each other through a magnetic field with a coefficient 11-1-, one can calculate the energy of the whole system by simple summation of the product of each spin state in the pairs and their coupling: E = _ 2131111019 (1) The Ising problem is to find for a given magnetic structure with particular coupling between elements, a unique configuration of spins (up or down) so that the whole magnetic structure has the lowest energy state. The Ising problems belong to the category of so-called NP-hard computational problems [2]. The term hard means that this problem is representative of the whole class of NP and can be mapped to an Ising machine and then solved with polynomial time.

An Ising machine can be implemented with a physical system that is made of an array of elements where each element has two well-defined and stable states that can be either +1 or -1 to represent an Ising spin and where each element can be connected to any other with a Variable continuous coefficient. By setting different coefficients one can program this machine to solve a particular Ising problem. An Ising machine when tumed on tends to go into the lowest energy state by flipping the states of individual spins. The final configuration of the spins coincides with the ground state of the Ising problem and represents its solution.

To date, Ising machines have been implemented with many physical paradigms, including quantum annealing [3-4], optical parametric oscillators [5-7], phase transition nano-oscillators [8], stochastic nanomagnets [9], electronic CMOS SRAM [l0], electronic LC oscillators [l l, l2], and spin-Hall nano-oscillators [l3-l4]. All the concepts are characterized by different speeds, power consumption, number of supported spins, physical dimension, etc. but can still be divided into two distinct groups - spatially distributed oscillator arrays and time-multiplexed soliton systems as shown in figure lA and lB.

The most important parameters for Ising Machines are the time to solution and the number of supported spins, and these parameters are strongly interconnected. In FIG. lC we present time to solution as a function of annealing time and problem size N. For Ising Machines based on physical arrays of oscillators the main problem is interconnectivity because as the number of oscillators N grows the number of intersections between coupling lines increases dramatically as - O(N2). This problem is solved by grouping the elements with a sparse connection into a so-called chimera graph [l6]. However, it trades off the computational time-to-solution which in the case of the chimera graph connection scheme increases as - O(N2) (FIG.lC). It is the same growth rate as the computation speed of classical computers based on a von-Neuman architecture which means that there is no principal computational advantage in using Ising Machines built with spatially distributed oscillators.

For Coherent Ising Machines based on propagating light pulses the interconnectivity for a problem of all-to-all connected spins is easily solved with the time- multiplexing method. Computational problems with all-to-all connected spins are characterized by the computational time that grows as - O(N) when solved on physical Ising Machines. Therefore, the time multiplexing method makes computational time to solution for Ising problems with a large (>5 0) number of spins reasonable. The first time- multiplexing Ising Machine was implemented with optical parametric oscillators (OPOs) that are in the form of propagating light pulses in optical waveguides. The interconnectivity is implemented electrically (FIG.lB) by consecutive measurements of each propagating light pulses-OPOs and then adding to the additional small (r << 1) contributions according to the coefficients in Ising problem: fi = Tzjfijcj (2) 111 111 1:1 = 5 5 (3) 1:1 fij For the moment, time-multiplexed CIMs seem to be the most promising Ising Machine configuration for combinatorial problems with large (> 50 number of elements) due to -O (N) computational speed. Even SHNO-based Ising Machines that are proj ected to provide unpresidential computational speed starting from tens of nanoseconds cannot compete with Coherent Ising Machines for a graph size above 50 as they belong to -O(N2) class of IM. The number of supported elements in Ising problems solved by optical CIMs was progressively growing from 100 [5] spins through 2000 [6] and recently a CIM supporting 100000 spins was reported [7]. However, despite clear progress in the number of supported spins, current CIMs are still not in the market as they have a considerable disadvantage of size, large power consumption, and, most importantly, costly optical infrastructure requiring optical tables, precise positioners, etc.

Summarv of the invention Accordingly, the present invention preferably seeks to mitigate, alleviate or eliminate one or more of the above-identified deficiencies in the art and disadvantages singly or in any combination and solves at least the above mentioned problems by providing a spinwave Coherent Ising computational machine, SCIM, comprising: a spinwave delay line that propagates a plurality of pseudo Ising spin pulses in the form of parametric spinwave oscillators; a spinwave-microwave transducer that transforrns between spinwaves and electrical pseudo spin pulses; an electronic phase-sensitive amplif1er conf1gured to cause phase degeneracy of a plurality of propagating spinwave pulses corresponding to a plurality of pseudo spins in the Ising model, and set an oscillation frequency corresponding to half a reference frequency; a measuring unit configured to temporarily measure pseudo spins of the plurality of pseudo spin pulses by temporarily measuring the phase of each of the plurality of pseudo spin pulses each time the plurality of pseudo spin pulses circularly propagates in a ring circuit and suspends measurement after one set of measurement is completed before another set of measurement is restarted; an interaction unit conf1gured to temporarily implement magnitudes and phases of all the interactions relating to the plurality of pseudo Ising spin pulses; and, an annealing and computing unit which perforrns annealing by restarting the circulation of the plurality of pseudo Ising spin pulses and perforrns statistical analysis on the plurality of steady-state solution.

This disclosure relates to a novel architecture of Coherent Ising Machines (CIM) for solving combinatorial optimization problems. The key element of Spinwave CIM (SCIM) is a YIG spinwave delay line that is used as a waveguide where spinwave RF pulses propagate and are stored. The advantage of spinwave devices is the exceptionally slow group velocity of propagating spinwaves that is several orders lower than the speed of light. It allows miniaturizing of the CIM waveguide size down to the mm scale while keeping a high number of supported spins. Another advantage of the proposed invention in contrast to optical CIMs is that spinwaves can also be easily excited by and transduced back to electrical RF signals, which allows performing the compensation of propagation losses via low-power and power-efficient RF phase-sensitive and linear amplifiers. The transformation between spinwave RF pulses to electric RF pulses and back is implemented with wideband thin-wire transducers. The rapid measurement of RF circulating pulses is done electrically via homodyne IQ-demodulator and CORDIC algorithm implemented in microcontroller or FPGA that allows independent deduction of the phase and the amplitude of the RF pulses within pulse circulation time.

Further advantageous embodiments are disclosed in the appended and dependent patent claims.

Brief description of the drawings These and other aspects, features and advantages of which the invention is capable will be apparent and elucidated from the following description of embodiments of the present inVention, reference being made to the accompanying drawings, in which Fig. 1A. is a spatially distributed oscillator array Ising Machine.

Fig. 1B is a time-multiplexed Ising Machine.

Fig. 1C is a graph showing time to solution as a function of annealing time and problem size N, showing a comparison between existing solutions to the present disclosure, which is shown as the dashed line labelled Coherent Ising Machine Spinwaves GU.

Fig. 2A is a block diagram of a Spinwave Coherent Ising Machine comprising microwave delay lines and a microcontroller according to an aspect.

Fig. 2B is a block diagram of a Spinwave Coherent Ising Machine comprising an FPGA according to an aspect.

Fig. 3A shows the S21-parameter of a YIG spinwave delay line.

Fig. SB shows the delay time of the YIG spinwave delay line.

Fig. 4 shows the amplification of the phase-sensitive amplification block as a function of phase difference between circulating RF pulses and a reference signal.

Fig. 5 shows time traces of circulating RF pulses, their instantaneous phase and a control signal for RF switches.

Fig. 6 is a schematic of the transition between a non-optimal initial state and an optimal solution for MAX-CUT4 problem computed with the SCIM after 12 circulations.

Fig. 7A shows time traces of the circulating RF pulses and their instantaneous phase for the 1st circulation period.

Fig. 7B shows time traces of the circulating RF pulses and their instantaneous phase for an optimal solution computed at the 12th circulation period.

Fig. 8 shows the evolution of time traces and instantaneous phase for the Std RF pulse within interrnediate circulation steps at top and middle panel, correspondingly. Bottom panel shows the evolution of the instantaneous phase at the center of 1st, 2ttd, Std and 4th RF pulse withing interrnediate circulation steps.

Fig. 9 shows a SCIM comprising a plurality of ring circuits, each ring Circuit connected to an FPGA.

Detailed description Figure 2A and 2B show a spinwaVe coherent Ising machine (SCIM). The key element of the SCIM is a spinwaVe delay line 21, 23 that supports up to N Ising spins in the forrn of propagating spinwaVe RF pulses. The spinwaVe delay line consists of a ferromagnetic material, such as YIG, spinwaVe waVeguide 21 with an input 20 and an output 22 electromagnetic-spinwaVe transducers. The operating frequency fo of SCIM is set by phase-sensitiVe amplification block 1 Via a reference frequency signal fref = fo 2. The losses in YIG spinwaVe waVeguide are compensated by the phase-sensitiVe block and a linear amplifier 7. The phase of each circulating RF pulse cj is measured by deflecting 10% of power after a filter 6 with a directional coupler 8, power diVider 9 and measuring the phase with an IQ-demodulator 10, 2-channel ADC 11, and a CORDIC algorithm implemented in microcontroller 13 as shown in f1g. 2A.

The cable delay lines delay each circulating RF pulse by multiple of pulse repetition time introducing the coupling between Ising spins. The microcontroller 13 perforrns the computation of the Ising matrix (Equation 2,3) and sets the coupling between each Ising spin by controlling RF switches 16 in every delay channel, a phase shifter 19 and a Variable amplifier 25 which change the phase and amplitude of additional coupling RF pulses. Coupling RF pulses after a Variable amplifier 25 are added to the circulating RF pulses Via power coupler 21. The microcontroller 13 controls RF switches Via control lines 12, perforrns an annealing procedure, and communicates with extemal systems via data port 14.

As shown in the SCIM in figure 2B, the CORDIC algorithm may be implemented in an FPGA. In such an architecture, the FPGA also perforrns the computation of the Ising matrix and sets the coupling between each Ising spin by controlling a phase shifter 17 and Variable amplifier 18 which change the phase and the amplitude of additional coupling RF pulses that are formed by the RF switch 16. Coupling RF pulses after an Variable amplifier 18 are added to the circulating RF pulses Via power coupler 19.

As shown in figure 2B, the SCIM comprising an FPGA comprises a microcontroller. The microcontroller ll controls the FPGA, clock frequencies 9, perforrns an annealing procedure, and optionally communicates with external systems.

The number N of supported spins is proportional to the total delay time rdelay in a spinwave delay line and inversely proportional to the minimum possible spinwave RF pulse width Ip: N = Td-:lay 17 The minimum possible spinwave RF pulse duration is limited by the largest value derived from the 3-dB bandwidth Bl/l/SW of the spinwave spectra and the delay time deviation Ardelay within BVVSW: 1 BWSW ' (5) ATdelay Ipä A particular Ising matrix (Equation 3) can be mapped into a Spinwave Coherent Ising Machine by setting an arbitrary time domain pattem of control signal for RF switches 16 so that each RF propagating pulse can be connected to any other RF pulses. In an architecture comprising cable delay lines, the number of cable delay channels should be M = N - 1 to support all-to-all connections. The delay time in cable delay lines CDL1 17 to CDL M is proportional to their index: TTL (6) TcDL m = ï pu ses where fpulses is a pulse repetition frequency.

FIG.3A shows the spinwave spectra of the YIG spinwave delay line in the form of Szl-parameter. The bandwidth of the spinwave generation spectrum is 60MHz measured at -3 dB level. FIG.3B shows the delay time of the YIG spinwave delay line [23]. The mean delay time Idelay is 270.62 ns. The limit of minimal pulse duration derived from bandwidth is 16.75 ns while the delay time deviation imposes a stronger limitation of 30.l9ns that results in 8 supported Ising spins , rounded down from 8.9.

Phase sensitivity is achieved by doubling the reference signal frequency via a frequency doubler 3 and combining it with an RF signal converted from spinwave RF pulses. The amplitude of the total signal after power divider/coupler 4 is set at a level that is close to the saturation point of an RF amplifier 5. The amplitudes of the signals at two inputs of a power divider/coupler 4 which have a relative phase close to 0° or 180° are added and the total signal amplification is affected by the saturation of the amplif1er 5 more than if signals have a phase difference is close to 90° or 360°. The signal after an amplif1er 5 is f1ltered by a highpass filter 6 with a cut of frequency fo to remove 2 fo signal after phase-sensitive amplif1cation.

FIG. 4 shows the amplification of the phase-sensitive amplif1cation block 1 as a function of a phase difference between the circulating RF pulses and the reference signal 2. The difference in amplification is denoted as APSA, which equals 6dB and represents the phase sensitivity of the block. The value of phase sensitivity can be adjusted by changing the amplitude of the reference signal at the input of the frequency multiplier.

FIG.5 shows time traces of control signal for RF switches, RF signal after RF amplif1er 7, and instantaneous phase that is calculated from RF signal Vampr RF switch 26 is used to counteract the spreading of propagating RF pulses due to the delay time deviation Ardelay of the YIG spinwave waveguide. tå?,% is the switching time from the rising front of control signal CLSW 26 to the moment when RF switch 16 is open for 90%. Similarly, tåg? is the switching time from CLSW 26 falling front to the moment when RF switch 16 is open for 10%. RF switch 16 is used to form coupling RF pulses that are added to the circulating RF pulses and has similar parameters tå?,% and tåg? Fig. 9 shows a SCIM comprising a plurality of ring circuits. Each ring circuit comprises a respective spinwave delay line 21, and a phase sensitive amplifier 1. Each ring circuit is connected to the FPGA 15 forrning at least the annealing and computing unit. Due to the reduced physical and power footprint of the present ring circuit, multiple ring circuits are practically feasible in a single SCIM. Each ring circuit may comprise a respective electric linear amplifier. The plurality of ring circuits enables an increased number of spins to be performed by the SCIM.

Examples 4 spin MAX-CUT solution The following illustrative example demonstrates the physical mechanisms and routes for obtaining a solution to a simple 4-spin MAX-CUT problem and is representatiVe of embodiments of the schematic design, the physical parameters, and methods described herein are not meant to be limiting.

FIG. 6 shows a 4-spin MAX-CUT problem with nearest-neighbour connections which is described by the following Ising matrix: 0-100 00-10 j" I 0 0 0 -1 (7) -1 0 0 0 The initial state of the system is chosen randomly to be cj = {+1 -1 -1 -1}. Value +1 corresponds to 180° of the phase difference between the RF pulse and the reference signal 2, while value -1 corresponds to 0° of phase difference.

The problem (7) is mapped into the SCIM using a single cable delay line with a delay time that according to (6) equals to a pulse repetition period. The RF switch 16 is always open so that each RF pulse c,- is directionally coupled to an RF pulse cl-fl Via a single cable delay. In initial state, additional coupling RF pulses have the following values f,- = {+1 of additional delayed coupling RF pulses according to table 1 below: -1 +1 +1}. Controlling components 19, 20 set the amplitude and phase f,- Amplitude Phase + 1 1 180° - 1 1 0° 0 0 1 80° or 0° TABLE 1 The measurement block composed of components 8, 9, 10, ll, 13 reads the instantaneous phase at the centre of each RF pulse, rounds its value, and processes the data as an interrnediate state of the Ising system.

FIG.6 shows schematically the initial state that corresponds to the edge cut number of 2 which is a non-optimal solution for a MAX-CUT problem defined by (Eq. 7). After 12 circulations SCIM evolves to a state 61- final = {+1 -1 +1 -1} which represents an optimal solution. The microcontroller detects it as a stable state and sends 61- final to the data port 14.

FIG.7 shows the time traces of circulating RF pulses and computed instantaneous phase for the 1st and 12th circulation periods. The phase of spin 3 has changed from 0° (63 = -1) to 180° (63 = +1).

FIG.8 shows the evolution of times traces and computed instantaneous phase for the 3fd RF pulse within interrnediate circulations steps. The instantaneous phase changes nonuniforrnly forrning a temporal "domain wall" whose interrnediate phase forces the formation of a dip in RF amplitude in the middle of the RF pulse which gradually propagates to the left and disappears leading to a uniforrnly distributed instantaneous phase of 180° (63 = +1).

Although, the present invention has been described above with reference to specific embodiments, it is not intended to be limited to the specific form set forth herein. Rather, the invention is limited only by the accompanying claims.

In the claims, the term "comprises/comprising" does not exclude the presence of other elements or steps. Furthermore, although individually listed, a plurality of means, elements or method steps may be implemented by e.g. a single unit or processor. Additionally, although individual features may be included in different claims, these may possibly advantageously be combined, and the inclusion in different claims does not imply that a combination of features is not feasible and/or advantageous. In addition, a: aa a: a, singular references do not exclude a plurality. The terms , an", "f1rst", "second" etc do not preclude a plurality. Reference signs in the claims are provided merely as a clarifying example and shall not be construed as limiting the scope of the claims in any way. 11 List of references used in the figures: Reference Numbers (Fig 2A & 2B unless otherwise indicated) Element 1 Phase sensitive aniplifier block 2 Reference signal 3 Frequency doubler 4 Power divider 1:1 RF aniplifier 6 Highpass filter 7 RF Aniplifier 8 RF coupler 1:10 9 (Fig 2A only) Power divider 1:1 9* (Fig 2B only) Clock line (Fig 2A only) I and Q den1odulator * (Fig 2B only) Control line 11 (Fig 2A only) Analog to Digital Converter 11* (Fig 2B only) Microcontroller 12 (Fig 2A only) Control line 12* (Fig 2B only) Data line 13 (Fig 2A only) Microcontroller 13* (Fig 2B only) I and Q den1odulator 14 (Fig 2A only) Data line 14* (Fig 2B only) Analog to Digital Converter (Fig 2A only) Multichannel RF splitter * (Fig 2B only) Field-Programmable Gate Array 16 RF switch 17 (Fig 2A only) Cable delay line 17* (Fig 2B only) Digital phase shifter 18 (Fig 2A only) Multichannel RF conibiner 18* (Fig 2B only) Digital attenuator 19 (Fig 2A only) Digital phase shifter 19* (Fig 2B only) RF coupler 1:10 Input electroniagnetic-to-spinwave transducer 21 YIG spinwave waveguide 22 Output spinwave-to-electroniagnetic transducer 23 Spinwave delay line 24 RF switch 25 Digital variable aniplifier 26 RF coupler 1:10 12 Literature References 1. S. Kirkpatrick, C. D. Gelatt, M. P. Vecchi, Science 220, 671 (1983). Publisher: American Association for the Advancement of Science 2. F. Barahona, Journal of Physics A: Mathematical and General 15, 3241 (1982). 3. Johnson, M., Amin, M., Gildert, S. et al. Quantum annealing With manufactured spins. Nature 473, 194-198 (2011). https://doi.org/10.1038/nature10012 4. Davide Venturelli, Salvatore Mandra, Sergey Knysh, Bryan O"Gorman, Rupak BisWas, and Vadim Smelyanskiy. Phys. Rev. X 5, 031040 - Published 18 September 2015 . McMahon, Peter L., Alireza Marandi, Yoshitaka Haribara, Ryan Hamerly, Carsten Langrock, Shuhei Tamate, Takahiro Inagaki et al. "A fully programmable 100-spin coherent Ising machine With all-to-all connections." Science 354, no. 6312 (2016): 614- 617. 6. Inagaki, Takahiro, Yoshitaka Haribara, Koji Igarashi, Tomohiro Sonobe, Shuhei Tamate, Toshimori Honjo, Alireza Marandi et al. "A coherent Ising machine for 2000- node optimization problems." Science 354, no. 6312 (2016): 603-606. 7. Honj o, Toshimori, Tomohiro Sonobe, Kensuke Inaba, Takahiro Inagaki, Takuya Ikuta, Yasuhiro Yamada, Takushi Kazama et al. "100,000-spin coherent Ising machine." Science advances 7, no. 40 (2021): eabh0952. 8. Dutta, S., A. Khanna, J. Gomez, K. Ni, Z. Toroczkai, and S. Datta. "Experimental demonstration of phase transition nano-oscillator based Ising machine." In 2019 IEEE Intemational Electron Devices Meeting (IEDM), pp. 37-8. IEEE, 2019. 9. Sutton, Brian, Kerem Yunus Camsari, Behtash Behin-Aein, and Supriyo Datta. "Intrinsic optimization using stochastic nanomagnets." Scientific reports 7, no. 1 (2017): 1-9.

. M. Yamaoka, C. Yoshimura, M. Hayashi, T. Okuyama, H. Aoki and H. Mizuno, "A 20k-Spin Ising Chip to Solve Combinatorial Optimization Problems With CMOS Annealing," in IEEE Joumal of Solid-State Circuits, vol. 51, no. 1, pp. 303-309, Jan. 2016, doi: 10.1109/JSSC.2015.2498601. 11. Tianshi Wang, Leon Wu, and Jaijeet Roychowdhury. 2019. New Computational Results and Hardware Prototypes for Oscillator-based Ising Machines. In Proceedings of the 56th Annual Design Automation Conference 2019 (DAC '19). Association for 13 Computing Machinery, New York, https://doi.org/10.1145/3316781.3322473. 12. Wang, Tianshi and J aij eet S. Roychowdhury. "OIM: Oscillator-based Ising Machines for Solving Combinatorial Optimisation Problems." UCNC (2019).

NY, USA, Article 239, 1-2. 13. Albertsson, Dagur Ingi, Mohammad Zahedinej ad, Afshin Houshang, Roman Khymyn, Johan Åkerman, and Ana Rusu. "Ultrafast Ising Machines using spin torque nano-oscillators." Applied Physics Letters 118, no. 11 (2021): 112404. 14. McGoldrick, Brooke C., Jonathan Z. Sun, and Luqiao Liu. "Ising machine based on electrically coupled spin Hall nano-oscillators." Physical Review Applied 17, no. 1 (2022): 014006.

Claims (10)

1.Claims A spinWaVe Coherent Ising computational machine, SCIM, comprising: a spinwave delay line (23) that propagates a plurality of pseudo Ising spin pulses in the form of parametric spinwave oscillators; a spinwave-microwave transducer (20, 22) that transforrns between spinwaves and electrical pseudo spin pulses; an electronic phase-sensitive amplifier (1) configured to cause phase degeneracy of a plurality of propagating spinwave pulses corresponding to a plurality of pseudo spins in the Ising model, and set an oscillation frequency corresponding to half a reference frequency; a measuring unit (8, 9, 10, 11, 11* 13, 13*,14*,15*) conf1gured to temporarily measure pseudo spins of the plurality of pseudo spin pulses by temporarily measuring the phase of each of the plurality of pseudo spin pulses each time the plurality of pseudo spin pulses circularly propagates in a ring circuit and suspends measurement after one set of measurement is completed before another set of measurement is restarted; an interaction unit (11, 16, 17, 18, 19, 15*,17*, 18*, 19*, 25, 26) configured to temporarily implement magnitudes and phases of all the interactions relating to the plurality of pseudo Ising spin pulses; and, an annealing and computing unit (11*, 13, 15*) Which perforrns annealing by restarting the circulation of the plurality of pseudo Ising spin pulses and perforrns statistical analysis on the plurality of steady-state solution.

2.The SCIM according to claim 1, Wherein the spinwave delay line (23) comprises a ferromagnetic Wave propagating material in Which the spin pulses propagate at a speed substantially less than the speed of light.

3.The SCIM according to claims 1 or 2, comprising an electric linear amplif1er (5) configured to compensate the amplitude losses of the plurality of propagating spinWaVe pulses.

4.The SCIM according to any of claims 1 to 3, Wherein the ferromagnetic Wave propagating material comprises the phase sensitive amplifier (1) or Wherein the SCIM comprises an additional ferromagnetic Wave propagating resonator forrning the phase sensitive amplif1er (1).

5.The SCIM according to any of claims 1 to 4, Wherein the machine comprises a FPGA (15*) conf1gured to map an Ising problem onto propagating spinWaVes by connecting them Via applying to each RF pulse, an additional coupling RF pulse With an amplitude and phase defined by a coupling matrix.

6.The SCIM according to claim 5, Wherein the measuring unit receives the propagated pseudo spin pulse after the measuring unit completes one set of measurement and before the measuring unit restarts another set of measurement.

7.The SCIM according to any of claims 1 to 6, Wherein the parametric oscillators are time-multiplexed in a single spinWaVe delay line.

8.The SCIM according to any of claims 1 to 4, Wherein the plurality of propagating spinwave pulses are connected to the interaction unit via a plurality of cable delay lines according to the Ising model.

9.The SCIM according to any of claim 1 to 8, Wherein the spinWaVe coherent Ising machine is configured to produce a set of numbers from a plurality of the Ising spins, Wherein the produced numbers contain information about statistics of stable states of the spinWaVe coherent Ising computational machine, and Wherein the produced numbers are generated from the phases of a plurality of pseudo spin pulses circularly propagating in the ring circuit.

10. The SCIM according to any of c1ain1s 1 to 9, Wherein the SCIM coniprises a p1ura1ity of ring circuits, each ring circuit coniprising a respective spinWaVe delay line (23) and phase-sensitive aniplifier (1).