WO2017033337A1 - 情報処理方法及び装置 - Google Patents

情報処理方法及び装置 Download PDF

Info

Publication number
WO2017033337A1
WO2017033337A1 PCT/JP2015/074263 JP2015074263W WO2017033337A1 WO 2017033337 A1 WO2017033337 A1 WO 2017033337A1 JP 2015074263 W JP2015074263 W JP 2015074263W WO 2017033337 A1 WO2017033337 A1 WO 2017033337A1
Authority
WO
WIPO (PCT)
Prior art keywords
analog
circuit
processing
spin
data
Prior art date
Application number
PCT/JP2015/074263
Other languages
English (en)
French (fr)
Inventor
地尋 吉村
拓哉 奥山
真人 林
山岡　雅直
亮仁 赤井
Original Assignee
株式会社日立製作所
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 株式会社日立製作所 filed Critical 株式会社日立製作所
Priority to PCT/JP2015/074263 priority Critical patent/WO2017033337A1/ja
Priority to JP2017536156A priority patent/JP6636529B2/ja
Publication of WO2017033337A1 publication Critical patent/WO2017033337A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F7/00—Methods or arrangements for processing data by operating upon the order or content of the data handled
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N99/00—Subject matter not provided for in other groups of this subclass

Definitions

  • the present invention relates to an analog computer that realizes a problem to be solved in an information processing apparatus as an apparatus having a structure corresponding thereto.
  • Emergence is the appearance of a property that is not just a simple sum of the properties of the part.
  • the behavior formed by local synaptic connections between neurons realizes some kind of information processing as a whole (for example, Patent Document 3).
  • Examples of emergent calculation models include genetic algorithms and swarm intelligence.
  • many calculation models based on local interactions can be said to be emergent.
  • the above-described analog computer has attracted attention again for realizing an emergent calculation model (for example, Patent Document 4).
  • Patent Document 5 describes a Ising model using the fact that a spin arrangement in a state with the minimum energy is regarded as a kind of optimization problem in a model that explains the behavior of a magnetic body by statistical physics called the Ising model.
  • This is a kind of analog computer in which a CMOS circuit having a similar structure is formed on a semiconductor substrate.
  • An analog computer is configured by arranging a number of elements corresponding to the components of a calculation model on a space (for example, on a semiconductor substrate when realized as a semiconductor element). For this reason, variations in the manufacture of the elements become a problem.
  • a neural network even if the elements simulate the same neuron, the characteristics differ from element to element due to variations in manufacturing.
  • the behavior is affected by individual differences in the input / output characteristics of the device. Therefore, a circuit that compensates for individual differences between elements has been added (for example, Patent Document 6).
  • one of information processing apparatuses executes a step of measuring and recording variation in characteristics of an analog computer in an information processing apparatus including one or more analog computers. It is an information processing apparatus characterized by having a function.
  • a numerical value storage memory that stores a numerical value to be processed or a result of the processing, a coefficient storage memory that stores a coefficient used for the processing, and at least a part of the processing using the numerical value and the coefficient
  • the information processing apparatus includes a characteristic data storage memory that stores characteristic data indicating characteristic variations of a plurality of analog processing circuits. This information processing apparatus stores problem data as numerical values or coefficients in a numerical value storage memory and a coefficient storage memory, and reads out the numerical values stored in the numerical value storage memory as solution data after one or more processes by the analog processing circuit.
  • the problem data before storing the problem data, the problem data is converted based on the characteristic data, and after the solution data is read, the solution data is converted based on the characteristic data.
  • a first numerical value storage memory a first coefficient storage memory, a first analog processing circuit, a second numerical value storage memory, and a second coefficient storage memory And a second analog processing circuit.
  • the second analog processing circuit performs processing using the numerical value stored in the first numerical value storage memory and the coefficient stored in the second coefficient storage memory, and the processing result is stored in the second numerical value storage memory.
  • Another example of the specific configuration of the present invention includes a processing device, a storage device, an input device, and an output device.
  • the input device receives the problem data, and stores the problem data received by the processing device in the storage device.
  • the processing device reads the problem data stored in the storage device, converts the read problem data based on the characteristic data, and stores it in the numerical value storage memory and coefficient storage memory of the calculation unit.
  • the processing device reads solution data from the numerical value storage memory after a plurality of processes by the analog processing circuit of the calculation unit, converts the solution data read based on the characteristic data, and outputs the solution data to the storage device or the output device.
  • Another example of the specific configuration of the present invention is to measure the characteristic variation of a plurality of analog processing circuits under the control of the processing device, and store the measurement result as characteristic data in the characteristic data storage memory.
  • the characteristic data storage memory may use a part of the above storage device.
  • a dedicated memory may be provided.
  • the position of the characteristic data storage memory in the device is not limited.
  • Another example of a specific configuration of the present invention includes an inspection circuit that inspects the processing contents of a plurality of analog processing circuits, and a digital processing circuit that digitally executes processing equivalent to the processing by the analog processing circuits. Then, the processing of the analog processing circuit is replaced with the digital processing circuit based on the inspection result of the inspection circuit.
  • the digital processing circuit may use a part of the above processing device.
  • a dedicated logic circuit may be provided. The position of the digital processing circuit within the device is not limited.
  • the inspection circuit As an example of the inspection by the inspection circuit, two analog values generated in the processing of the analog processing circuit are inspected, and when the difference between the analog values is smaller than the first threshold and larger than the second threshold, the analog There is one that replaces the processing of the processing circuit with the digital processing circuit.
  • the problem data is problem data for obtaining the spin arrangement of the Ising model or problem data for obtaining the state of the neurons of the neural network.
  • Another aspect of the present invention is a numerical value storage memory that stores a numerical value that is a processing target or a processing result, a coefficient storage memory that stores a coefficient used for processing, and at least a part of processing using the numerical value and the coefficient Is an information processing apparatus including a plurality of analog processing circuits that execute analog processing.
  • the information processing apparatus includes an inspection circuit that inspects the processing contents of the plurality of analog processing circuits, and a digital processing circuit that digitally executes processing equivalent to the processing by the analog processing circuit. Based on the result, the processing of the analog processing circuit is replaced with a digital processing circuit.
  • the inspection by the inspection circuit inspects the difference between two analog values (for example, voltage values) generated in the processing of the analog processing circuit, and the difference between the analog values is predetermined.
  • the processing of the analog processing circuit is replaced with the digital processing circuit.
  • the digital processing circuit a general-purpose microprocessor may be used, or a dedicated logic circuit may be used.
  • the method in an information processing method using an information processing apparatus including one or more analog computers, includes measuring a characteristic variation of the analog computer and recording it as characteristic variation data. Is an information processing method.
  • the analog computer obtains the spin arrangement of the Ising model, and the conversion of the input data inverts the sign of the interaction coefficient of the Ising model or the value of the external magnetic field coefficient.
  • the interpretation of the output data is the inversion of the spin value.
  • the analog computer has an inspection unit that performs the evaluation for each structural unit for evaluating the influence of the characteristic variation, and malfunctions by reading the inspection result from the inspection unit. Detect possible structural units.
  • a logic circuit included in the information processing apparatus temporarily substitutes the processing performed by the structural unit.
  • FIG. 3 is a block diagram for explaining an example of the configuration of unit components constituting the analog computer shown in FIG. 2. It is a model figure explaining an example of the calculation model which the analog computer shown in FIG. 2 simulates.
  • FIG. 5 is a model diagram for explaining a correspondence relationship between unit components shown in FIG. 3 and a calculation model shown in FIG. 4.
  • FIG. 4 is a circuit diagram illustrating an example of a configuration of a circuit that is used in the unit components of FIG. 3 and that causes manufacturing variations.
  • FIG. 3 is a block diagram explaining an example of a configuration of a circuit that is used in the unit components of FIG. 3 and that causes manufacturing variations.
  • FIG. 7 is a circuit diagram illustrating an example of a configuration of a circuit used in combination with the circuit of FIG. 6. It is a time chart explaining an example of the operation
  • 2 is a flowchart illustrating an example of operations performed by the information processing apparatus of FIG. 1 in order to realize information processing while avoiding manufacturing variations.
  • FIG. 2 is a conceptual diagram illustrating an example of data recorded with characteristic variations held by the information processing apparatus of FIG. 1.
  • 2 is a flowchart for explaining an example of an operation performed by the information processing apparatus of FIG. 1 in order to realize measurement of characteristic variation.
  • FIG. 3 is a model diagram for explaining an example of data input to the analog computer of FIG. 2 in order to realize measurement of characteristic variation.
  • FIG. 14 It is a block diagram explaining an example of a structure of the analog computer to which this invention is applied. It is a circuit diagram explaining an example of a structure of the test
  • inspection circuit shown in FIG. 14 is a flowchart illustrating an example of a procedure for the information processing apparatus of FIG. 1 to control the analog computer of FIG. It is a circuit diagram explaining another example of a structure of the test
  • notations such as “first”, “second”, and “third” are attached to identify the constituent elements, and do not necessarily limit the number or order.
  • a number for identifying a component is used for each context, and a number used in one context does not necessarily indicate the same configuration in another context. Further, it does not preclude that a component identified by a certain number also functions as a component identified by another number.
  • the information processing apparatus 100 includes an Ising model computer 200 (FIG. 2) having a structure similar to the Ising model.
  • the information processing apparatus 100 realizes information processing by avoiding manufacturing variations of MOSFET (Metal-Oxide-Semiconductor Field-Effect-Transistor) groups (FIG. 6) constituting the Ising model computer 200. .
  • MOSFET Metal-Oxide-Semiconductor Field-Effect-Transistor
  • the Ising model that the Ising model computer 200 is trying to simulate is briefly described.
  • the Ising model is a statistical mechanics model for explaining the behavior of magnetic materials.
  • the Ising model is defined by a spin that takes a binary value of + 1 / -1 (or 0/1, up / down), an interaction coefficient indicating the interaction between the spins, and an external magnetic field coefficient for each spin. .
  • the Ising model can calculate the energy at that time from a given spin arrangement, interaction coefficient, and external magnetic field coefficient.
  • the energy function of the Ising model is generally expressed by the following equation. Where ⁇ i and ⁇ j are the values of the i-th and j-th spins, J i and j are the interaction coefficients between the i-th and j-th spins, h i is the external magnetic field coefficient for the i-th spin, ⁇ I, j> represents a combination of two adjacent sites, and ⁇ represents a spin arrangement.
  • Calculating the ground state of the Ising model is an optimization problem for obtaining an array of spins that minimizes the energy function of the Ising model.
  • FIG. 1 is a diagram illustrating an example of the configuration of an information processing apparatus 100 that includes an Ising model computer 200 (details will be described later) and controls the Ising model computer 200 to perform information processing.
  • an Ising model computer 200 (details will be described later) and controls the Ising model computer 200 to perform information processing.
  • an accelerator configured with Ising model computers 200A and 200B is mounted on a personal computer that is generally used at present.
  • Ising model computers 200A and 200B they are simply referred to as reference numeral 200.
  • the information processing apparatus 100 has a central processing unit (CPU) 110, a temporary storage device (RAM) 120, and a fixed storage device (HDD) 160, and these are connected by a system bus 130.
  • a network interface controller (NIC) 140 is connected via the system bus 130, and the information processing apparatus 100 can communicate with other information processing apparatuses via the inter-device network 190.
  • the input / output controller 150 is connected to the system bus 130, and one or more Ising model computers (two Ising model computers 200A and 200B in the example of FIG. 1) are connected to the system bus 130.
  • the input / output controller 150 converts the protocol on the system bus 130 according to the interface of the Ising model computer 200.
  • Software for example, the control program 121 operating on the CPU 110 of the information processing apparatus 100 generally reads / writes to a specific address (so-called MMIO, MemoryeMapped I / O), and an input / output controller
  • the Ising model computer 200 can be controlled via 150.
  • the information processing apparatus 100 is typically considered to be a personal computer that is generally used at present. Therefore, the information processing apparatus 100 includes a display device such as a display and an input device such as a keyboard, and can be operated interactively by the user.
  • An analog computer such as the Ising model computer 200 is an arithmetic accelerator (accelerator) for a personal computer, and is added as a kind of peripheral device.
  • a graphics accelerator also referred to as GPU
  • PCI Express registered trademark
  • the Ising model computer 200 is also implemented in a similar form.
  • the Ising model computer 200 can be configured as, for example, one or more semiconductor integrated circuit devices (Ising chips).
  • the problem data 161 is stored in the HDD 160.
  • the problem data 161 may be input directly from the input device.
  • the problem data 161 is an Ising model that can be solved by the Ising model computer, and is stored on the HDD 160 as a list of interaction coefficients and external magnetic field coefficients that define the Ising model, for example. That is, the problem data 161 lists the interaction coefficients J i and j and the external magnetic field coefficients h i shown in Equation 1.
  • the control program 121 is a program executed by the CPU 110 and is a program for controlling the Ising model computer 200 based on the problem data 161 and the characteristic data 162.
  • the execution result of the control program 121 is an answer to the problem defined by the problem data 161 generated by the Ising model computer 200.
  • the characteristic measurement program 122 is a program executed by the CPU 110, and is a program for measuring manufacturing variations of the Ising model computer 200.
  • the characteristic measurement program 122 records the measurement result as characteristic data 162.
  • FIG. 2 is an example of a configuration diagram of the Ising model computer 200 of the present embodiment.
  • the Ising model computer 200 includes a spin array 210, an I / O driver 220, an I / O address decoder 230, an interaction address decoder 240, a random number generator 250, and a bit adjuster 251.
  • the Ising model computer in the present embodiment will be described assuming that it is implemented as a CMOS (Complementary Metal-Oxide-Semiconductor) integrated circuit that is widely used at present, but it can also be realized by other solid-state devices.
  • the spin array 210 is constituted by a semiconductor memory array. There is no limitation on the type of memory, but here, an example using SRAM (Static Random Access Memory) capable of high speed operation will be described.
  • the Ising model computer 200 can be driven by an SRAM compatible interface 222 for reading / writing the spin array 210.
  • the SRAM compatible interface includes an address bus 290, a data bus 291, an R / W control line 293, and an I / O clock 292.
  • the Ising model computer 200 may include another memory 270 such as a nonvolatile memory for control and the like.
  • the spin ⁇ i , the interaction coefficient J i , j and the external magnetic field coefficient h i of the Ising model are all expressed by information stored in the memory cells in the spin array 210. Reading / writing of these values to the memory cells of the spin array is performed by the I / O address decoder 230 and the I / O driver 220 via the SRAM compatible interface 222.
  • an address is given to the memory cell in the spin array 210. Further, the spin ⁇ i , the interaction coefficient J i , j , and the external magnetic field coefficient h i of the spin array 210 are associated with the address of the memory cell.
  • a technique for reading / writing information to / from a memory cell at an arbitrary address is the same as a known semiconductor memory control technique.
  • the address bus 290, the data bus 291 and the R / W control line 293 constituting the SRAM compatible interface operate in synchronization with a clock input to the I / O clock 292.
  • the interface does not have to be synchronous, and may be an asynchronous interface.
  • the Ising model computer 200 realizes an interaction between spins in the spin array 210 in order to perform a ground state search.
  • the interaction is to read the value of the memory cell of the spin array 210, calculate the local energy in the vicinity of the spin corresponding to the memory cell from the read value, and calculate the energy of the spin that decreases the energy.
  • This interaction is controlled from the outside by an interaction control interface. Specifically, an address for specifying a spin group to perform the interaction is input as an interaction address 280, and the clock input to the interaction clock 281 is used as a clock. Synchronize and interact. Internally, the interaction address decoder 240 generates a clock for updating a group of spins corresponding to the interaction address 280 in the spin array 210 from the interaction address 280 and the interaction clock 281.
  • the interaction is not necessarily realized by a clock synchronous circuit, and may be an asynchronous circuit.
  • the role of the interaction clock 281 is not an input of the clock but an enable signal that permits execution of the interaction.
  • the ground state search of the applied Ising model can be realized by energy minimization due to the interaction between spins, which will be described later, but this alone may lead to a local optimal solution. Since the interaction circuit 310 described below basically has only a movement in the direction of reducing the local energy, once it falls into the local optimum solution, it cannot escape from it and does not reach the global optimum solution. Therefore, in order to escape from the local optimal solution, a process that allows a transition that increases energy with a predetermined probability is introduced.
  • the spin value is inverted with a probability determined by the virtual temperature T.
  • the temperature T is obtained by comparing the ground state search process with physical annealing. At the initial stage of the ground state search, the temperature is set to a high temperature, a local search is performed while gradually decreasing the temperature, and cooling is finally performed until the temperature becomes zero.
  • a random number generator 250 and a random number generation clock 260 are used to invert the spin value with a predetermined probability.
  • the random number generator 250 generates a random number every clock of the random number generation clock 260.
  • the bit adjuster 251 adjusts the output bit from the random number generator 250 so as to invert the spin value with a high probability at the initial stage of the ground state search and to invert the spin value with a low probability at the final stage. .
  • 1 is large at the beginning of the ground state search, and 0 is at the end of the ground state search. Adjust the output to generate more.
  • the spin array 210 in the Ising model computer 200 corresponds to the Ising model and simulates it.
  • the spin array 210 includes, as a basic structural unit, a spin unit 300 (FIG. 3) that realizes the retention of one spin, the associated interaction coefficient, and the external magnetic field coefficient, and ground state search processing.
  • FIG. 3 is a block diagram showing an example of the configuration of unit components of the spin array 210.
  • FIG. 4 is an example of the Ising model, and the unit component in FIG. 3 simulates one spin unit.
  • FIG. 4 By arranging a plurality of spin units 300 shown in FIG. 3, an Ising model having a three-dimensional lattice topology as shown in FIG. 4 is formed.
  • the example of FIG. 4 is a three-dimensional lattice (18 spins in total) having a size of 3 (X-axis direction) ⁇ 3 (Y-axis direction) ⁇ 2 (Z-axis direction).
  • the definition of the coordinate axis is the X-axis in the right direction of the drawing, the Y-axis in the downward direction of the drawing, and the Z-axis in the depth direction of the drawing (the same applies to the following drawings).
  • This coordinate axis is only necessary for the convenience of description of the embodiment, and is not related to the invention.
  • a topology other than a three-dimensional lattice for example, a tree-like topology, it is expressed by the number of stages of the tree separately from the coordinate axes.
  • One spin unit 300 shown in FIG. 3 includes adjacent spins (for example, when there are five adjacent spins on the upper side, the lower side, the left side, the right side, and the depth side) ⁇ j , ⁇ k , ⁇ l , ⁇ m , ⁇
  • the value of n is input.
  • the spin unit 300 has J j, i , J k, i , J l, i , J m, i , which are the interaction coefficients with the adjacent spins. It has a memory cell that holds J n, i (an interaction coefficient with 5 adjacent spins). In the example of FIG. 3, the memory cell holding the interaction coefficient is placed on the input side, but can be placed on the output side.
  • the Ising model generally has an interaction expressed by an undirected graph.
  • J i, j ⁇ ⁇ i ⁇ ⁇ j exists as a term representing the interaction, and this indicates the interaction from the i-th spin to the j-th spin.
  • the general Ising model does not distinguish between the interaction from the i-th spin to the j-th spin and the interaction from the j-th spin to the i-th spin. That is, J i, j and J j, i are the same.
  • J i, j and J j, i are the same.
  • the Ising model computer 200 of this embodiment is extended to a directed graph, and the interaction from the i-th spin to the j-th spin and the interaction from the j-th spin to the i-th spin are made asymmetric. Is shown.
  • the interaction coefficients held by this spin unit are J j, i , J k, i , J l, i , J m, i 1 , J n, i are mutual reciprocals from adjacent j-th, k-th, l-th, m-th and n-th spins ⁇ j , ⁇ k , ⁇ l , ⁇ m , ⁇ n to i-th spin ⁇ i It determines the action.
  • an arrow (interaction) corresponding to an interaction coefficient included in the spin unit 300 changes from a spin outside the spin unit 300 shown in the figure to a spin inside the spin unit 300.
  • an arrow (interaction) corresponding to an interaction coefficient included in the spin unit 300 changes from a spin outside the spin unit 300 shown in the figure to a spin inside the spin unit 300.
  • the spin unit 300 holds the Ising model spin ⁇ i , interaction coefficients J j, i , J k, i , J l, i , J m, i , J n, i and external magnetic field coefficient h i.
  • a plurality of 1-bit memory cells are provided. This 1-bit memory cell is shown in FIG. 3 as SPIN holding spin ⁇ i , IS0 and IS1 holding external magnetic field coefficient h i , IU0, IU1, IL0, IL1, IR0, IR1, holding interaction coefficients. These are illustrated as ID0, ID1, IF0, and IF1.
  • IS0 and IS1, IU0 and IU1, IL0 and IL1, IR0 and IR1, ID0 and ID1, and IF0 and IF1 each play a role in one set, and therefore ISx, IUx, and ILx, respectively. , IRx, IDx, and IFx.
  • each memory cell included in the spin unit 300 a known SRAM memory cell can be used.
  • the memory cell has a data holding unit composed of two CMOS inverters. By controlling the gate of the pass gate transistor connected in series with the inverter with the word line and connecting the source / drain path to the bit line, data reading / writing to the data holding unit is realized.
  • the memory cell structure is not limited to this, and any structure that can store at least two values may be used.
  • other memories such as DRAM and flash memory can be used.
  • each memory cell is an SRAM memory cell that can store 1 bit.
  • multi-valued memory cells may be used.
  • the multi-valued memory cell one memory cell may be physically used, or a plurality of 1-bit memory cells may be used to store a plurality of bits.
  • the degree of freedom of the coefficient values that can be used increases.
  • the circuit scale and processing time increase. Therefore, an appropriate memory may be selected according to the type of problem to be solved and the application.
  • the spin unit 300 represents the i-th spin.
  • the memory cell SPIN is a memory cell for expressing the spin ⁇ i and holds the spin value.
  • the spin value is + 1 / ⁇ 1 (+1 is also expressed as “up” and ⁇ 1 is also expressed as “bottom”). In this specification, +1 corresponds to 1 and ⁇ 1 corresponds to 0, but the opposite may be possible.
  • FIG. 5 shows the correspondence between the memory cells of the spin unit 300 and the topology of the Ising model shown in FIG. ISx represents an external magnetic field coefficient.
  • IUx, ILx, IRx, IDx, and IFx each represent an interaction coefficient.
  • IUx is the upper spin (-1 in the Y axis direction)
  • ILx is the left spin (-1 in the X axis direction)
  • IRx is the right spin (+1 in the X axis direction)
  • IDx is the lower spin (Y axis).
  • IFx indicates an interaction coefficient with the spin (+1 or ⁇ 1 in the Z-axis direction) connected in the depth direction.
  • the Ising model when viewed from a certain spin, it has a coefficient of the influence of other spins on the own spin.
  • the coefficient of influence of the own spin on other spins belongs to each other spin. That is, the spin unit 300 is connected to a maximum of 5 spins.
  • One spin unit determines the spin value by interaction with five spins and the influence of a total of six external magnetic fields.
  • each of the memory cells ISx, IUx, ILx, IRx, IDx, and IFx is a combination of two memory cells whose end numbers are 0 and 1 (for example, IS0 and IS1 in the case of ISx). , + 1/0 / ⁇ 1.
  • + 1 / -1 is expressed by IS1
  • +1 is expressed when IS1 is 1
  • -1 is expressed when 0.
  • the spin unit 300 is driven, controlled, and read by the I / O address decoder 230 and the I / O driver 220, so that the memory cells in the spin unit 300 are compatible with the SRAM of the Ising model computer 200, as in a general SRAM. It can be read / written by the interface.
  • next state of the spin is determined so as to minimize the energy between the adjacent spins, which is a positive value when the product of the adjacent spins and the interaction coefficient and the external magnetic field coefficient are evaluated. Equivalent to judging which negative value is dominant.
  • the i-th spin sigma i as a spin ⁇ j, ⁇ k, ⁇ l , ⁇ m, the sigma n are adjacent, next state of the spin sigma i is determined as follows.
  • the external magnetic field coefficient may always be read as an interaction coefficient with a spin having a value of +1.
  • the local energy between the i-th spin and the adjacent spin is obtained by multiplying the above-described coefficient by the value of the i-th spin and further inverting the sign.
  • the local energy with respect to the j-th spin is -1 when the i-th spin is +1, and +1 when the i-th spin is -1. Therefore, the i-th spin is +1. It works in the direction to reduce the local energy here.
  • the interaction circuit 310 illustrated in the spin unit 300 of FIG. 3 shows an outline when a circuit for performing the above-described operation is expressed as a logic circuit.
  • XNOR exclusive OR
  • the next state of the spin that minimizes the energy when only the interaction is seen is (+1 is encoded as 1 and -1 is encoded as 0).
  • the interaction coefficient is only + 1 / ⁇ 1
  • the next state of the spin can be determined by determining which of the outputs is + 1 / ⁇ 1 by majority logic.
  • the value of the external magnetic field coefficient is simply a value to be input to the majority logic circuit 311 that determines the next state of the spin.
  • the output of the majority logic circuit 311 is input to the memory cell SPIN as the next state.
  • the random pulse train that is output in the bit adjustment period 251 is to randomly change the next state input to the memory cell SPIN according to the condition, and to prevent falling into a local solution during the calculation. Is.
  • the coefficient 0 is realized by not voting + 1 / ⁇ 1 as described above in the majority logic (represented by a switch in the interaction circuit 310). If the interaction circuit 310 is configured as a logic circuit, it is necessary to configure a majority logic circuit having a disable terminal corresponding to this switch corresponding to each input.
  • a circuit having an analog operation is a circuit that handles an electric signal that changes continuously.
  • the force that pulls the spin upward (+1) and the force that pulls the spin downward ( ⁇ 1) are expressed by the amount of current drawn, and the difference is amplified to determine the spin direction. That is, it can be said that the magnitude of the interaction between the spins is expressed in a similar manner by the flow rate of the current on the electronic circuit.
  • FIG. 6 is a circuit constituting the interaction circuit 310 based on an analog operation, in which a force that pulls the spin toward the +1 side is expressed as a current draw amount (FIG. 6A, CT side circuit), and ⁇ 1.
  • 6 is a circuit (referred to as a CB side circuit in FIG. 6B) that expresses the pulling force to the side as the amount of current drawn.
  • FIG. 7 is a diagram for explaining an example of the configuration of a circuit used in combination with the circuit of FIG. 6 in the interaction circuit 310 based on an analog operation.
  • FIG. 8 is a time chart for explaining an example of the operation realized by the circuits of FIGS.
  • the two circuits shown in FIG. 6 and the circuit of FIG. 7 that amplifies the output from these two circuits are operated according to the timing chart shown in FIG. 8, thereby realizing a function corresponding to the interaction circuit 310. I can do it.
  • a signal for supplying a value stored in each memory cell in FIG. 3 is indicated by a signal name with T added to the end of the name of the memory cell (for example, IS0).
  • the value 0/1 stored in is referred to as IS0T).
  • a signal obtained by inverting the signal is indicated by a signal name obtained by adding B to the end of the name of the memory cell (for example, a signal obtained by inverting IS0T is referred to as IS0B).
  • adjacent spin values are expressed as NUT (upper side), NLT (left side), NRT (right side), and NDT (lower side).
  • NFT there is an NFT in the depth direction, which is omitted due to the paper width of the drawing.
  • the circuit diagram of FIG. 6 it can be seen that the same configuration is repeated in the circuits related to NUT to NDT, but the portion related to NFT can be configured in the same manner.
  • the inverted value is expressed as NUB, NLB, NRB, NDB (there is NFB corresponding to NFT, but it is omitted in the figure).
  • CT is a line which makes a pair with CB which will be described later.
  • CT and CB are P-type MOSFETs driven by a clock signal CLKI and connected to the power supply voltage. When CLKI is LOW, the P-type MOSFET is turned on and CT and CB are precharged.
  • CT and CB are grounded to the ground by the N-type MOSFET driven by the clock signal CLKI at the bottom stage after passing through several N-type MOSFETs. That is, CT and CB are precharged when CLKI is LOW, and when CLKI is HIGH, the precharged charges are extracted to the ground through some of the N-type MOSFETs described above.
  • the number of N-type MOSFET groups in the path from the power supply voltage to the ground represents the force that pulls the spin to the +1 side in CT, and the spin in the CB represents ⁇ 1. Express the pulling force to the side.
  • the explanation will focus on the CT side that pulls the spin to the +1 side.
  • the external magnetic field coefficient is +1
  • the interaction coefficient regarding the interaction with the upper adjacent spin, the upper adjacent spin is +1 (NUT is 1) and the interaction coefficient is +1 (IU1T is 1) or the upper adjacent spin.
  • the basic idea of the circuit is the same as that on the CT side, except for the conditions for current flow. That is, the signal input to the gate of the N-type MOSFET is different.
  • FIG. 8 illustrates an operation realized by the circuits of FIGS. 6 and 7 as a timing chart.
  • the timing chart of FIG. 8 shows the voltage of the main signal lines of FIG. 6 and FIG. 7 in the interaction operation (the period during which this operation is performed is called an interaction period) performed by the circuits of FIG. 6 and FIG.
  • the time is shown on the axis and the time is shown on the horizontal axis.
  • the interaction period is divided into two phases, ie, a current period and a determination period.
  • CT and CB are precharged, and at the beginning of the current period, both signal levels are HIGH.
  • CLKI becomes HIGH
  • precharge is completed, and current flows from CT and CB to the ground through the N-type MOSFET group.
  • a voltage drop occurs in CT and CB in accordance with the force pulling to the +1 side (CT side) and the force pulling to the ⁇ 1 side (CB side).
  • CT side +1 side
  • CB side the force pulling to the ⁇ 1 side
  • the configuration examples of the circuits in FIGS. 6 to 8 are examples. For example, not limited to FIGS. 7 to 8, various known circuit configurations for measuring the voltage difference between CT and CB can be employed.
  • this state is referred to as an antagonistic state.
  • the antagonistic state there is no change in energy regardless of whether the spin direction is +1 or -1.
  • the spin direction is often randomly selected in the antagonistic state.
  • the previous spin value is often held in the sense of maintaining the current state.
  • the degree of voltage drop between CT and CB is different due to the above manufacturing variation problem, and one of them always has a lower voltage, that is, it always wins the majority. . Therefore, it can be said that the value of spin in the antagonistic state (hereinafter referred to as the antagonistic value) is determined by the manufacturing variation.
  • the problem of the antagonistic state is not limited to the examples shown in FIGS. 6 to 8, but can be said to be an essential problem for analog computers in particular. This is because even if the manufacturing variation as described above is suppressed as much as possible, in the antagonistic state, a minute variation appears as a phenomenon that the spin value is fixedly determined as either +1 or ⁇ 1 in the output result of the sense amplifier. is there.
  • Non-Patent Document 1 As manufacturing variation factors of transistors (mainly MOSFETs are used here), mainly dimensional errors in the manufacturing process, RDF (Random Doppant Fluctuation, see Non-Patent Document 1 for details, for example), and RTN (Random Brass Noise), For details, see Non-Patent Document 2, for example.
  • the dimensional error in the manufacturing process is a problem that the characteristics of the transistor differ due to an error in forming a geometric figure constituting a circuit on a semiconductor substrate by a technique such as photolithography. Error.
  • RDF is a random variation in the position and density of impurity atoms that occurs when doping impurities (Donor or Acceptor, generically Doppant) during semiconductor manufacturing, and this causes variations in the threshold voltage of the transistor. Variations due to RDF are determined at the time of manufacture, and then hardly change during the lifetime of the transistor. In addition, when the size of the transistor itself is large, the effect appears to be relatively small due to random variations in the position and density of the impurity atoms, and the effect appears to be large as miniaturization progresses. come. Among the miniaturization processes, especially for a transistor that needs to keep the threshold voltage variation low, measures such as increasing the gate length are taken to reduce the influence of RDF by increasing the gate area.
  • RTN is a phenomenon in which variations in threshold voltage occur dynamically when carriers are trapped and released by an insulating film (such as an oxide film) at the gate of a transistor.
  • an insulating film such as an oxide film
  • RTN phenomena with various time constants are combined and observed as threshold voltage variations.
  • the time constant is considered to be determined by some factor at the time of manufacture. However, since the time constant is extremely short to long, it becomes a factor that makes screening during manufacture difficult. Differences from RDF include the fact that RTN has a stronger degree of influence with miniaturization than RDF and that the threshold voltage is dynamically changed.
  • the Ising model computer 300 shown in the present embodiment is manufactured by a process of a generation of about 65 nm, a static influence that is estimated to be caused by a dimensional error or RDF.
  • the dynamic variation considered to be the influence of RTN has not yet been realized.
  • the present inventors measure the characteristics of each individual in advance, and modify the problem that is input to each individual based on the measurement result, thereby adapting to the characteristic variations of the elements that constitute each individual. Devised.
  • the characteristics of each individual need only be measured once, it is possible to measure at the time of shipment from the factory and to provide the individual corresponding to the information to the customer.
  • the handling of the characteristic variation in the present embodiment is realized by the control program 121 and the characteristic measurement program 122 of the information processing apparatus 100 that controls the Ising model computer 200.
  • FIG. 9 is a flowchart showing the operation of the control program 121.
  • the operation based on this flowchart is as follows.
  • information processing by the control program 121 that is software is realized by using hardware resources of the information processing apparatus 100 such as the CPU 110, the RAM 120, and the HDD 160.
  • the CPU 110 performs overall control and calculation, and also controls read / write to the RAM 120 and HDD 160. These are the same in other flowcharts described below.
  • Step 901 This flowchart is a procedure for solving the Ising model by the Ising model computer 200. Prior to the execution of this flowchart, it is assumed that the user has prepared the problem data 161 on the HDD 160 by inputting it via an input device.
  • the information processing apparatus 100 includes a plurality of Ising model computers 200A and 200B. However, the CPU 110 selects an Ising model computer that allows the user to solve the problem data 161, or the CPU 110 appropriately selects in step 901. Shall.
  • the Ising model computer 200 refers to the Ising model computer selected in step 901.
  • Step 902 It is checked whether or not the characteristic data 162 corresponding to the Ising model computer 200 included in the information processing apparatus 100 exists on the HDD 160. If not, the characteristic data 162 is generated in step 903. If it exists, step 903 is not executed and the processing transitions to step 904. Details of the characteristic data 162 will be described later.
  • Step 903 and Step 904 The characteristic measurement program 162 is executed to generate characteristic data 162 corresponding to the Ising model computers 200A and 200B. Details of the characteristic measurement program 162 will be described later.
  • Step 905 Based on the characteristic data 162, the problem data 161 is converted so as to avoid characteristic variation. Details of the conversion performed in this step will be described later.
  • the converted problem data is referred to as post-conversion problem data.
  • Step 906 The converted problem data is input to the Ising model computer 200. Specifically, the interaction coefficient and the external magnetic field coefficient constituting the post-conversion problem data are written into the memory cell group constituting the Ising model computer 200. This writing is realized by the input / output controller 150 controlling the SRAM compatible interface 222 of the Ising model computer 200. Therefore, in step 906, the CPU 110 controls the input / output controller 150 by means such as MMIO.
  • Step 907 The solution is instructed to the Ising model computer 200. Specifically, the interaction address 280, the interaction clock 281 and the random number generation clock 260 of the Ising model computer 200 are controlled. Since this control is also directly performed by the input / output controller 150, the CPU 110 controls the input / output controller 150 by means such as MMIO.
  • Step 908 The result obtained by the Ising model computer 200, that is, the spin array is read out.
  • the CPU 110 controls the input / output controller 150 by means such as MMIO.
  • Step 909 The spin array read in step 908 is converted to obtain a solution to the problem data 161 (referred to as a spin array after inverse conversion). Since the problem solved by the Ising model computer 200 is the converted problem data, the solution obtained in step 908 is the solution of the converted problem data. It is necessary to interpret this solution as a solution of the problem data 161 before conversion.
  • the spin array obtained by this interpretation is a spin array after inverse transformation. Details of the interpretation method will be described later.
  • Step 910 Present the spin array after inverse transformation to the user as a solution result.
  • the solution result may be stored in the HDD 160 of the information processing apparatus 100 or may be output to the outside via the output device as it is.
  • FIG. 10 shows an example of the configuration of the characteristic data 161.
  • the characteristic data 161 is a list showing the antagonistic value of each spin in this embodiment. Since the antagonistic value is different for each individual Ising model computer, the characteristic data 161 has an entry for each individual (individual information entries 1010 and 1020). For example, the characteristics of the Ising model computer 200A are stored in the individual information entry 1010, and the characteristics of the Ising model computer 200B are stored in the individual information entry 1020.
  • Each entry has a configuration such as a serial number portion 1011, a time stamp 1012, and an antagonistic value recording portion 1013, for example, in the example of the individual information entry 1010. Since the individual information entry 1010 corresponds to an individual of the Ising model computer, information for identifying the individual of the Ising model computer, that is, a serial number is included in the individual information entry, thereby associating the individual information entry 1010 with the individual. Further, according to the study by the present inventors, the characteristic of the antagonistic value is almost determined at the time of manufacture, but it cannot be said that there is no secular change. Therefore, the information constituting the individual information entry 1010 in the time stamp 1012 Stores the date and time of the measurement. The antagonistic value is recorded in the antagonistic value recording unit 1013 with the antagonistic value of each spin of the Ising model computer 200 together with the coordinates of each spin (in this example, the X coordinate, the Y coordinate, and the Z coordinate).
  • FIG. 11 shows the operation of the characteristic measurement program 162 for generating the characteristic data 161 as a flowchart. The operation based on this flowchart is as follows.
  • Step 1101 This flowchart is a procedure for measuring the antagonistic value of each spin of the Ising model computer 200.
  • the information processing apparatus 100 includes a plurality of Ising model computers 200A and 200B, but this procedure is executed for each Ising model computer (may be executed in parallel).
  • Step 1102 Create a problem (measurement model) for measuring the antagonistic value.
  • An example of a model for causing an antagonistic state is shown in FIG.
  • the model of FIG. 12 has +1 and ⁇ 1 interactions with spins at both ends (Y-axis direction + 1 and Z-axis direction + 1), respectively, in order to know the antagonistic value of a certain spin, and the spin values at both ends. Are each set to +1 to generate an antagonistic state.
  • this model when an interaction is performed once at a temperature of 0 (a state in which random spin inversion is not performed), the spin of interest shifts to +1 or ⁇ 1 depending on manufacturing variations.
  • the antagonistic value of one spin for example, as shown in FIG. 12, when one spin adjacent to each other in the Y direction and the Z direction is used, the other spins excluding those spins are By including the same interaction in one measurement model, it can be measured simultaneously.
  • Step 1103 The measurement model created in step 1102 is input to the Ising model computer 200.
  • the input method is the same as in step 906.
  • Step 1104 The Ising model computer 200 is instructed to execute a solution so as to execute an interaction for each of the spins to be measured in the measurement model.
  • the spin to be measured is set in an antagonistic state by the measurement model.
  • the temperature is 0. That is, the random spin inversion is suppressed by stopping the random number generation clock 260 or changing the setting of the random number generator 250 or the bit adjuster 251.
  • Step 1105 Read the spin arrangement from the Ising model computer 200. The method is similar to step 908.
  • Step 1106 Of the spin arrangement read in Step 1105, the value of the spin to be measured in the measurement model generated in Step 1102 is recorded in the characteristic data 161 as the antagonistic value of the spin.
  • Step 1107 It is determined whether all spin antagonistic values have been measured. If measurement of all the spins is not completed, a measurement model for measuring the remaining spins is generated, and the process returns to step 1102 to perform measurement again. When the measurement of all spins is completed, the process proceeds to step 1108.
  • Step 1108 After completing the generation of the characteristic data 161, the characteristic measurement program 162 ends.
  • the characteristic data 161 is held in the RAM 120 or the HDD 160 so that the control program 121 can be used. Alternatively, it may be stored in a nonvolatile memory 270 such as a flash memory in the Ising model computer 200. If the Ising model computer 200 has its own characteristic data 161, the characteristic data 161 can be used even when the Ising model computer 200 is used in another information processing apparatus.
  • the Ising model (problem data) is converted based on the antagonistic value shown in the characteristic data 161 and the spin value (hereinafter referred to as desired value) that should be at the time of the antagonistic, and in combination with step 909, the spin at the time of the antagonistic Convert the model so that the value is the desired value.
  • the antagonistic value of the i-th spin is u i
  • the desired value of the i-th spin is e i .
  • the interaction coefficient J i, j is converted into J i, j ⁇ u i ⁇ u j ⁇ e i ⁇ e j .
  • the external magnetic field coefficient hi is converted to h i ⁇ u i ⁇ e i .
  • step 909 as a conversion paired with step 905, the spin value read in step 908 is converted into a solution of the model before conversion.
  • step 905 and step 909 the conversion performed in step 905 and step 909 is such that when the antagonistic value is determined by the manufacturing variation of the element, the value at the time of the antagonistic state can be converted into the desired value by post-processing in step 909.
  • step 905 it can be said that the sign of the coefficient is switched.
  • step 909 the desired value can be obtained by changing the sign of the spin that is the solution.
  • the desired value needs to be determined.
  • the desired value is +1 if the external magnetic field coefficient h i of each spin is positive, and the desired value is ⁇ 1 if it is negative. This is based on the idea of giving priority to the external magnetic field when an antagonistic state occurs as the sum of the interaction and the external magnetic field (or the idea that the external magnetic field is a minute amount larger than the interaction).
  • the value may be stored in the HDD 160 in association with the problem data 161, for example. Alternatively, it may be input separately from the problem data 161 and stored in the RAM 120 so that the control program 121 can be used.
  • the same problem is generally solved several times to obtain the best solution among them.
  • a desired value is determined at random, and an attempt is made to solve the problem multiple times.
  • random data may be stored in advance or may be generated each time.
  • the desired value is not determined in advance
  • the basic processing flow is the same as the example shown in FIG. 9, but the desired value is randomly changed during the processing 904 for converting the problem data.
  • the process 909 for converting the spin arrangement a conversion process corresponding to the problem data is performed.
  • the timing for changing the desired value is arbitrary, but if the change cycle is shortened, the effect becomes large, but the processing becomes complicated.
  • the desired value is randomly changed.
  • a method of changing the value alternately may be considered.
  • the Ising model computer 200 having a structure similar to the Ising model is shown as an example.
  • the present invention can be similarly applied to a neuro computer having a structure similar to a neural network.
  • the neural network receives the firing state of other neurons (for example, expressed by the numerical values 0 and 1), and multiplies the firing state of other neurons by the synapse connection coefficient between the neurons and whether or not the own neuron fires. The sum is calculated and the sum is passed through a threshold function and determined by its output.
  • this operation is implemented with an analog electronic circuit, it has a structure similar to the Ising model computer of this embodiment (if the ignition state is read as the spin value and the synapse coupling coefficient is interpreted as the interaction coefficient, the structure is almost the same as the Ising model. Become).
  • the configuration of the interaction circuit can basically use the configuration shown in FIG. 6, the firing state instead of the spin value NUT and the like, and the synapse instead of the interaction coefficient IU1T. What is necessary is just to input a coupling coefficient.
  • Example 2 considers the case where variations other than the antagonistic state are also considered in the Ising model computer 200.
  • the first embodiment an implementation corresponding to the problem that the spin value that appears when the interaction is antagonized falls due to variations in the manufacturing characteristics of the MOSFET in the interaction circuit 310 of the Ising model computer 200 falls in a fixed direction.
  • the method shown in the first embodiment is a method that corresponds to the control of the Ising model computer 200 and the information processing apparatus 100 including the same without adding an additional circuit for compensating for the interaction circuit 310. It was.
  • Example 2 discloses an example of a configuration that can cope with the case where the influence of manufacturing variation is further increased and the influence affects not only the antagonistic state but also other interactions.
  • the variation in the current (I DS ) between the drains of the MOSFET depends on the channel length (L) and the channel width (W), and is known to be proportional to 1 / ⁇ (W ⁇ L). . Therefore, in order to suppress the variation of 1 / n of the I DS must doubled n the (W / L). This indicates that the suppression of variation and miniaturization (high integration) are in a contradictory relationship.
  • the value ranges of the interaction coefficient and the external magnetic field coefficient were +1 to ⁇ 1.
  • the interaction and the external magnetic field are combined, and six types of force are applied to one spin (total of six types in the external magnetic field in addition to five types in the vertical and horizontal depth directions of the interaction).
  • the force directed to +1 and the force directed to -1 were half, the antagonistic state was reached, and the response to the manufacturing variation in this antagonistic state was Example 1.
  • the spin is directed in a direction where the force is larger.
  • the force directed to +1 and the force directed to -1 are not halved, it is naturally expected that the spin is directed in a direction where the force is larger.
  • the force is pulled by a force of 3 on the +1 side and 2 on the -1 side in the 5 directions of interaction.
  • the circuit since the current difference between the CT side circuit and the CB side circuit is 1/6, in this case, if a current of 1/6 or more is further passed as a result of the current variation of the MOSFET of the CB side circuit, the circuit does not operate correctly. Alternatively, there is a problem in that the amount of current flowing on the CT side is small due to variations and the amount on the CB side is large.
  • the interaction coefficient and the external magnetic field coefficient are +1 to ⁇ 1 and there are six types of forces including the interaction coefficient and the external magnetic field coefficient, the above problem does not occur.
  • the relationship between the aforementioned variation and the channel width / channel length shows how much area efficiency is reduced in a complicated design.
  • the coefficient width is +3 to -3, three times as many transistors are required for each of the CT side circuit and the CB side circuit as compared with the current design.
  • an antagonistic state occurs, there are 9 to 9 cases in which the coefficient + 3 / -3 is used in each of three directions.
  • the difference is reduced to 1/3. That is, the characteristic variation of the MOSFET must be reduced to 1/3, and for that purpose, the size of the MOSFET must be increased 9 times.
  • the number of MOSFETs must be increased in order to suppress the variation in characteristics of the MOSFETs, and the number of MOSFETs is increased in order to realize originally complex interactions.
  • the number of transistors in addition to the size 9 times to suppress the characteristic variation, the number of transistors is originally 3 times to realize complex interaction, so a total area of 27 times is required. . For this reason, it becomes a factor which remarkably inhibits scalability.
  • this embodiment shows an example corresponding to the case of obvious variations in I DS Besides the antagonistic state. For example, even though a 3 to 2 cases not antagonistic state, the I DS variations in MOSEFT the CT side or CB side, so that the side draw at second force than a side draw in the third force would be misidentified as strong An example is shown to deal with various cases.
  • FIG. 13 is a configuration diagram of an analog computer 1300 according to the embodiment of the present embodiment.
  • This analog computer 1300 is assumed to be incorporated and used in a part of an apparatus such as the information processing apparatus 100, similarly to the Ising model computer 200A of the first embodiment.
  • FIG. 13 is described as an analog computer that abstracts the Ising model computer 200 of the first embodiment.
  • the analog computer 1300 includes analog calculation units 1310A to 1310A (referred to simply as 1310 when it is not necessary to distinguish between A and E), an analog arithmetic circuit 1311, a storage element 1312, an inspection circuit 1313, and an inter-unit interface 1320. , And an inspection interface 1330.
  • the Ising model computer 200 corresponds to the analog computer 1300
  • the spin unit 300 corresponds to the analog calculation unit 1310.
  • the interaction circuit 310 corresponds to, for example, the analog arithmetic circuit 1311 whose specific examples are shown in FIGS.
  • a memory cell array in the spin unit 300 corresponds to the storage element 1312
  • a transmission line that transmits the spin value of the adjacent spin unit between the spin units corresponds to the inter-unit interface 1320.
  • the inter-unit interface 1320 and the inspection interface 1330 have a bus-like topology in FIG. 13, but this topology is arbitrary. In the case of the Ising model computer 200 of the first embodiment, it can be said that the inter-unit interface has a three-dimensional lattice-like topology.
  • the inspection circuit 1313 detects a case where the processing result of the analog arithmetic circuit 1311 may be erroneous, and outputs the result to the inspection interface 1330. There are two cases where the processing result of the analog arithmetic circuit 1311 is incorrect. The first is the already described antagonistic state. Second, although it is not in an antagonistic state (that is, there is no error if there is no manufacturing variation of the analog arithmetic circuit 1311), the pulling force on the +1 side and the pulling force on the -1 side are close, so depending on the manufacturing variation This is a case where there is a possibility of erroneous determination. It is possible to estimate to some extent how much manufacturing variation will occur from process generation and MOSFET design parameters.
  • a combination of a simple inspection circuit 1313 incorporated in the analog computer 1300 and the overall control of the information processing apparatus 100 including the analog computer 1300 can suppress an increase in the amount of the analog computer 1300.
  • FIG. 14 shows an example of the configuration of the inspection circuit 1313.
  • the inspection circuit 1313 detects whether or not the voltage difference between the CT side circuit and the CB side circuit is within a certain range in order to detect the above-described two cases at once. If it falls within a certain range, it is determined that there is a possibility of malfunction. As a specific example, when the voltage difference between the CT side circuit and the CB side circuit is smaller than a predetermined threshold, it is determined that there is a possibility of malfunction.
  • the range in which the voltage difference is allowed can be set by estimating how much manufacturing variation occurs.
  • the sense amplifier 14 has a basic structure of the latch type sense amplifier of FIG.
  • the sense amplifier functions to amplify the voltage difference between the two types of signals (CT and CB).
  • CT and CB the driving force of the MOSFET 1401 shown in FIG. 14
  • the voltage difference is gradually amplified and the gradient has a small voltage difference. It becomes moderate. By detecting this gentle gradient, it is detected whether the voltage difference is within an allowable range.
  • FIG. 15 is a time chart when the inspection circuit 1313 is operated.
  • the structure and operation of the inspection circuit 1313 will be described in detail with reference to FIG. Similar to the time chart of FIG. 8, current is drawn from the pre-charged CT / CB by the MOSFET group whose ON / OFF is determined by the adjacent spin value, interaction coefficient, and external magnetic field coefficient, causing a voltage drop in CT / CB. Let After that, the sense amplifier is operated by setting CLKIB to HIGH. However, since the driving force of the MOSFET 1401 is suppressed as described above, the slope of the voltage appearing at CLT and CLB becomes gentle.
  • the inspection circuit 1313 is realized by detecting that the slope is less than a constant.
  • a buffer having a threshold value and an exclusive OR (XOR) circuit are used in this embodiment.
  • the output OUTA from the XOR is observed at a timing when CLKIB is delayed by a minute time ⁇ t.
  • the possibility of malfunction can be inspected depending on whether the observed level of OUTA is HIGH or LOW.
  • an edge trigger type flip-flop (FF) may be driven with CLKIB delayed by ⁇ t.
  • the structure of the inspection circuit 1313 is similar to the sense amplifier of FIG. 7, it can also be used as a sense amplifier.
  • FIG. 14 is merely an example, and a circuit for measuring the voltage difference between CT and CB may be provided completely separate from the sense amplifier. In this case, various known circuit configurations for measuring voltage can be employed.
  • FIG. 16 is a flowchart showing the operation of the information processing apparatus 100 corresponding to the malfunction of the analog computer 1300 using the test circuit 1313.
  • Step 1601 In this flow, the processing performed in Step 907 in the control of the Ising model computer 200 (control of the analog computer 1300 in this embodiment) by the information processing apparatus 100 shown in FIG. A detailed flow of control using the circuit 1313 is shown.
  • Step 1602 Each of the analog calculation units 1310A to 1310B included in the analog computer 1300 is operated once to execute the calculation.
  • Step 1603 As a result of the calculation performed in Step 1602, the inspection circuit 1313 checks whether the possibility of malfunction of the analog arithmetic unit 1311 is detected. If the possibility of malfunction is not detected, the process proceeds to step 1607. If the possibility of malfunction is detected, the process proceeds to step 1604.
  • Step 1604 If there is a malfunction, the information obtained from the test circuit 1313 via the test interface 1330 is used to identify the analog calculation unit 1310 that may have malfunctioned, and the unit and the memory elements around the unit (For example, in the case of the Ising model computer 200 of the first embodiment, only the spin value of the adjacent spin needs to be read).
  • Step 1605 The calculation that should be performed by the analog calculation unit 1310 that has malfunctioned is performed on the CPU 110 of the information processing apparatus 100 by using the information obtained in Step 1604.
  • the CPU of the information processing apparatus 100 may be used as the alternative CPU, but a dedicated CPU or a logic circuit may be provided in or in association with the analog computer 1300.
  • a dedicated CPU or a logic circuit may be provided in or in association with the analog computer 1300.
  • the dedicated CPU may be programmable by software, or can be realized by hardware such as FPGA (Field Programmable Gate Array) and ASIC (Application Specific Integrated Circuit). Moreover, you may comprise with exclusive hardware.
  • a circuit corresponding to the interaction circuit 310 shown in FIG. 3 can be used as a digital arithmetic circuit.
  • the analog arithmetic circuit 1311 may be replaced with this digital arithmetic circuit.
  • the digital arithmetic circuit can be provided outside or inside the analog computer 1300. If one digital arithmetic circuit is provided for a plurality of analog arithmetic circuits 1311, an effect of reducing the circuit scale can be obtained.
  • Step 1606 The calculation result of step 1605 is written back to the storage device 1312 of the corresponding analog calculation unit 1310.
  • Step 1607 It is confirmed whether a predetermined number of operations to be performed in this flow is satisfied. If not, return to Step 1602 and repeat.
  • Step 1608 When the predetermined number of operations is satisfied, this flow is finished.
  • the CPU 110 of the information processing apparatus 100 temporarily plays the role of the analog arithmetic circuit 1311 only at that time.
  • the information processing apparatus 100 as a whole can perform processing without malfunction.
  • the CPU 110 serving as a substitute for the analog computer unit 1310 may be provided in the analog computer 1300.
  • FIG. 17 is a diagram illustrating another example of the configuration of the inspection circuit 1313.
  • the inspection circuit 1313 may be mounted like an inspection circuit 1700 illustrated in FIG.
  • an OUTA signal (output from the XOR) is taken into a flip-flop with clock signals with different timings generated by two types of delay circuits (first delay circuit and second delay circuit).
  • first delay circuit generates a clock signal whose clock rises at time t
  • second delay circuit generates a clock signal whose clock rises at time t + ⁇ t.
  • the magnitude of the voltage difference can be detected by appropriately setting the time t of the first delay circuit and the time t + ⁇ t of the second delay circuit.
  • the inspection circuit 1313 detects all when the voltage difference is small, but the inspection circuit 1700 determines that there is no problem when the voltage difference is smaller than the predetermined value and larger than the predetermined value. If there is a potential difference, it can be detected that there is a problem.
  • the embodiment is shown by way of an example of the Ising model computer 200 having a similar structure to the Ising model.
  • the present invention is similarly applied to a neuro computer having a similar structure to the neural network. This is possible as in the first embodiment.
  • the present invention is not limited to the above-described embodiment, and includes various modifications.
  • a part of the configuration of one embodiment can be replaced with the configuration of another embodiment, and the configuration of another embodiment can be added to the configuration of one embodiment.
  • Information processing apparatus 200 Ising model computer 300: Spin unit 1300: Analog computer

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Logic Circuits (AREA)

Abstract

　素子の製造ばらつきの影響を回避して情報処理を実現する手段を提供することを目的とすることを目的とする。　１個以上のアナログ計算機を含む情報処理装置において、アナログ計算機の特性ばらつきを計測して記録するステップを有することを特徴とする情報処理装置が開示される。

Description

情報処理方法及び装置
　本発明は情報処理装置において、求解すべき問題をそれに対応した構造の装置として実現したアナログ計算機に関するものである。
　現在主流の計算機は、トランジスタをスイッチ素子として捉えて構成したデジタル計算機である。一方で、トランジスタなどの素子が持つ特性をより直接的に利用して情報処理を実現するアナログ計算機がある。アナログ計算機の代表的なものに演算増幅器を用いて微分方程式を解くものがかつて盛んに用いられた（例えば、特許文献１）。しかし、アナログ電子回路の実装に起因する精度の制限から、デジタル計算機上での数値計算法に取って替わられた（例えば、特許文献２）。
　ところで、近年、ニューラルネットワークに代表される創発型の計算モデルが広く用いられている。創発とは部分の性質の単純な総和にとどまらない性質が全体として現れることである。ニューラルネットワークであればニューロン間の局所的なシナプス結合によって成される振舞いが、全体としてはある種の情報処理を実現している（例えば、特許文献３）。創発型の計算モデルとして、例えば他に遺伝的アルゴリズムや群知能があげられる。また、局所的な相互作用に基づく計算モデルであれば創発型と言えるものは多い。創発型の計算モデルの特徴として、一般的にそれを実現する個々の操作の精度が多少低くても、大局的には影響を及ぼさないことがあげられる。この特徴に基づいて、前述したアナログ計算機が創発型の計算モデルを実現する上で再度注目を浴びている（例えば、特許文献４）。また、特許文献５はイジングモデルと呼ばれる統計物理学で磁性体の振舞いを説明するモデルにおいて、エネルギー最小の状態のスピン配列を求めることが一種の最適化問題として捉えられることを用いて、イジングモデルと相似な構造のＣＭＯＳ回路を半導体基板上の形成したアナログ計算機の一種である。
　アナログ計算機は計算モデルの構成要素と対応する素子を多数個空間上（例えば半導体素子として実現する場合には半導体基板上）に並べて構成する。そのため、素子の製造時のばらつきが問題となる。ニューラルネットワークを実現する場合、同じニューロンを模擬する素子であったとしても、製造時のばらつきによって素子毎に特性が異なる。アナログ計算機ではその振舞いが素子の入出力特性の個体差に影響を受ける。そこで、素子の個体差を補償するような回路を付加することが行われている（例えば、特許文献６）。
特開昭６１－２８５５８３号 特開平２－２２６１８６号 特開２００２－８０３２号 特開２００７－６６２５８号 ＰＣＴ／ＪＰ２０１４／０７３２２６ 特開平１０－１１６３０９号
Stolk， P． A．， Widdershoven， F． P． & Klaassen， D． B． M． Modeling statistical dopant fluctuations in MOS transistors． IEEE Trans． on Electron Devices 45， 1960－1971 (1998)． Miki， H． et． al． Understanding short－term BTI behavior through comprehensive observation of gate－voltage dependence of RTN in highly scaled high－κ/metal－gate pFETs． In Symp． on VLSI Tech． 148－149 (IEEE， 2011)．
　本発明では、アナログ計算機を備える情報処理装置において、アナログ計算機を構成する素子の特性ばらつきの影響を回避することのできる情報処理装置を提供することを目的とする。
　上記課題を解決するために、代表的な本発明の情報処理装置の一つは、１個以上のアナログ計算機を含む情報処理装置において、アナログ計算機の特性ばらつきを計測して記録するステップを実行する機能を有することを特徴とする情報処理装置である。
　本発明の他の一側面は、処理の対象もしくは処理の結果となる数値を格納する数値格納メモリと、処理に用いる係数を格納する係数格納メモリと、数値と係数を用いて処理の少なくとも一部をアナログ的に実行するアナログ処理回路と、を複数組備える情報処理装置である。この情報処理装置は、複数のアナログ処理回路の特性ばらつきを示す特性データを格納する特性データ格納メモリを備える。この情報処理装置は、問題データを数値または係数として数値格納メモリと係数格納メモリに格納し、アナログ処理回路による１または複数回の処理の後、数値格納メモリに格納された数値を解データとして読み出すように動作する。
　本発明の好ましい態様では、問題データを格納する前に、特性データに基づいて問題データを変換し、解データを読み出した後に、特性データに基づいて解データを変換する。
　本発明の具体的な構成の一例としては、第１の数値格納メモリと、第１の係数格納メモリと、第１のアナログ処理回路と、第２の数値格納メモリと、第２の係数格納メモリと、第２のアナログ処理回路と、を有する計算ユニットを備る。第２のアナログ処理回路が、第１の数値格納メモリに格納された数値と、第２の係数格納メモリに格納された係数を用いて処理を行い、第２の数値格納メモリに処理の結果を格納する。
　本発明の具体的な構成の他の一例とは、処理装置、記憶装置、入力装置、出力装置を備える。入力装置が問題データを受付け、処理装置が受付けた問題データを記憶装置に格納する。処理装置が記憶装置に格納されている問題データを読み出し、特性データに基づいて読み出した問題データを変換し、計算ユニットの数値格納メモリと係数格納メモリに格納する。処理装置は、計算ユニットのアナログ処理回路による複数回の処理の後、数値格納メモリから解データ読み出し、特性データに基づいて読み出した解データを変換し、記憶装置または出力装置に出力する。
　本発明の具体的な構成の他の一例は、処理装置の制御により、複数のアナログ処理回路の特性ばらつきを測定し、測定結果を特性データとして特性データ格納メモリに格納する。特性データ格納メモリは、上記の記憶装置の一部を利用してもよい。また、専用のメモリを設けてもよい。特性データ格納メモリの装置内の位置は限定されない。
　本発明の具体的な構成の他の一例は、複数のアナログ処理回路の処理内容を検査する検査回路と、アナログ処理回路による処理と同等の処理をデジタル的に実行するデジタル処理回路と、を有し、検査回路の検査結果に基づいて、アナログ処理回路の処理をデジタル処理回路に代替させる。デジタル処理回路は、上記の処理装置の一部を利用してもよい。また、専用の論理回路を設けてもよい。デジタル処理回路の装置内の位置は限定されない。
　検査回路による検査の例としては、アナログ処理回路の処理において生成される２つのアナログ値を検査して、アナログ値の差が第１の閾値より小さく、第２の閾値より大きい場合に、前記アナログ処理回路の処理を前記デジタル処理回路に代替させるものがある。
　本発明の具体的な適用例としては、問題データは、イジングモデルのスピン配列を求めるための問題データ、あるいは、ニューラルネットワークのニューロンの状態を求めるための問題データである。
　本発明の他の一側面は、処理の対象もしくは処理の結果となる数値を格納する数値格納メモリと、処理に用いる係数を格納する係数格納メモリと、数値と係数を用いて少なくとも一部の処理をアナログ的に実行するアナログ処理回路と、を複数組備える情報処理装置である。この情報処理装置は、記複数のアナログ処理回路の処理内容を検査する検査回路と、アナログ処理回路による処理と同等の処理をデジタル的に実行するデジタル処理回路と、を有し、検査回路の検査結果に基づいて、アナログ処理回路の処理をデジタル処理回路に代替させる。
　本発明の具体的な一例を示すと、検査回路による検査は、アナログ処理回路の処理において生成される２つのアナログ値（例えば電圧値）の差を検査するものであり、アナログ値の差が所定の条件下にある場合、アナログ処理回路の処理をデジタル処理回路に代替させる。デジタル処理回路としては、汎用のマイクロプロセッサを用いてもよいし、専用の論理回路を用いてもよい。
　本発明の他の一側面は、１個以上のアナログ計算機を含む情報処理装置を用いた情報処理方法において、アナログ計算機の特性ばらつきを計測して特性ばらつきデータとして記録するステップ、を有することを特徴とする情報処理方法である。
　本発明のさらに具体的な構成例では、特性ばらつきデータを参照して、アナログ計算機に入力するデータを変換するステップと、特性ばらつきデータを参照して、アナログ計算機から出力されたデータを解釈するステップとを有する。
　本発明のさらに具体的な構成例では、アナログ計算機はイジングモデルのスピン配列を求めるものであり、入力するデータの変換は、イジングモデルの相互作用係数乃至は外部磁場係数の値の符号を反転させるものであり、出力されたデータの解釈は、スピン値の反転である。
　本発明のさらに具体的な構成例では、アナログ計算機は特性ばらつきの影響を評価する構成単位毎に、その評価を行う検査部を有しており、検査部からの検査結果を読み出すことで、誤動作の可能性のある構成単位を検出する。
　本発明のさらに具体的な構成例では、誤動作の可能性のある構成単位を検出したときに、情報処理装置が有する論理回路が当該構成単位の担う処理を一時的に代替する。
　本発明によれば、素子の製造ばらつきがあるアナログ計算機を利用して、回路物量の大幅な増加なく、その影響を回避した情報処理を実現することが出来る。
本発明を適用した情報処理装置の構成の一例を説明するブロック図である。 本発明によって製造ばらつきを回避して情報処理を実現するアナログ計算機の構成の一例を説明するブロック図である。 図２に示すアナログ計算機を構成する単位構成要素の構成の一例を説明するブロック図である。 図２に示すアナログ計算機が模擬する計算モデルの一例を説明するモデル図である。 図３に示す単位構成要素と、図４に示す計算モデルとの対応関係を説明するモデル図である。 図３の単位構成要素中で用いられている回路であり、製造ばらつきが生じる回路の構成の一例を説明する回路図である。 図６の回路と合わせて用いられる回路の構成の一例を説明する回路図である。 図６及び図７の回路によって実現される動作の一例を説明するタイムチャートである。 製造ばらつきを回避して情報処理を実現するために、図１の情報処理装置が行う動作の一例を説明する流れ図である。 図１の情報処理装置が保持する特性ばらつきを記録したデータの一例を説明する概念図である。 特性ばらつきの計測を実現するために、図１の情報処理装置が行う動作の一例を説明する流れ図である。 特性ばらつきの計測を実現するために、図２のアナログ計算機に入力されるデータの一例を説明するモデル図である。 本発明を適用したアナログ計算機の構成の一例を説明するブロック図である。 図１３に示したアナログ計算機の検査回路の構成の一例を説明する回路図である。 図１４に示した検査回路の動作の一例を説明するタイムチャートである。 図１の情報処理装置が図１３のアナログ計算機を制御する手順の一例を説明する流れ図である。 図１３に示したアナログ計算機の検査回路の構成の他の一例を説明する回路図である。
　以下、実施例を図面を用いて説明する。ただし、本発明は以下に示す実施の形態の記載内容に限定して解釈されるものではない。本発明の思想ないし趣旨から逸脱しない範囲で、その具体的構成を変更し得ることは当業者であれば容易に理解される。
　以下に説明する発明の構成において、同一部分又は同様な機能を有する部分には同一の符号を異なる図面間で共通して用い、重複する説明は省略することがある。
　本明細書等における「第１」、「第２」、「第３」などの表記は、構成要素を識別するために付するものであり、必ずしも、数または順序を限定するものではない。また、構成要素の識別のための番号は文脈毎に用いられ、一つの文脈で用いた番号が、他の文脈で必ずしも同一の構成を示すとは限らない。また、ある番号で識別された構成要素が、他の番号で識別された構成要素の機能を兼ねることを妨げるものではない。
　図面等において示す各構成の位置、大きさ、形状、範囲などは、発明の理解を容易にするため、実際の位置、大きさ、形状、範囲などを表していない場合がある。このため、本発明は、必ずしも、図面等に開示された位置、大きさ、形状、範囲などに限定されない。
　本明細書において単数形で表される構成要素は、特段文脈で明らかに示されない限り、複数形を含むものとする。
　実施例１では、本発明を適用した情報処理装置１００（図１）を説明する。情報処理装置１００はイジングモデルに相似した構造を持つイジングモデル計算機２００（図２）を備える。情報処理装置１００は本発明を適用することにより、イジングモデル計算機２００を構成するＭＯＳＦＥＴ（Ｍｅｔａｌ―Ｏｘｉｄｅ―Ｓｅｍｉｃｏｎｄｕｃｔｏｒ　Ｆｉｅｌｄ―Ｅｆｆｅｃｔ―Ｔｒａｎｓｉｓｔｏｒ）群（図６）の製造ばらつきを回避して情報処理を実現する。
　なお、本実施例を説明する上で、イジングモデル計算機２００が模擬しようとしているイジングモデルについて簡単に説明する。
　イジングモデルは磁性体の振舞いを説明するための統計力学のモデルである。イジングモデルは＋１／－１（ないしは，０／１，上／下）の２値をとるスピンと，スピン間の相互作用を示す相互作用係数，および，スピン毎にある外部磁場係数で定義される。
　イジングモデルは与えられたスピン配列，相互作用係数，および，外部磁場係数から，その時のエネルギーを計算することが出来る。イジングモデルのエネルギー関数は一般的に次式で表わされる。なお，σｉ，σｊはそれぞれｉ番目とｊ番目のスピンの値，Ｊｉ，ｊはｉ番目とｊ番目のスピンの間の相互作用係数，ｈｉはｉ番目のスピンに対する外部磁場係数，＜ｉ，ｊ＞は隣接する２つのサイトの組合せ，σはスピンの配列を表わすものとする。
Figure JPOXMLDOC01-appb-M000001
　イジングモデルの基底状態を求めるというのは，イジングモデルのエネルギー関数を最小化するスピンの配列を求める最適化問題である。
　図１はイジングモデル計算機２００（詳細は後述）を備え、イジングモデル計算機２００を制御して情報処理を行う情報処理装置１００の構成の一例を説明する図である。
　情報処理装置１００の典型例としては、現在一般的に用いられているパーソナルコンピュータに、イジングモデル計算機２００Ａ、Ｂで構成されたアクセラレータを装着したものである。なお、２個のイジングモデル計算機２００Ａ、Ｂを特に区別する必要が無い時には、単に符号２００として参照する。
　情報処理装置１００は、中央処理装置（ＣＰＵ）１１０、一時記憶装置（ＲＡＭ）１２０、固定記憶装置（ＨＤＤ）１６０、を有し、これらがシステムバス１３０で結合されている。また、システムバス１３０経由でネットワークインタフェースコントローラ（ＮＩＣ）１４０が接続され、情報処理装置１００は装置間ネットワーク１９０を経由して他の情報処理装置と通信を行うことが可能である。
　加えて、システムバス１３０に入出力コントローラ１５０が接続され、その先に１個ないしは複数個のイジングモデル計算機（図１の例では、２個のイジングモデル計算機２００Ａ、Ｂ）が接続される。
　入出力コントローラ１５０はシステムバス１３０上のプロトコルを、イジングモデル計算機２００のインタフェースに合せて変換するものである。情報処理装置１００のＣＰＵ１１０上で動作するソフトウェア（例えば制御プログラム１２１）は、一般的には特定のアドレスに対してリード／ライトを行うことで（いわゆるＭＭＩＯ、Ｍｅｍｏｒｙ Ｍａｐｐｅｄ Ｉ／Ｏ）、入出力コントローラ１５０を経由して、イジングモデル計算機２００を制御することが出来る。
　情報処理装置１００は、典型的には現在一般的に用いられているパーソナルコンピュータのようなものが一実装形態として考えられる。よって、情報処理装置１００はディスプレイなどの表示装置とキーボードなどの入力装置を備え、ユーザが対話的に操作することが可能なものである。そして、イジングモデル計算機２００などのアナログ計算機は、パーソナルコンピュータに対する演算加速装置（アクセラレータ）であり、一種の周辺機器として付加される。例えば、現在のパーソナルコンピュータでは表示装置での描画処理を高速化するためのアクセラレータとしてグラフィックスアクセラレータ（ＧＰＵとも呼ばれる）をＰＣＩ　Ｅｘｐｒｅｓｓ（登録商標）などのインタフェースで接続する拡張ボードとして付加することが多いが、イジングモデル計算機２００もこれに似た形態で実装される。イジングモデル計算機２００は例えば、１つまたは複数の半導体集積回路装置（イジングチップ）として構成することができる。
　情報処理装置１００でイジングモデル計算機２００に問題を解かせようとする場合には、例えばＨＤＤ１６０に問題データ１６１を格納しておく。あるいは、問題データ１６１は入力装置から直接入力してもよい。問題データ１６１とは、イジングモデル計算機が求解可能なイジングモデルであり、例えばイジングモデルを定義する相互作用係数と外部磁場係数のリストとしてＨＤＤ１６０上に保持される。つまり、問題データ１６１は数１に示した相互作用係数Ｊｉ，ｊと外部磁場係数ｈｉを列挙したものである。
　制御プログラム１２１はＣＰＵ１１０が実行するプログラムであり、問題データ１６１と特性データ１６２に基づいてイジングモデル計算機２００を制御するプログラムである。制御プログラム１２１の実行結果はイジングモデル計算機２００が生成した問題データ１６１で定義される問題に対する答えである。
　特性計測プログラム１２２はＣＰＵ１１０が実行するプログラムであり、イジングモデル計算機２００の製造ばらつきを計測するプログラムである。特性計測プログラム１２２は計測結果を特性データ１６２として記録する。
　図２は、本実施例のイジングモデル計算機２００の構成図の例である。イジングモデル計算機２００はスピンアレイ２１０、Ｉ／Ｏドライバ２２０、Ｉ／Ｏアドレスデコーダ２３０、および、相互作用アドレスデコーダ２４０、乱数発生器２５０、ビット調整器２５１から構成される。
　本実施例でのイジングモデル計算機は、現在広く用いられているＣＭＯＳ（Ｃｏｍｐｌｅｍｅｎｔａｒｙ　Ｍｅｔａｌ―Ｏｘｉｄｅ―Ｓｅｍｉｃｏｎｄｕｃｔｏｒ）集積回路として実装されることを想定して説明するが、他の固体素子でも実現可能である。本実施例では、スピンアレイ２１０は半導体メモリアレイで構成されている。メモリの種類は問わないが、ここでは高速動作可能なＳＲＡＭ（Ｓｔａｔｉｃ Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）を用いた例で説明する。イジングモデル計算機２００は、スピンアレイ２１０にリード／ライトを行うための、ＳＲＡＭ互換インタフェース２２２によって駆動可能である。ＳＲＡＭ互換インタフェースは、アドレスバス２９０、データバス２９１、Ｒ／Ｗ制御線２９３、および、Ｉ／Ｏクロック２９２を含む。また、イジングモデル計算機２００は、制御等のために不揮発メモリ等の他のメモリ２７０を備えていてもよい。
　イジングモデル計算機２００では、イジングモデルのスピンσｉ、相互作用係数Ｊｉ，ｊ、および、外部磁場係数ｈｉを、全てスピンアレイ２１０内のメモリセルに記憶する情報で表現する。これらの値のスピンアレイのメモリセルへのリード／ライトは、ＳＲＡＭ互換インタフェース２２２を介して、Ｉ／Ｏアドレスデコーダ２３０及びＩ／Ｏドライバ２２０で行う。
　基底状態を探索すべきイジングモデルをイジングモデル計算機２００に設定するために、スピンアレイ２１０内のメモリセルにはアドレスが付与されている。また、スピンアレイ２１０のスピンσｉ、相互作用係数Ｊｉ，ｊ、および、外部磁場係数ｈｉは、メモリセルのアドレスに対応付けられている。任意のアドレスのメモリセルに情報をリード／ライトする手法は、公知の半導体メモリの制御技術と同様である。
　なお、ＳＲＡＭ互換インタフェースを構成するアドレスバス２９０、データバス２９１、および、Ｒ／Ｗ制御線２９３は、Ｉ／Ｏクロック２９２に入力されるクロックに同期して動作する。但し、本発明においてインタフェースが同期式である必要性は無く、非同期式のインタフェースでも良い。
　また、イジングモデル計算機２００は基底状態探索を行うために、スピンアレイ２１０の内部でスピン間の相互作用を実現する。ここで、相互作用とは、スピンアレイ２１０のメモリセルの値を読み出し、読み出した値から当該メモリセルに対応するスピン近傍での局所的なエネルギーを計算し、エネルギーが小さくなるような当該スピンの値を決定し、決定した当該スピンの値をメモリセルに書き込む一連の処理をいう。
　この相互作用を外部から制御するのが相互作用制御インタフェースであり、具体的には相互作用を行うスピン群を指定するアドレスを相互作用アドレス２８０で入力し、相互作用クロック２８１に入力されるクロックに同期して相互作用を行う。内部的には、相互作用アドレスデコーダ２４０が、相互作用アドレス２８０と相互作用クロック２８１から、スピンアレイ２１０内の相互作用アドレス２８０に対応する一群のスピンを更新するためのクロックを生成する。
　なお、相互作用は必ずしもクロック同期回路で実現される必要はなく、非同期回路であっても構わない。この場合、相互作用クロック２８１の役割はクロックの入力ではなく、相互作用の実行を許可するイネーブル信号であるものとする。
　後述するスピン間の相互作用によるエネルギー最小化で、適用されたイジングモデルの基底状態探索を実現することが出来るが、これだけでは局所最適解に陥ってしまう可能性がある。後述の相互作用回路３１０では基本的に、局所的なエネルギーを小さくする方向の動きしかないため、一旦局所最適解に陥るとそこから抜け出すことが出来ず、大域最適解に到達しない。そのため、局所最適解から脱出するために、所定の確率でエネルギーを大きくする遷移を許容する処理を取り入れる。
　本実施例では、スピンの値を仮想的な温度Ｔで決まる確率で反転させることにした。ここで温度Ｔは基底状態探索の過程を物理的なアニーリングに例えたものである。基底状態探索の初期では高い温度とし、徐々に温度を下げながら局所的な探索を行い、最終的に温度がゼロとなる状態まで冷却していく。
　本実施例では、スピンの値を所定確率で反転させるために、乱数発生器２５０と乱数発生クロック２６０を用いる。乱数発生器２５０は乱数発生クロック２６０の１クロック毎に乱数を発生させる。ビット調整器２５１は、基底状態探索の初期では高い確率でスピンの値を反転させ、終期では低い確率でスピンの値を反転させるように、乱数発生器２５０からの出力ビットを調整するものである。具体的には、乱数発生器２５０の出力から所定ビット数を取り出し、多入力のＡＮＤ回路またはＯＲ回路で演算することにより、基底状態探索の初期では１が多く、基底状態探索の終期では０が多く発生するように出力を調整する。
　イジングモデル計算機２００内のスピンアレイ２１０は、イジングモデルに対応し、これを模擬するものである。スピンアレイ２１０は、１個のスピンとそれに付随する相互作用係数、及び、外部磁場係数の保持と、基底状態探索処理を実現するスピンユニット３００（図３）を基本構成単位としている。
　図３は、スピンアレイ２１０の単位構成要素の構成の一例を示すブロック図である。
　図４は、イジングモデルの一例であり、図３の単位構成要素が１つのスピンユニットを模擬する。
　図３に示すスピンユニット３００を複数個並べることで、図４に示すような３次元格子状のトポロジを持つイジングモデルを構成する。図４の例は、３（Ｘ軸方向）×３（Ｙ軸方向）×２（Ｚ軸方向）の大きさの３次元格子（合計１８スピン）である。
　座標軸の定義は図示した通り、図面右方向をＸ軸、図面下方向をＹ軸、図面奥行き方向をＺ軸としている（以下の図面でも同様）。この座標軸は実施例の説明上便宜的に必要なだけであり、発明とは関係しない。３次元格子以外のトポロジ、例えばツリー状のトポロジなどを利用する場合には、座標軸とは別にツリーの段数等で表現することになる。
　図４の３次元格子状のトポロジにおいて、スピン間の相互作用をグラフとしてとらえると、最大で次数５のスピン（頂点）が必要となる。なお、外部磁場係数の接続も含めて考えると、最大で次数６が必要となる。
　図３に示す１個のスピンユニット３００は隣接するスピン（例えば隣接するスピンが上側、下側、左側、右側、奥行側の５個の場合）σj，σｋ，σｌ，σｍ，σｎの値が入力される。また、スピンユニット３００はスピンσｉと外部磁場係数ｈｉに加え、前記した隣接するスピンとの相互作用係数であるＪj，i，Ｊk，i，Ｊl，i，Ｊm，i，Ｊn，i（隣接する５スピンとの相互作用係数）を保持するメモリセルを有している。なお、図３の例では相互作用係数を保持するメモリセルは入力側に置いているが、出力側に置くこともできる。
　ところで、イジングモデルは一般的に無向グラフで表現される相互作用を有している。前記した（数１）では、相互作用を表わす項として、Ｊi，j×σｉ×σｊがあるが、これはｉ番目スピンからｊ番目スピンへの相互作用を示している。この時、一般的なイジングモデルではｉ番目スピンからｊ番目スピンへの相互作用と、ｊ番目スピンからｉ番目スピンへの相互作用を区別することはない。つまり、Ｊi，jとＪj，iは同一である。また、一般的なイジングモデルでは、イジングモデル全体のエネルギーを数１に従って計算するときには、どちらか一方の方向しかエネルギーの計算に含めなくてよい。
　しかし、本実施例のイジングモデル計算機２００では、このイジングモデルを有向グラフに拡張し、ｉ番目スピンからｊ番目スピンへの相互作用と、ｊ番目スピンからｉ番目スピンへの相互作用を非対称にする例を示している。
　そのため、１個のスピンユニット３００をｉ番目スピンσｉを含むユニットと考えた時に、このスピンユニットが保持する相互作用係数であるＪj，i，Ｊk，i，Ｊl，i，Ｊm，i，Ｊn，iは、隣接するｊ番目、ｋ番目、ｌ番目、ｍ番目、ｎ番目のスピンσj，σｋ，σｌ，σｍ，σｎから、ｉ番目スピンσｉへの相互作用を決めるものである。
　このことは、図４において、スピンユニット３００に含まれている相互作用係数が対応する矢印（相互作用）が、図示されているスピンユニット３００の外部のスピンから、スピンユニット３００の内部のスピンに向かっていることに対応している。
　スピンユニット３００は、イジングモデルのスピンσｉ、相互作用係数Ｊj，i，Ｊk，i，Ｊl，i，Ｊm，i，Ｊn，iおよび、外部磁場係数ｈｉを保持するために、１ビットのメモリセルを複数個備えている。この１ビットのメモリセルは図３中に、スピンσｉを保持するＳＰＩＮ，外部磁場係数ｈｉを保持するＩＳ０，ＩＳ１，相互作用係数を保持するＩＵ０，ＩＵ１，ＩＬ０，ＩＬ１，ＩＲ０，ＩＲ１，ＩＤ０，ＩＤ１，ＩＦ０，ＩＦ１として図示されているものである。なお、ＩＳ０とＩＳ１、ＩＵ０とＩＵ１、ＩＬ０とＩＬ１、ＩＲ０とＩＲ１、ＩＤ０とＩＤ１、および、ＩＦ０とＩＦ１はそれぞれ２個１組で役割を果たすものであるため、それぞれまとめてＩＳｘ，ＩＵｘ，ＩＬｘ，ＩＲｘ，ＩＤｘ，および、ＩＦｘと略記する。
　スピンユニット３００が有する各メモリセルの構造の一例としては、公知のＳＲＡＭメモリのセルを利用することができる。すなわち、メモリセルはＣＭＯＳインバータ２個で構成されるデータ保持部を有する。インバータと直列に接続されるパスゲートトランジスタのゲートをワード線で制御し、ソース・ドレイン経路をビット線に接続することで、データ保持部へのデータ読み書きを実現する。もっとも、メモリセル構造はこれに限られず、少なくとも２値を記憶できる構成であればよい。例えば、ＤＲＡＭやフラッシュメモリのような他のメモリを用いることができる。
　図３の例では、各メモリセルは１ビット格納できるＳＲＡＭメモリセルとしている。しかし、多値化したメモリセルを用いてもよい。多値化したメモリセルは、物理的に１つのメモリセルを用いてもよいし、複数の１ビットメモリセルを用いて複数ビットを格納できるように構成してもよい。メモリセルを多値化すると、使用できる係数の値の自由度が広がる。一方、回路規模や処理時間は拡大する。したがって、解くべき問題の種類や用途に応じて、適切なメモリを選択すればよい。
　ここで、スピンユニット３００はｉ番目のスピンを表現するものとして説明を行う。メモリセルＳＰＩＮはスピンσｉを表現するためのメモリセルでありスピンの値を保持する。スピンの値はイジングモデルでは＋１／－１（＋１を上、－１を下とも表現する）であるが、これをメモリセルの２値である１／０に対応させる。本明細書では、＋１を１、－１を０に対応させることにするが、逆の対応でもかまわない。
　図５を用いて、スピンユニット３００が有するメモリセルと、図４に示したイジングモデルのトポロジとの対応関係を示す。ＩＳｘは外部磁場係数を表現する。また、ＩＵｘ，ＩＬｘ，ＩＲｘ，ＩＤｘ，ＩＦｘはそれぞれ相互作用係数を表現する。ＩＵｘは上側のスピン（Ｙ軸方向で－１）、ＩＬｘは左側のスピン（Ｘ軸方向で－１）、ＩＲｘは右側のスピン（Ｘ軸方向で＋１）、ＩＤｘは下側のスピン（Ｙ軸方向で＋１）、ＩＦｘは奥行き方向に接続するスピン（Ｚ軸方向で＋１ないしは－１）との相互作用係数を示している。また、イジングモデルを有向グラフとして捉えた場合に、あるスピンから見ると他のスピンが自スピンに及ぼす影響の係数を持つことになる。自スピンが他のスピンに与える影響の係数は、それぞれの他のスピンに属する。すなわち、このスピンユニット３００は最大で５個のスピンと接続される。１個のスピンユニットは５個のスピンとの相互作用と、１個の外部磁場の合計６個の影響によって、スピンの値を決める。
　本実施例のイジングモデル計算機２００では、外部磁場係数、及び、相互作用係数は、＋１／０／－１の３値に対応する。そのため、外部磁場係数、および、相互作用係数を表わすためには、少なくともそれぞれ２ビットのメモリセルが必要となる。図３に示すスピンユニット３００において、メモリセルＩＳｘ，ＩＵｘ，ＩＬｘ，ＩＲｘ，ＩＤｘ，ＩＦｘは、末尾の数字が０と１の２つのメモリセル（例えばＩＳｘの場合にはＩＳ０とＩＳ１）の組合せで、＋１／０／－１の３値を表現する。例えば、ＩＳｘの場合には、ＩＳ１で＋１／－１を表現し、ＩＳ１が１の時は＋１，０の時には－１を表現する。これに加えて、ＩＳ０が０の時には外部磁場係数を０と見なし、ＩＳ０が１の時にはＩＳ１で決まる＋１／－１のいずれかを外部磁場係数とする。外部磁場係数が０の時は外部磁場係数をディセーブルしていると考えれば、ＩＳ０は外部磁場係数のイネーブルビットであると言うことができる（ＩＳ０＝１の時に、外部磁場係数がイネーブルされる）。相互作用係数を表現するＩＵｘ，ＩＬｘ，ＩＲｘ，ＩＤｘ，ＩＦｘも同様に係数とビットの値を対応させている。
　スピンユニット３００は、Ｉ／Ｏアドレスデコーダ２３０とＩ／Ｏドライバ２２０で駆動、制御、読み出しすることで、一般的なＳＲＡＭと同様に、スピンユニット３００内のメモリセルをイジングモデル計算機２００のＳＲＡＭ互換インタフェースでリード／ライトすることが出来る。
　スピンユニット３００では隣接スピンとの間でエネルギーを最小化するようにスピンの次状態を決定するが、それは隣接スピンと相互作用係数の積、及び、外部磁場係数を評価した時に、正の値と負の値のどちらが支配的か判断することと等価である。例えば、ｉ番目スピンσｉに、スピンσj，σｋ，σｌ，σｍ，σｎが隣接しているとして、スピンσｉの次状態は以下のように決まる。まず、隣接スピンの値はσj＝＋１，σｋ＝－１，σｌ＝＋１，σｍ＝－１，σｎ＝＋１とし、相互作用係数はＪj，i＝＋１，Ｊk，i＝＋１，Ｊl，i＝＋１，Ｊm，i＝－１，Ｊn，i＝－１、外部磁場係数ｈｉ＝＋１とする。このとき、相互作用係数と隣接スピンの積、及び、外部磁場係数をそれぞれ並べると、σj×Ｊj，i＝＋１，σｋ×Ｊk，i＝－１，σｌ×Ｊl，i＝＋１，σｍ×Ｊm，i＝＋１，σｎ×Ｊn，i＝－１，ｈｉ＝＋１となる。外部磁場係数は、常に値が＋１のスピンとの相互作用係数と読み替えて良い。
　ここで、ｉ番目のスピンと隣接スピンとの間での局所的なエネルギーは、前述した係数にそれぞれｉ番目スピンの値を乗じて、さらに符号を反転させたものになる。例えば、ｊ番目スピンとの間での局所的なエネルギーは、ｉ番目スピンを＋１とした時には－１、ｉ番目スピンを－１としたときには＋１となるので、ｉ番目スピンを＋１にするほうが、ここでの局所的なエネルギーを小さくする方向に働く。
　このような局所的なエネルギーを全ての隣接スピン間と外部磁場係数について考えた時に、ｉ番目スピンを＋１／－１のどちらにしたほうがエネルギーを小さくできるかを計算する。これは、先程示した相互作用係数と隣接スピンの積、及び、外部磁場係数をそれぞれ並べたものにおいて、＋１と－１のどちらが多いか数えれば良い。先程の例では、＋１が４個、－１が２個である。仮に、ｉ番目スピンを＋１とするとエネルギーの総和は－２、ｉ番目スピンを－１とするとエネルギーの総和は＋２になる。よって、＋１の個数が多い時にはｉ番目スピンの次状態を＋１とし、－１の個数が多い時にはｉ番目スピンの次状態を－１にするという多数決で、エネルギーを最小化するｉ番目スピンの次状態を決定することが出来る。
　図３のスピンユニット３００に図示する相互作用回路３１０は、前記した動作を行うための回路を、論理回路として表現した場合の概略を示している。まず、隣接スピンの状態と相互作用係数の＋１／－１を示すメモリセルとの排他的論理和の否定（ＸＮＯＲ）で、その相互作用だけを見た時にエネルギーを最小化するスピンの次状態を計算することができる（＋１は１、－１は０にエンコードされているものとする）。もし、相互作用係数が＋１／－１だけであれば、この出力のうち＋１／－１のどちらが多いか多数決論理で判定すればスピンの次状態を決定することができる。外部磁場係数に関しては、常に状態＋１のスピンとの相互作用係数に相当するものと考えれば、単に外部磁場係数の値がスピンの次状態を決定する多数決論理回路３１１に入力すべき値となる。
　多数決論理回路３１１の出力は、次状態としてメモリセルＳＰＩＮに入力される。なお、前述のように、ビット調整期２５１の出力となるランダムパルス列は、条件に応じてメモリセルＳＰＩＮに入力される次状態をランダムに変更し、演算途中で局所解に陥ることを防ぐためのものである。
　次に、係数０は多数決論理に上記のような＋１／－１の投票をしない（相互作用回路３１０内のスイッチで表現される）ことで実現する。もし相互作用回路３１０を論理回路として構成する場合には、このスイッチに対応するようなディセーブル端子を各入力に対応して持つ多数決論理回路を構成する必要がある。
　相互作用回路３１０のような回路をＡＮＤ，ＯＲ，ＮＯＴのような論理ゲートの組み合わせだけで作ることは回路的なコストが大きい。よって、本実施例では図６～図８に示すような一部アナログ的な動作を持つ回路で、相互作用回路３１０を実現する例を示している。アナログ的な動作を持つ回路とは、連続的に変化する電気信号を取り扱う回路である。ここでは、スピンを上向き（＋１）に引っ張る力と、スピンを下向きに（－１）に引っ張る力を電流の引き込み量で表現し、その差を増幅することで、スピンの向きを決定する。つまり、スピン間の相互作用の大きさを電子回路上の電流の流量で相似的に表現していると言うことができる。
　図６は、アナログ的動作による相互作用回路３１０を構成する回路であり、スピンを＋１側に引っ張る力を電流の引き込み量として表現する回路（図６(a)、ＣＴ側回路）と、－１側に引っ張る力を電流の引き込み量として表現する回路（図６(b)、ＣＢ側回路と称する）である。
　図７は、アナログ的動作による相互作用回路３１０で、図６の回路と合わせて用いられる回路の構成の一例を説明する図である。
　図８は図６及び図７の回路によって実現される動作の一例を説明するタイムチャートである。
　図６記載のこの二つの回路と、この二つの回路からの出力を増幅する図７の回路を、図８に示すタイミングチャートに沿って動作させることによって、相互作用回路３１０に相当する機能を実現することが出来る。
　図６の回路の動作を説明する。説明においては、図６に示されるように、図３の各メモリセルに格納されている値を供給する信号を、そのメモリセルの名前の末尾にＴを付加した信号名で示す（例えば、ＩＳ０に格納されている値０／１はＩＳ０Ｔとして参照される）。また、その信号を反転した信号を、メモリセルの名前の末尾にＢを付加した信号名で示す（例えば、ＩＳ０Ｔの反転した信号をＩＳ０Ｂとする）。
　また、隣接するスピンの値をＮＵＴ（上側）、ＮＬＴ（左側）、ＮＲＴ（右側）、ＮＤＴ（下側）で表現する。さらに、奥行方向のＮＦＴがあるが、図面の紙幅の都合で割愛する。図６の回路図上はＮＵＴ～ＮＤＴに関する回路は同様の構成が繰り返し現れていることが見て取れるが、同様にＮＦＴに関する部分も構成できる。なお、それぞれの信号について、その反転した値をＮＵＢ，ＮＬＢ，ＮＲＢ，ＮＤＢと表現する（ＮＦＴに対応してＮＦＢもあるが、図では省略している）。
　図６において、ＣＴは、後で説明するＣＢと対をなす線である。ＣＴ及びＣＢはクロック信号ＣＬＫＩによって駆動されるＰ型ＭＯＳＦＥＴで電源電圧に接続され、ＣＬＫＩがＬＯＷの時にＰ型ＭＯＳＦＥＴがＯＮとなって、ＣＴ及びＣＢはプリチャージされる。
　ＣＴ及びＣＢは、いくつかのＮ型ＭＯＳＦＥＴを経由した後に、最下段において、クロック信号ＣＬＫＩによって駆動されるＮ型ＭＯＳＦＥＴでグラウンドに接地する。つまり、ＣＴ，ＣＢはＣＬＫＩがＬＯＷの時にプリチャージされ、ＣＬＫＩがＨＩＧＨの時に前述したいくつかのＮ型ＭＯＳＦＥＴを通って、プリチャージされた電荷がグラウンドに引き抜かれる。
　ＣＴ，ＣＢそれぞれにおいて、電源電圧からグラウンドに至る経路にあるＮ型ＭＯＳＦＥＴ群がＯＮになっている数が、ＣＴにおいてはスピンを＋１側に引っ張る力を表現し、ＣＢにおいてはスピンを－１が側に引っ張る力を表現する。
　まずスピンを＋１側に引っ張るＣＴ側に注目して説明する。例えば、外部磁場係数が＋１のとき、それはスピンを＋１側に引っ張る力であるから、このことを回路で表現するために、ＩＳ１Ｔで駆動するＮ型ＭＯＳＦＥＴをＣＴからグラウンドに電流を流すように配置する。同様に、相互作用係数について考えると、上側に隣接するスピンとの相互作用については、上側の隣接スピンが＋１（ＮＵＴが１）で相互作用係数が＋１（ＩＵ１Ｔが１）、ないしは、上側の隣接スピンが＋１（ＮＵＢが１）で相互作用係数が－１（ＩＵ１Ｂが１）の時にスピンを＋１側に引っ張る力となるので、この条件の時にＣＴからグラウンドに向けて電流が流れるように、Ｎ型ＭＯＳＦＥＴを配置して、ゲートに適切な信号を入力する。なお、相互作用係数が０の場合には、この電流の流れを阻止すればよいので、下段にＩＳ０Ｔで駆動されるＮ型ＭＯＳＦＥＴを設けている。ここでは上側に隣接するスピンとの相互作用について説明したが、それ以外の向きで隣接しているスピンとの相互作用も同様である。
　次にスピンを－１側に引っ張るＣＢ側について説明する。回路の基本的な考え方はＣＴ側と同じで、電流を流す条件だけが異なっている。つまり、Ｎ型ＭＯＳＦＥＴのゲートに入力される信号が違う。
　上記の構成、及び、動作によって、スピンを＋１側に引っ張る力の総和に比例してＣＴからグラウンドに電流が流れ、スピンを－１側に引っ張る力の総和に比例してＣＢからグラウンドに電流が流れる。その結果、ＣＴとＣＢの電圧はそれぞれ降下し、その降下量もスピンを＋１側（ＣＴ）ないしは－１側（ＣＢ）に引っ張る力に比例する。ここで、ＣＴとＣＢの電圧差を図７に示すラッチ型センスアンプで増幅し、その出力をスピンの次の値とする。なお、ＣＬＫＩＢはクロック信号ＣＬＫＩの反転である。
　図８に図６と図７の回路で実現される動作を、タイミングチャートとして説明する。図８のタイミングチャートは、図６と図７の回路によってになわれる相互作用の動作（この動作を行う期間を相互作用期間と称する）における、図６及び図７の主要信号線の電圧を縦軸に示し、横軸に時間経過を示すものである。説明上、相互作用期間を電流期間と判定期間の２フェーズに分けて説明する。
　まず、相互作用期間に入る前に、ＣＴ及びＣＢはプリチャージされており、電流期間の始点に置いて、両者とも信号のレベルはＨＩＧＨにある。この電流期間ではＣＬＫＩがＨＩＧＨになることにより、プリチャージは終了し、ＣＴ及びＣＢからＮ型ＭＯＳＦＥＴ群を通してグラウンドに電流が流れる。その結果、＋１側に引っ張る力（ＣＴ側）と－１側に引っ張る力（ＣＢ側）それぞれに応じて、ＣＴとＣＢに電圧降下が生じる。例えば、図８の例ではＣＢの方がより電圧降下が大きいと仮定する。つまり、－１側に引っ張る力の方が強いというケースである。
　その後、判定期間ではＣＬＫＩがＬＯＷになることで、次の電流期間に向けたプリチャージが行われることに加えて、判定期間の始点におけるＣＴ及びＣＢの電圧が図７に示すラッチ型センスアンプで増幅される。センスアンプはＣＬＫＩ信号がＨＩＧＨの時にＣＴ及びＣＢと接続しており、ＣＬＫＩがＬＯＷの時（つまりＣＬＫＩＢがＨＩＧＨの時）にアンプとして動作を行う。このセンスアンプによる増幅結果はＣＬＴとＣＬＢとして現れ、この場合にはＣＬＴがＨＩＧＨ，ＣＬＢがＬＯＷという信号により相互作用の結果を得る。図８のタイムチャートの場合には、電圧降下が相対的に大きいＣＬＢが増幅されてＬＯＷとなり、電圧降下が相対的に少ないＣＬＴが増幅されてＨＩＧＨになる。ＣＬＢ＝ＬＯＷ（ＣＬＴ＝ＨＩＧＨ）の時にスピンの値を＋１側にすべき（＋１側に引っ張る力が強い）であり、ＣＬＴ＝ＬＯＷ（ＣＬＢ＝ＨＩＧＨ）の時にスピンの値を－１側にすべき（－１側に引っ張る力が強い）ということになる。
　なお、図６～図８の回路の構成例は一例である。例えば図７～図８に限らず、ＣＴとＣＢの電圧差を測定する公知の種々の回路構成を採用することができる。
　発明者らが上記のアナログ回路を検討したところ、図６に示すＮ型ＭＯＳＦＥＴ群の製造ばらつきが問題となることがわかった。後述する理由によって、各ＭＯＳＦＥＴの閾値電圧はばらつきを持ち、その結果としてソース－ドレイン間の電流特性にもばらつきを生じる。このことは、図６の回路においてスピンを＋１側ないしは－１側に引っ張る力に、ＭＯＳＦＥＴ毎の差が出てしまうということである。例えば、仮に図６のＣＴ側において、ＮＵＴとＩＵ１Ｔがゲートに接続されているＭＯＳＦＥＴの方が他のＭＯＳＦＥＴと比較してＯＮ時のドレイン間電流（ＩＤＳ）が多く流れるように製造ばらつきがあった場合、ＮＵＴとＩＵ１Ｔに対応する相互作用（上側に隣接するスピンと＋１の係数で接続されているときに、スピンを＋１側に向けようとする力）は、他の相互作用と比較して大きく影響することになってしまう。
　特にこの問題が顕在化するのは、＋１側に向かう力と－１側に向かう力が拮抗した時（これ以降、この状態を拮抗状態と呼ぶ）である。拮抗状態では、スピンの向きを＋１にしても－１にしてもエネルギーに変化は無い。一般的には、探索空間を広げる観点から、拮抗状態ではスピンの向きをランダムに選ぶことも多い。もしくは、エネルギーに変化が無いため現状維持という意味で、前のスピンの値を保持するようにすることも多い。ところが、上記の製造ばらつきの問題から、ＣＴとＣＢの電圧降下の度合いが異なり、常にどちらか一方がより電圧が低くなる、つまり常に多数決の勝者になってしまうことが実験的にも分かっている。よって、拮抗状態のスピンの値（これ以降、これを拮抗値と呼ぶ）は製造ばらつきによって決まるといえる。
　拮抗状態の問題は図６～図８の例に限らず、特にアナログ計算機の本質的な課題であると言うことが出来る。なぜならば、前述したような製造ばらつきを極力抑えたとしても、拮抗状態では微小なばらつきがセンスアンプの出力結果ではスピン値が＋１もしくは－１のどちらかに固定的に決まるという現象として現れるからである。
　トランジスタ（ここでは主にＭＯＳＦＥＴが用いられる）の製造ばらつき要因として、主に製造プロセス上の寸法誤差、ＲＤＦ（Ｒａｎｄｏｍ　Ｄｏｐａｎｔ　Ｆｌｕｃｔｕａｔｉｏｎ、詳細は例えば非特許文献１を参照）、及びＲＴＮ（Ｒａｎｄｏｍ　Ｔｅｌｅｇｒａｐｈ　Ｎｏｉｓｅ、詳細は例えば非特許文献２を参照）がある。
　製造プロセス上の寸法誤差はフォトリソグラフィ等の技術によって半導体基板上に回路を構成する幾何学的図形を形成する際の誤差により、トランジスタの特性が異なってしまうという問題であり、製造時に決定する静的な誤差である。
　ＲＤＦは半導体製造時に不純物（ＤｏｎｏｒないしはＡｃｃｅｐｔｏｒ、総称してＤｏｐａｎｔ）をドーピングするときに生じる、不純物原子の位置と密度のランダムなばらつきのことで、これによってトランジスタの閾値電圧にばらつきが生じる。ＲＤＦによるばらつきは製造時に決まり、その後トランジスタの寿命の中で変化することはほとんど無い。また、不純物原子の位置と密度のランダムなばらつきに起因するという点から、トランジスタ自体の大きさが大きい場合には、その影響は相対的に小さく見え、微細化の進展と共にその影響が大きく見えてくる。微細化プロセスの中でも、特に閾値電圧ばらつきを低く抑える必要があるトランジスタについては、ゲート長を長くするなど、ゲートの面積を大きく取って、ＲＤＦの影響を緩和する策がとられる。
　ＲＴＮはトランジスタのゲートにある絶縁膜（酸化膜など）にキャリアが捕縛、解放されることにより、動的に閾値電圧ばらつきが生じる現象である。キャリアの捕縛、解放が起こる時定数は様々であり、様々な時定数のＲＴＮ現象が合わさって閾値電圧ばらつきとして観測される。その時定数は製造時に何らかの要因で決まると考えられているが、時定数が極めて短いものから長いものまであるため、製造時のスクリーニングを困難にする要因となる。ＲＤＦとの違いとして、微細化に伴って影響を及ぼす度合いがＲＤＦよりＲＴＮの方が強いこと、閾値電圧の動的な変化を及ぼすことがあげられる。
　本発明者らが検討したところによれば、本実施例で示したイジングモデル計算機３００を、６５ｎｍ程度の世代のプロセスで製造した場合には、寸法誤差やＲＤＦによると推定される静的な影響が主体的であり、ＲＴＮの影響と考えられるような動的なばらつきはまだ顕在化していない。
　そこで、本発明者らは、各個体の特性を事前に計測し、その計測結果に基づいて各個体に入力する問題を改変することで、各個体を構成する素子の特性ばらつきに適応するという技術を考案した。また、この方法では各個体の特性の計測は基本的に一度だけ行えばよいため、工場出荷時に計測して、その情報と対応する個体を共に顧客に提供するといった形態も可能である。
　本実施例における特性ばらつきへの対処は、イジングモデル計算機２００を制御する情報処理装置１００の制御プログラム１２１及び特性計測プログラム１２２によって実現される。
　図９は制御プログラム１２１の動作を示すフローチャートである。このフローチャートに基づく動作は以下の通りである。これらの動作は、ソフトウェアである制御プログラム１２１による情報処理が、ＣＰＵ１１０、ＲＡＭ１２０、ＨＤＤ１６０等、情報処理装置１００のハードウェア資源を用いて実現している。全体の制御や計算はＣＰＵ１１０が行い、ＲＡＭ１２０、ＨＤＤ１６０へのリード・ライトも制御する。これらは、以下の他のフローチャートでも同様である。
　ステップ９０１：本フローチャートはイジングモデルをイジングモデル計算機２００で求解するための手順である。本フローチャートの実行に先立ち、ユーザは、入力装置を経由して入力する等して、問題データ１６１をＨＤＤ１６０上に用意したものとする。なお、情報処理装置１００は複数のイジングモデル計算機２００Ａ，２００Ｂを有しているが、ユーザが問題データ１６１を解かせるイジングモデル計算機を自分で選択するか、ステップ９０１の段階でＣＰＵ１１０が適宜選択するものとする。これ以降の説明でイジングモデル計算機２００は、ステップ９０１で選択されたイジングモデル計算機のことを称するものとする。
　ステップ９０２：情報処理装置１００が備えるイジングモデル計算機２００に対応する特性データ１６２がＨＤＤ１６０上に存在するかを検査する。存在しない場合、ステップ９０３で特性データ１６２を生成する。存在する場合、ステップ９０３を実行せず、ステップ９０４に遷移する。なお、特性データ１６２の詳細は後述する。
　ステップ９０３及びステップ９０４：特性計測プログラム１６２を実行して、イジングモデル計算機２００Ａ、２００Ｂに対応する特性データ１６２を生成する。なお、特性計測プログラム１６２の詳細は後述する。
　ステップ９０５：特性データ１６２を元にして、特性ばらつきを回避できるように問題データ１６１を変換する。このステップで行われる変換の詳細は後述する。変換した問題データは変換後問題データと称する。
　ステップ９０６：イジングモデル計算機２００に変換後問題データを入力する。具体的には、変換後問題データを構成する相互作用係数と外部磁場係数を、イジングモデル計算機２００を構成するメモリセル群に書き込む。この書込みは入出力コントローラ１５０がイジングモデル計算機２００のＳＲＡＭ互換インタフェース２２２を制御することで実現される。そのため、ステップ９０６ではＣＰＵ１１０がＭＭＩＯ等の手段によって、入出力コントローラ１５０を制御する。
　ステップ９０７：イジングモデル計算機２００に対して、求解を指示する。具体的には、イジングモデル計算機２００の相互作用アドレス２８０、相互作用クロック２８１、及び、乱数発生クロック２６０を制御する。この制御も直接は入出力コントローラ１５０によって行われるので、ＣＰＵ１１０はＭＭＩＯ等の手段によって、入出力コントローラ１５０を制御する。
　ステップ９０８：イジングモデル計算機２００が求解した結果、すなわちスピン配列を読み出す。ステップ９０６と同様にＳＲＡＭ互換インタフェース２２２を制御するために、ＣＰＵ１１０がＭＭＩＯ等の手段で、入出力コントローラ１５０を制御する。
　ステップ９０９：ステップ９０８で読みだしたスピン配列を変換し、問題データ１６１に対する解を得る（逆変換後スピン配列と称する）。イジングモデル計算機２００で解いた問題は変換後問題データであるから、ステップ９０８で得られる解は変換後問題データの解である。この解を、変換前の問題データ１６１の解として解釈する作業が必要である。この解釈によって得られたスピン配列が逆変換後スピン配列である。解釈の方法の詳細は後述する。
　ステップ９１０：逆変換後スピン配列を求解結果としてユーザに呈示する。また求解結果は、情報処理装置１００のＨＤＤ１６０に格納してもよいし、そのまま出力装置を介して外部に出力してもよい。
　図１０は特性データ１６１の構成の一例である。特性データ１６１は本実施例の場合、各スピンの拮抗値を示したリストである。拮抗値はイジングモデル計算機の個体毎に異なるため、特性データ１６１は個体毎のエントリ（個体情報エントリ１０１０，１０２０）を持つ。例えばイジングモデル計算機２００Ａの特性は個体情報エントリ１０１０に、イジングモデル計算機２００Ｂの特性は個体情報エントリ１０２０に格納される。
　各エントリは、例えば個体情報エントリ１０１０の例で説明すると、シリアル番号部１０１１、タイムスタンプ１０１２、拮抗値記録部１０１３のよう構成をとる。個体情報エントリ１０１０はイジングモデル計算機の個体に対応することからイジングモデル計算機の個体を識別する情報、すなわちシリアル番号を個体情報エントリに含めることで、個体情報エントリ１０１０と個体との紐付けを行う。また、本発明者らの検討したところによれば、拮抗値の特性は製造時にほぼ決定されるが、経年変化が皆無とは言い切れないので、タイムスタンプ１０１２に個体情報エントリ１０１０を構成した情報を計測した日付及び時刻を格納する。拮抗値は拮抗値記録部１０１３に、イジングモデル計算機２００の各スピンの拮抗値を、各スピンの座標（本例ではＸ座標、Ｙ座標、Ｚ座標）と共に記録している。
　図１１に特性データ１６１を生成するための特性計測プログラム１６２の動作をフローチャートとして示す。このフローチャートに基づく動作は以下の通りである。
　ステップ１１０１：本フローチャートはイジングモデル計算機２００の有する各スピンの拮抗値を計測するための手順である。なお、情報処理装置１００は複数のイジングモデル計算機２００Ａ，２００Ｂを有しているが、本手順はイジングモデル計算機それぞれに対して実行される（並行して実行しても良い）。
　ステップ１１０２：拮抗値を測定するための問題（測定用モデル）を作成する。拮抗状態を引き起こさせるためのモデルの例を図１２に示す。
　図１２のモデルは、あるスピンの拮抗値を知るために、その両端（Ｙ軸方向＋１とＺ軸方向＋１）のスピンとそれぞれ＋１及び－１の相互作用を設け、かつ、両端のスピンの値はそれぞれ＋１にして拮抗状態を生成するものである。このモデルにおいて、温度０（ランダムなスピン反転を行わない状態）の元に１回相互作用を行うと、注目しているスピンはその製造ばらつきによって＋１ないしは－１に遷移する。なお、１個のスピンの拮抗値を計測するために、例えば図１２に示すようにＹ方向とＺ方向に隣接するスピンをそれぞれ１個ずつ用いると、それらのスピンを除いた他のスピンは、同様の相互作用を１つの測定用モデルに含めることで同時に計測することが出来る。例えば、図４に示したイジングモデルであれば、２個の測定用モデルを順次用いることで全てのスピンを測定できる。つまり１個目の測定用モデルでは、座標（Ｘ，Ｙ，Ｚ）＝（０，０，０），（１，０，０），（２，０，０），（０，１，１），（１，１，１），（２，１，１），（０，２，０），（１，２，０），（２，２，０）を測定対象とする。そして、２個目の測定用モデルでは、座標（Ｘ，Ｙ，Ｚ）＝（０，０，１），（１，０，１），（２，０，１），（０，１，０），（１，１，０），（２，１，０），（０，２，２），（１，２，２），（２，２，２）を測定対象とする。
　ステップ１１０３：ステップ１１０２で作成した測定用モデルをイジングモデル計算機２００に入力する。入力の方法はステップ９０６と同様である。
　ステップ１１０４：イジングモデル計算機２００に、当該測定用モデルで測定対象となっているスピンそれぞれについて、相互作用を実行するように求解を指示する。測定対象となっているスピンは、図１２で説明したように、測定用モデルにより拮抗状態に設定されている。なお、温度は０とする。すなわち、乱数発生クロック２６０を止めるか、乱数発生器２５０ないしはビット調整器２５１の設定を変更するなどして、ランダムなスピン反転は抑止する。
　ステップ１１０５：イジングモデル計算機２００からスピン配列を読み出す。その方法はステップ９０８と同様である。
　ステップ１１０６：ステップ１１０５で読みだしたスピン配列のうち、ステップ１１０２で生成した測定用モデルで測定対象となっているスピンの値を当該スピンの拮抗値として特性データ１６１に記録する。
　ステップ１１０７：全てのスピンの拮抗値を測定したかを判定する。もし、全てのスピンの測定が完了していなければ、その残りのスピンを測定するための測定モデルを生成して、再度測定を行うために、ステップ１１０２にも戻る。全てのスピンの測定が完了したら、ステップ１１０８に遷移する。
　ステップ１１０８：特性データ１６１の生成を完了した後に、特性計測プログラム１６２は終了する。特性データ１６１は、制御プログラム１２１が利用できるように、ＲＡＭ１２０あるいはＨＤＤ１６０に保持する。あるいは、イジングモデル計算機２００内のフラッシュメモリ等の不揮発メモリ２７０に格納しておいてもよい。イジングモデル計算機２００に自分の特性データ１６１を持たせておけば、イジングモデル計算機２００を他の情報処理装置で使用する場合でも、特性データ１６１を利用することができる。
　図９のステップ９０５に示した、問題データの変換について説明する。特性データ１６１に示された拮抗値と、拮抗時にあるべきスピン値（以降、希望値と呼ぶ）に基づいて、イジングモデル（問題データ）を変換し、ステップ９０９との組合せで、拮抗時のスピン値が希望値になるようにモデルを変換する。ここで、ｉ番目スピンの拮抗値をｕｉ、ｉ番目スピンの希望値をｅｉとする。この時、相互作用係数Ｊｉ，ｊはＪｉ，ｊ×ｕｉ×ｕｊ×ｅｉ×ｅｊに変換する。また、外部磁場係数ｈｉはｈｉ×ｕｉ×ｅｉに変換する。
　ステップ９０５と対になる変換としてステップ９０９では、ステップ９０８で読みだしたスピン値から、変換前のモデルの解に変換する。読みだしたスピン値をａｉ、変換前の問題データの解としてのスピン値をσｉとした時、σｉ＝ａｉ×ｕｉ×ｅｉとする。
　具体的な例を説明する。図１２に示すようなモデルで、上記のステップに基づいて、製造ばらつきに関わらず拮抗時には希望値が得られることを示す。３個のスピンσ１、σ２、σ３があり、σ１の値を決める相互作用を考える。σ２とσ３は両方とも＋１になっているものとする。拮抗状態を作るために相互作用係数Ｊ１，２は＋１、Ｊ１，３は－１とする。測定によって得られた３個のスピンそれぞれの拮抗値はｕ１＝＋１，ｕ２＝＋１，ｕ３＝＋１，希望値はｅ１＝－１，ｅ２＝＋１，ｅ３＝＋１とする。このモデルをそのままイジングモデル計算機に解かせると、σ１＝ｕ１＝＋１という答えが得られる。そこで、このモデルをステップ９０５に示すように変換すると、Ｊ１，２＝＋１×＋１×＋１×－１×＋１＝－１、Ｊ１，３＝－１×＋１×＋１×－１×＋１＝＋１というように変換される。この変換したモデルをイジングモデル計算機３００に入力すると、やはり拮抗するのでａ１＝＋１という結果が得られる。この結果をステップ９０９に示すようにσ１＝＋１×＋１×－１と変換してσ１＝－１を得る。
　つまり、ステップ９０５とステップ９０９で行う変換は、素子の製造ばらつきによって拮抗値が決まっているときに、拮抗状態になった時の値をステップ９０９での後処理によって希望値に変換できるように、ステップ９０５で係数の正負を入れ替えていると言うことができる。そして、ステップ９０９で解となるスピンの正負を入れ替えることにより希望値を得ることができる。
　この方式では、希望値が決まっている必要がある。実用上は、各スピンの外部磁場係数ｈｉが正であれば希望値を＋１、負であれば希望値を－１とする。これは、相互作用及び外部磁場の総和として拮抗状態が生じた場合に、外部磁場を優先するという考え方に基づいている（もしくは、外部磁場が相互作用より微小量大きいという考え方もできる）。希望値が決まっている場合には、その値を例えば問題データ１６１と関連付けてＨＤＤ１６０に保持しておけばよい。あるいは、問題データ１６１とは別に入力して、制御プログラム１２１が利用できるようにＲＡＭ１２０に保持しておいてもよい。
　もしくは、イジングモデル計算機２００のようにランダム性を利用する解法では、一般的に同じ問題を複数回解いて、その中の最良解を得ることが多い。そのような場合には問題を解く毎に希望値をランダムに決めて、複数回問題を解くことを試行する。この場合も、本変換を用いず拮抗状態で拮抗値に固定されてしまう場合と比べて、大きな空間を探索することができるため、製造ばらつきの影響を回避して解の品質向上に寄与することができる。ランダム性を利用する場合には、予めランダムなデータを保持しておいてもよいし、その都度生成してもよい。
　また、希望値を予め決めない別の方式として、問題を１回解いている間に、その過程において希望値をランダムに変更する方式がある。このように求解の途中において、異なる希望値を与えて変換した問題データを入力すれば、一回の求解において、その過程において拮抗状態が生じることで探索すべき解空間が狭まってしまうことを防ぐことができる。この場合、基本的な処理の流れは図９に示す例と同様であるが、問題データを変換する処理９０４の間、途中で希望値がランダムに変更される。スピン配列を変換する処理９０９でも、問題データに対応した変換処理が行われる。希望値を変更するタイミングは任意であるが、変更の周期を短くすれば効果が大きくなる反面、処理が複雑になる。
　以上では希望値をランダムに変えることを想定したが、希望値が２値を取る場合は、交互に変更する方式も考えられる。
　なお、本実施例ではイジングモデルと相似な構造を持つイジングモデル計算機２００を例として実施例を示したが、ニューラルネットワークと相似な構造を持つニューロ計算機についても同様に実施することができる。ニューラルネットワークでは、他のニューロンの発火状態（例えば、数値の０と１で表現される）を受け取り、自ニューロンが発火するかどうかを、他のニューロンの発火状態とニューロン間のシナプス結合係数を乗じたものの総和を計算し、総和を閾値関数に通してその出力で決定する。この動作をアナログ電子回路で実装しようとすると、本実施例のイジングモデル計算機と類似した構造を持つ（発火状態をスピン値、シナプス結合係数を相互作用係数と読みかえればイジングモデルとほぼ同じ構造となる）。
　ニューロ計算機の場合も相互作用回路（シナプス結合回路）の構成は、基本的に図６に示した構成が使用可能であり、スピン値ＮＵＴ等の代わりに発火状態、相互作用係数ＩＵ１Ｔの代わりにシナプス結合係数を入力すればよい。
　実施例２は、イジングモデル計算機２００で拮抗状態以外のばらつきも考慮する場合も考える。実施例１では、イジングモデル計算機２００の相互作用回路３１０におけるＭＯＳＦＥＴの製造特性ばらつきに起因する、相互作用が拮抗したときに顕在化するスピン値が固定的な方向に倒れてしまう問題に対応する実施例を示した。実施例１で示した方式は、相互作用回路３１０にそれを補償するための付加的な回路を加えること無く、イジングモデル計算機２００とそれを含んだ情報処理装置１００の制御で対応する方法であった。
　実施例２では、さらに製造ばらつきの影響が大きくなり、その影響が拮抗状態のみならず、それ以外の相互作用にも影響を及ぼすような場合に、対応できる構成の例を開示する。
　一般的にＭＯＳＦＥＴのドレイン間電流（ＩＤＳ）のばらつきは、チャネル長（Ｌ）とチャネル幅（Ｗ）に依存しており、１／√（Ｗ×Ｌ）に比例することが知られている。そこで、ＩＤＳのばらつきを１／ｎに抑制するためには、（Ｗ／Ｌ）をｎ２倍にしなければならない。このことは、ばらつきの抑制と微細化（高集積化）が相反関係になることを示している。
　一方で、実施例１に示したイジングモデル計算機２００よりも、複雑なイジングモデル（相互作用の関係を増やす、ないしは、相互作用係数・外部磁場係数の値域を広げる）に対応しようとすると、ばらつきの抑制に対してより厳しい要求が科される。
　例えば、実施例１のイジングモデル計算機２００は、相互作用係数及び外部磁場係数の値域は＋１～－１であった。また、相互作用と外部磁場を併せて、１個のスピンには６種類の力がかかっている（相互作用の上下左右奥行方向で５種に加えて外部磁場で合計６種）。この時、＋１に向ける力と－１に向ける力が半々の時に拮抗状態となり、この拮抗状態における製造ばらつきへの対応が実施例１であった。
　ここで、例えば、＋１に向ける力と－１に向ける力が半々で無い時は、当然より力が大きい方向にスピンが向くことが期待される。例えば、６種の力のうち外部磁場が無く、相互作用５方向で、＋１側に３、－１側に２の力で引っ張られるという３対２のケースを考える。この時は、３の力がある＋１側になるべきである。つまり、相互作用回路３１０ではＣＴ側回路でより多くの電流が流れて、ＣＴの電圧降下の幅が大きいはずである。
　しかし、そのためには、ＭＯＳＦＥＴの特性ばらつきが十分小さい範囲に収まっている必要がある。例えば相互作用回路３１０のＣＴ側回路とＣＢ側回路で最大に流れる電流（同時にＯＮになりうる全てのＭＯＳＦＥＴをＯＮにした時の電流）を１とする。その時に、３対２であれば、理想的にはＣＴ側回路は（１／６）×３＝１／２の電流が、ＣＢ側回路は（１／６）×２＝１／３の電流が流れる。つまり、ＣＴ側回路とＣＢ側回路の電流差は１／６であるため、このケースではＣＢ側回路のＭＯＳＦＥＴの電流ばらつきの結果、１／６以上の電流をさらに流してしまうと正しく動作しない。もしくは、ＣＴ側の電流がばらつき起因で流れる量が少なく、ＣＢ側の方は逆に多いというような場合も問題になる。
　そのため、上記の問題が生じないように、ＩＤＳのばらつきを抑制できるようにＭＯＳＦＥＴの大きさ（Ｗ，Ｌ）を決めて設計するわけだが、それも相互作用や外部磁場を複雑化すると問題になってくる。例えば、実施例１のように、相互作用係数と外部磁場係数が＋１～－１で、相互作用係数と外部磁場係数を併せて６種の力がある場合に、上記の問題が生じない程度にＩＤＳのばらつきをおさえることに成功した設計があるとしよう。この設計を元に、前述したばらつきとチャネル幅・チャネル長の関係で、複雑化させた設計で面積効率がどの程度落ちるかを示す。
　例えば、係数の幅を＋３～－３にすると、現状の設計と比べてＣＴ側回路、ＣＢ側回路それぞれ３倍の数のトランジスタが必要になる。現状では６種の力それぞれに４個のＭＯＳＦＥＴ（但し、同時にＯＮになるのはそのうち２個）を用いているが、１種に４×３＝１２個のＭＯＳＦＥＴ（但し、同時にＯＮになるのは最大６個）を必要とする。この構成では、拮抗状態が起きる時には係数＋３／－３をそれぞれ３方向ずつで使う９対９のケースになる。そこで、拮抗状態一歩手前の９対８のケースについて考えると、ＣＴ側回路は（１／１８）×９＝９／１８の電流が、ＣＢ側回路は（１／１８）×８＝８／１８の電流が流れることになるので、その差は１／１８である。係数＋１～－１の場合と比べると、差が１／３に縮小している。つまり、ＭＯＳＦＥＴの特性ばらつきも１／３に収めなければならず、そのためにはＭＯＳＦＥＴのサイズを９倍にしなければならない。
　このことから、相互作用が複雑になると、ＭＯＳＦＥＴの特性ばらつきを抑えるためにＭＯＳＦＥＴを多くしなければならない上に、元々複雑な相互作用を実現するためにＭＯＳＦＥＴの数が増えている。例えば上記の場合だと特性ばらつきを抑えるためのサイズ９倍に加えて、元々複雑な相互作用を実現するために３倍の数のトランジスタを用いているので、合計２７倍の面積を必要とする。このため、スケーラビリティを著しく阻害する要因となる。
　そのため、本実施例では拮抗状態以外でもＩＤＳのばらつきが顕在化する場合に対応する例を示す。例えば拮抗状態ではない３対２のケースであったとしても、ＣＴ側ないしはＣＢ側のＭＯＳＥＦＴのＩＤＳばらつきによって、３の力で引く側よりも２の力で引く側が強いと誤認されてしまうようなケースに対応するため例を示す。
　図１３は本実施例の実施の形態であるアナログ計算機１３００の構成図である。このアナログ計算機１３００は実施例１のイジングモデル計算機２００Ａと同様に、情報処理装置１００のような装置の一部に組み込まれて利用することを想定している。
　図１３は説明の都合上、実施例１のイジングモデル計算機２００を抽象化したアナログ計算機として説明する。そのアナログ計算機１３００は、アナログ計算ユニット１３１０Ａ～Ｅ（Ａ～Ｅを特に区別する必要が無い時は、単に１３１０として参照する）、アナログ演算回路１３１１、記憶素子１３１２、検査回路１３１３、ユニット間インタフェース１３２０、検査インタフェース１３３０から構成されている。
　これは実施例１に照らし合わせれば、イジングモデル計算機２００がアナログ計算機１３００に対応し、スピンユニット３００がアナログ計算ユニット１３１０に対応する。相互作用回路３１０が、例えば図６と図７で具体例を示したアナログ演算回路１３１１に対応する。また、スピンユニット３００内のメモリセルアレイが記憶素子１３１２に対応し、スピンユニット間で隣接スピンユニットのスピン値を伝送する伝送線がユニット間インタフェース１３２０に対応する。なお、ユニット間インタフェース１３２０及び検査インタフェース１３３０は図１３上ではバス状のトポロジになっているが、このトポロジは任意である。実施例１のイジングモデル計算機２００の場合には、ユニット間インタフェースが３次元格子状のトポロジを持っていると言うことが出来る。
　検査回路１３１３は、アナログ演算回路１３１１の処理結果が誤る可能性がある場合を検出して、その結果を検査インタフェース１３３０に出力する。アナログ演算回路１３１１の処理結果が誤る場合とは、２種類のケースがある。第一には、既に述べた拮抗状態である。第二には、拮抗状態ではないものの（つまりアナログ演算回路１３１１の製造ばらつきが一切なければ誤ることはないものの）、＋１側に引く力と－１側に引く力が近いため、製造ばらつきによっては誤判定する可能性があるケースである。どの程度の製造ばらつきが生じるかは、プロセスの世代やＭＯＳＦＥＴの設計パラメータからある程度推定することが可能である。
　一般的に、アナログ演算回路１３１１のようなアナログ電子回路で生じる各種ばらつきは、それを打ち消すための回路を付加して補償することが多く行われているが、このような方法では、元の回路の２倍近い物量が必要となることも多く、スケーラビリティを下げる原因となる。
　そこで、本実施例では後述するように、アナログ計算機１３００に組み込まれた簡単な検査回路１３１３と、アナログ計算機１３００を含む情報処理装置１００の全体制御の組み合わせで、アナログ計算機１３００の物量増加を抑えて、前述したような状況に対処する。
　図１４に検査回路１３１３の構成の一例を示す。検査回路１３１３では前述した２ケースを一括して検出するために、ＣＴ側回路とＣＢ側回路の電圧差が一定の範囲内に収まっているかどうかを検出する。一定の範囲内に収まった場合には、誤動作の可能性ありと判断する。具体例としては、ＣＴ側回路とＣＢ側回路の電圧差が所定閾値より小さい場合に誤動作の可能性ありと判定する。電圧差を許容する範囲は、どの程度の製造ばらつきが生じるかの見積によって設定できる。
　図１４の検査回路１３１３は図７のラッチ型センスアンプを基本構造としている。センスアンプは２種類の信号（ＣＴとＣＢ）の電圧差を増幅する役割だが、図１４に示すＭＯＳＦＥＴ１４０１の駆動力を下げると、電圧差の増幅が緩やかとなり、かつ、その勾配は電圧差が小さいほど緩やかなものとなる。この緩やかな勾配を検出することで、電圧差が許容範囲内かどうかを検出する。
　図１５は検査回路１３１３を動作させた時のタイムチャートであるが、これを参照しながら検査回路１３１３の構造と動作を詳説する。図８のタイムチャートと同様に、プリチャージしたＣＴ／ＣＢから隣接スピン値、相互作用係数、及び、外部磁場係数でＯＮ／ＯＦＦが決まるＭＯＳＦＥＴ群で電流を引き抜き、ＣＴ／ＣＢに電圧降下を生じさせる。その後、ＣＬＫＩＢをＨＩＧＨとしてセンスアンプを作動させるが、前述したようにＭＯＳＦＥＴ１４０１の駆動力をおさえているため、ＣＬＴ、ＣＬＢに現れる電圧の傾きは緩やかなものとなる。さらに、ＣＴとＣＢの電圧差が大きい場合と小さい場合を比較すると、小さい場合の方がより緩やかな電圧変化を見せる。よって、ＣＴとＣＢの電圧差はＣＬＴとＣＬＢの電圧変化の傾きとして現れ、その傾きが一定より緩いことを検出することで、検査回路１３１３を実現する。
　傾きが一定より緩いことを検出するために、本実施例では閾値を持つバッファと排他的論理和（ＸＯＲ）回路を用いる。ＸＯＲからの出力ＯＵＴＡを、ＣＬＫＩＢを微小時間Δｔ遅延させたタイミングで観測する。観測されたＯＵＴＡのレベルがＨＩＧＨかＬＯＷかで誤動作の可能性を検査することができる。このタイミングで観測するためには、例えばエッジトリガ型のフリップフロップ（ＦＦ）をΔｔ遅延させたＣＬＫＩＢで駆動すれば良い。
　なお、検査回路１３１３の構造が図７のセンスアンプと類似していることから分かるように、センスアンプと兼用することも可能である。なお、図１４の例は一例であり、センスアンプと全く別個にＣＴとＣＢの電圧差を測定する回路を設けてもよい。この場合は、電圧を測定する公知の種々の回路構成を採用することができる。
　図１６にこの検査回路１３１３を用いて、アナログ計算機１３００の誤動作に対応する情報処理装置１００の動作をフローチャートとして示す。
　ステップ１６０１：本フローは、実施例１の図９に示した情報処理装置１００によるイジングモデル計算機２００の制御（本実施例ではアナログ計算機１３００の制御）のうち、ステップ９０７で行われる処理を、検査回路１３１３を用いて制御する詳細なフローを示すものである。
　ステップ１６０２：アナログ計算機１３００が有する各アナログ計算ユニット１３１０Ａ～１３１０Ｂをそれぞれ一回ずつ動作させて、演算を実行する。
　ステップ１６０３：ステップ１６０２で行った演算の結果、検査回路１３１３でアナログ演算ユニット１３１１の誤動作の可能性が検出されていないかを検査する。誤動作の可能性が検出されていなければ、ステップ１６０７に進む。誤動作の可能性が検出されていれば、ステップ１６０４に進む。
　ステップ１６０４：誤動作があった場合、検査回路１３１３から検査インタフェース１３３０を介して得られた情報で、誤動作があった可能性のあるアナログ計算ユニット１３１０を特定し、当該ユニット及び当該ユニット周辺の記憶素子の値を読み出す（例えば、実施例１のイジングモデル計算機２００であれば、隣接スピンのスピン値だけを読み出せばよい）。
　ステップ１６０５：誤動作があったアナログ計算ユニット１３１０が行うべき演算を、ステップ１６０４で得た情報を用いて、情報処理装置１００のＣＰＵ１１０上で代替して行う。
　なお代替するＣＰＵは、情報処理装置１００のＣＰＵを使用してもよいが、アナログ計算機１３００の中あるいはこれに付随して、専用のＣＰＵあるいは論理回路を備えてもよい。アナログ計算機１３００の中に配置する場合には、複数のアナログ演算回路１３１１に対して１つのＣＰＵを備えるようにすれば、回路規模の縮小効果を得ることができる。
　専用のＣＰＵはソフトウェアによりプログラム可能なものでもよいし、FPGA（Field Programmable Gate Array）、ASIC（Application Specific Integrated Circuit）などのハードウェアでも実現できる。また、専用のハードウェアで構成してもよい。専用のハードウェアとしては、例えば図３に示した相互作用回路３１０相当の回路を、デジタル演算回路として使用できる。アナログ演算回路１３１１に誤動作可能性のある場合には、このデジタル演算回路でアナログ演算回路１３１１を代替させればよい。デジタル演算回路は、アナログ計算機１３００の外部あるいは内部に設けることができる。複数のアナログ演算回路１３１１に対して１つのデジタル演算回路を備えるようにすれば、回路規模の縮小効果を得ることができる。
　ステップ１６０６：ステップ１６０５の演算結果を、対応するアナログ計算ユニット１３１０の記憶装置１３１２に書き戻す。
　ステップ１６０７：本フローで行うべき所定の動作回数を満たしたかを確認する。満たしていなければ、ステップ１６０２に戻って繰り返す。
　ステップ１６０８：所定の動作回数を満たしたら、本フローを終了する。
　この動作フローによって、アナログ計算機ユニット１３１０が備えるアナログ演算回路１３１１において、製造ばらつきに起因する誤動作の可能性がある時、その時に限ってアナログ演算回路１３１１の役割を情報処理装置１００のＣＰＵ１１０が一時的に代替することが出来るようになり、情報処理装置１００全体として誤動作の無い処理が可能となる。
　なお、アナログ計算機ユニット１３１０の代替を担うＣＰＵ１１０は、アナログ計算機１３００の中に備えても良い。
　図１７は検査回路１３１３の構成の他の一例を説明する図である。検査回路１３１３を図１７に示す検査回路１７００のように実装しても良い。この検査回路１７００では、二種類の遅延回路（第１遅延回路、第２遅延回路）で生成した異なるタイミングのクロック信号でＯＵＴＡ信号（ＸＯＲからの出力）をフリップフロップに取り込む。例えば、第１遅延回路が時刻ｔにおいてクロックが立ち上がるクロック信号を生成し、第２遅延回路が時刻ｔ＋Δｔにおいてクロックが立ち上がるクロック信号を生成するものとする。
　そして、時刻ｔ＋Δｔ以降に２個のフリップフロップの出力（Ｌ１，Ｌ２）を観測すると、（Ｌ１，Ｌ２）＝（０，０），（Ｌ１，Ｌ２）＝（０，１），（Ｌ１，Ｌ２）＝（１，１）の３通りがある。ここで、（０，０）のケースが最もＣＴ側／ＣＢ側の電位差が小さく、逆に（１，１）はＣＴ側／ＣＢ側の電位差が大きい。よって、Ｌ１とＬ２のＸＯＲ（２入力が不一致の時に１を出力）を取ると（０，１）のケースが抽出される（ＸＯＲの出力に１が出る）が、このケースは電位差が（０，０）と（１，１）の中間程度にあると言える。
　よって、第１遅延回路の時刻ｔと、第２遅延回路の時刻ｔ＋Δｔを適切に設定することで、電圧差の大小を検出することが出来る。これにより、検査回路１３１３では電圧差が小さい場合全てを検出していたが、検査回路１７００では電圧差が所定より小さい場合、及び、所定より大きい場合は問題が無いと判断し、ある範囲内の電位差だった場合に問題があると検出することができる。
　この検査回路１７００を用いると、実施例１と実施例２を組み合わせた実施の形態をとることもできる。つまり、検査回路１７００で（Ｌ１，Ｌ２）＝（０，０）が拮抗状態の時に現れるように時刻ｔとｔ＋Δｔを設定することで、拮抗状態（０，０）、誤判定の可能性がある状態（０，１）、誤判定の可能性が無い状態（１，１）の３状態を分けて検出することが可能で、イジングモデル計算機２００に入力する問題は実施例１と同様にして生成しておき、拮抗状態に対する対処は実施例１と同様に行い、拮抗状態を超える電位差であり、かつ、その電位差が誤判定が無いと言いきれない程度に小さい場合には、実施例２記載の動作が行われるようにする。
　なお、本実施例も実施例１と同様にイジングモデルと相似な構造を持つイジングモデル計算機２００を例として実施例を示したが、ニューラルネットワークと相似な構造を持つニューロ計算機についても同様に実施することができることは、実施例１と同様である。
　本発明は上記した実施形態に限定されるものではなく、様々な変形例が含まれる。例えば、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることが可能である。また、各実施例の構成の一部について、他の実施例の構成の追加・削除・置換をすることが可能である。
　イジングモデルやニューラルネットワーク等の計算を行う技術分野に適用が可能である。
１００：情報処理装置
２００：イジングモデル計算機
３００：スピンユニット
１３００：アナログ計算機

Claims (15)

  1. 　処理の対象もしくは処理の結果となる数値を格納する数値格納メモリと、
    　前記処理に用いる係数を格納する係数格納メモリと、
    　前記数値と前記係数を用いて前記処理の少なくとも一部をアナログ的に実行するアナログ処理回路と、
    　を複数組備え、
    　前記複数のアナログ処理回路の特性ばらつきを示す特性データを格納する特性データ格納メモリを備え、
    　問題データを前記数値または前記係数として前記数値格納メモリと前記係数格納メモリに格納し、前記アナログ処理回路による１または複数回の処理の後、前記数値格納メモリに格納された数値を解データとして読み出す、
    　情報処理装置。
  2. 　前記問題データを格納する前に、前記特性データに基づいて前記問題データを変換し、
    　前記解データを読み出した後に、前記特性データに基づいて前記解データを変換する、
    　請求項１記載の情報処理装置。
  3. 　第１の数値格納メモリと、第１の係数格納メモリと、第１のアナログ処理回路と、
    　第２の数値格納メモリと、第２の係数格納メモリと、第２のアナログ処理回路と、
    を有する計算ユニットを備え、
    　前記第２のアナログ処理回路が、前記第１の数値格納メモリに格納された数値と、前記第２の係数格納メモリに格納された係数を用いて処理を行い、前記第２の数値格納メモリに処理の結果を格納する、
    　請求項２記載の情報処理装置。
  4. 　処理装置、記憶装置、入力装置、出力装置を備え、
    　前記入力装置が、前記問題データを受付け、
    　前記処理装置が、前記受付けた問題データを前記記憶装置に格納し、
    　前記処理装置が、前記記憶装置に格納されている前記問題データを読み出し、前記特性データに基づいて前記読み出した問題データを変換し、前記計算ユニットの数値格納メモリと前記係数格納メモリに格納し、
    　前記処理装置が、前記アナログ処理回路による複数回の処理の後、前記計算ユニットの数値格納メモリから前記解データを読み出し、前記特性データに基づいて前記読み出した解データを変換し、前記記憶装置または前記出力装置に出力する、
    　請求項３記載の情報処理装置。
  5. 　前記処理装置の制御により、前記複数のアナログ処理回路の特性ばらつきを測定し、測定結果を前記特性データとして前記特性データ格納メモリに格納する、
    　請求項４記載の情報処理装置。
  6. 　前記複数のアナログ処理回路の処理内容を検査する検査回路と、
    　前記アナログ処理回路による処理と同等の処理をデジタル的に実行するデジタル処理回路と、
    　を有し、
    　前記検査回路の検査結果に基づいて、前記アナログ処理回路の処理を前記デジタル処理回路に代替させる、
    　請求項１記載の情報処理装置。
  7. 　前記検査回路による検査は、前記アナログ処理回路の処理において生成される２つのアナログ値を検査するものであり、
    　前記アナログ値の差が第１の閾値より小さく、第２の閾値より大きい場合に、前記アナログ処理回路の処理を前記デジタル処理回路に代替させる、
    　請求項６記載の情報処理装置。
  8. 　前記問題データは、
    　イジングモデルのスピン配列を求めるための問題データ、あるいは、ニューラルネットワークのニューロンの状態を求めるための問題データである、
    　請求項１記載の情報処理装置。
  9. 　処理の対象もしくは処理の結果となる数値を格納する数値格納メモリと、
    　前記処理に用いる係数を格納する係数格納メモリと、
    　前記数値と前記係数を用いて少なくとも一部の処理をアナログ的に実行するアナログ処理回路と、
    　を複数組備え、
    　前記複数のアナログ処理回路の処理内容を検査する検査回路と、
    　前記アナログ処理回路による処理と同等の処理をデジタル的に実行するデジタル処理回路と、
    　を有し、
    　前記検査回路の検査結果に基づいて、前記アナログ処理回路の処理を前記デジタル処理回路に代替させる、
    　情報処理装置。
  10. 　前記検査回路による検査は、前記アナログ処理回路の処理において生成される２つのアナログ値の差を検査するものであり、
    　前記アナログ値の差が所定の条件下にある場合、前記アナログ処理回路の処理を前記デジタル処理回路に代替させる、
    　請求項９記載の情報処理装置。
  11. 　１個以上のアナログ計算機を含む情報処理装置を用いた情報処理方法において、
    　前記アナログ計算機の特性ばらつきを計測して特性ばらつきデータとして記録するステップ、
    　を有することを特徴とする情報処理方法。
  12. 　前記特性ばらつきデータを参照して、前記アナログ計算機に入力するデータを変換するステップと、
    　前記特性ばらつきデータを参照して、前記アナログ計算機から出力されたデータを解釈するステップと、
    　を有することを特徴とする請求項１１記載の情報処理方法。
  13. 　前記アナログ計算機はイジングモデルのスピン配列を求めるものであり、
    　前記入力するデータの変換は、イジングモデルの相互作用係数乃至は外部磁場係数の値の符号を反転させるものであり、
    　前記出力されたデータの解釈は、スピン値の反転であることを特徴とする請求項１２記載の情報処理方法。
  14. 　前記アナログ計算機は特性ばらつきの影響を評価する構成単位毎に、その評価を行う検査部を有しており、
    　前記検査部からの検査結果を読み出すことで、誤動作の可能性のある前記構成単位を検出することを特徴とする請求項１１記載の情報処理方法。
  15. 　前記誤動作の可能性のある構成単位を検出したときに、前記情報処理装置が有する論理回路が当該構成単位の担う処理を一時的に代替することを特徴とする請求項１４記載の情報処理方法。
PCT/JP2015/074263 2015-08-27 2015-08-27 情報処理方法及び装置 WO2017033337A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/JP2015/074263 WO2017033337A1 (ja) 2015-08-27 2015-08-27 情報処理方法及び装置
JP2017536156A JP6636529B2 (ja) 2015-08-27 2015-08-27 情報処理方法及び装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2015/074263 WO2017033337A1 (ja) 2015-08-27 2015-08-27 情報処理方法及び装置

Publications (1)

Publication Number Publication Date
WO2017033337A1 true WO2017033337A1 (ja) 2017-03-02

Family

ID=58101177

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2015/074263 WO2017033337A1 (ja) 2015-08-27 2015-08-27 情報処理方法及び装置

Country Status (2)

Country Link
JP (1) JP6636529B2 (ja)
WO (1) WO2017033337A1 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02232786A (ja) * 1989-03-07 1990-09-14 Fujitsu Ltd ニューロコンピュータにおける可変積分パルスによる誤差吸収方式
JPH064691A (ja) * 1992-06-23 1994-01-14 Nippon Telegr & Teleph Corp <Ntt> アナログニューラルネットワーク素子の試験方法
JPH0626858A (ja) * 1992-07-10 1994-02-04 Olympus Optical Co Ltd 測距装置
WO2014192153A1 (ja) * 2013-05-31 2014-12-04 株式会社日立製作所 半導体装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02232786A (ja) * 1989-03-07 1990-09-14 Fujitsu Ltd ニューロコンピュータにおける可変積分パルスによる誤差吸収方式
JPH064691A (ja) * 1992-06-23 1994-01-14 Nippon Telegr & Teleph Corp <Ntt> アナログニューラルネットワーク素子の試験方法
JPH0626858A (ja) * 1992-07-10 1994-02-04 Olympus Optical Co Ltd 測距装置
WO2014192153A1 (ja) * 2013-05-31 2014-12-04 株式会社日立製作所 半導体装置

Also Published As

Publication number Publication date
JPWO2017033337A1 (ja) 2018-08-02
JP6636529B2 (ja) 2020-01-29

Similar Documents

Publication Publication Date Title
CN109426698B (zh) 预测半导体集成电路良率的装置和半导体器件的制造方法
JP6006411B2 (ja) 半導体装置
US10796068B2 (en) Standard cell design system, standard cell design optimization method thereof, and semiconductor design system
Niki et al. A digitized replica bitline delay technique for random-variation-tolerant timing generation of SRAM sense amplifiers
JP6295325B2 (ja) 半導体装置及び情報処理システム
JP2010160787A (ja) パラメータ情報作成システム、歩留まり算出システム、プログラム及び記録媒体
US20160064053A1 (en) Semiconductor device and information processing device
Agbo et al. Quantification of sense amplifier offset voltage degradation due to zero-and run-time variability
US8813006B1 (en) Accelerated characterization of circuits for within-die process variations
US9666252B2 (en) Semiconductor device for calculating an interaction model
Agbo et al. BTI impact on SRAM sense amplifier
TW201928986A (zh) 測試系統及其操作方法
Amrouch et al. Special session: Machine learning for semiconductor test and reliability
Agbo et al. Comparative BTI analysis for various sense amplifier designs
Agbo et al. BTI analysis of SRAM write driver
WO2016030938A1 (ja) シミュレーションシステム、及びシミュレーション方法
WO2017033337A1 (ja) 情報処理方法及び装置
JP2008066598A (ja) ソフトエラー率の計算方法、プログラム、半導体集積回路の設計方法及び設計装置、並びに半導体集積回路
Tsiligiannis et al. SRAM soft error rate evaluation under atmospheric neutron radiation and PVT variations
Cheng et al. A new march test for process-variation induced delay faults in srams
US8988957B2 (en) Sense amplifier soft-fail detection circuit
US20160064080A1 (en) Semiconductor device and information processing device
Ghosh et al. Bridge-Defect Prediction in SRAM Circuits Using Random Forest, XGBoost, and LightGBM Learners
Selvam et al. Low power and low voltage SRAM design for LDPC codes hardware applications
Barceló et al. An SET propagation EDA tool based on analytical glitch propagation model

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15902306

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2017536156

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 15902306

Country of ref document: EP

Kind code of ref document: A1