JP2020046784A - 計算装置、計算プログラム、記録媒体及び計算方法 - Google Patents

Abstract

【課題】最適化問題を高速に計算できる計算装置、計算プログラム、記録媒体及び計算方法を提供する。【解決手段】実施形態によれば、計算装置は、処理手順を繰り返す処理部を含む。前記処理手順は、第１〜第３変数更新を含む。第１変数更新は、第１関数を加えてｉ番目の第１変数ｘｉを更新する。第２変数更新は、ｉ番目の第２関数、ｉ番目の第３関数及びｉ番目の第１要素関数の演算結果を加えて、ｉ番目の第２変数ｙｉを更新する。第３変数更新は、ｉ番目の第２要素関数を加えて前記第３変数ｚを更新する。第２要素関数の変数は、第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含む。処理部は、処理手順を繰り返した後に得られる前記ｉ番目の前記第１変数ｘｉを出力する出力する。【選択図】図１

Description

本発明の実施形態は、計算装置、計算プログラム、記録媒体及び計算方法に関する。

様々な社会課題において、最適化問題が生じる。最適化問題の１つの例として、イジング問題がある。大規模な最適化問題を高速に解くことが求められる。

T. Inagaki et al., Science 354, 603 (2016). H. Goto, Sci. Rep. 6, 21686 (2016). Y. Haribara et al., Quantum Sci. Technol. 2, 044002 (2017).

本発明の実施形態は、最適化問題を高速に計算できる計算装置、計算プログラム、記録媒体及び計算方法を提供する。

本発明の実施形態によれば、計算装置は、処理手順を繰り返す処理部を含む。前記処理手順は、第１変数更新、第２変数更新及び第３変数更新を含む。前記第１変数更新は、前記第１変数更新前の前記ｉ番目（ｉは１以上Ｎ以下の整数であり、前記Ｎは２以上の整数）の第１変数ｘ_ｉにｉ番目の第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することを含む。前記ｉ番目の前記第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つである。前記ｉ番目の前記第１関数の変数は、前記ｉ番目の第２変数ｙ_ｉ含む。前記ｉ番目の前記第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つである。前記ｉ番目の前記第１関数は、第１関数群の１つである。第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉに、前記ｉ番目の第２関数、前記ｉ番目の第３関数及び前記ｉ番目の第１要素関数の演算結果を加えて、前記ｉ番目の前記第２変数ｙ_ｉを更新することを含む。前記ｉ番目の前記第２関数の変数は、前記ｉ番目の前記第１変数ｘ_ｉを含む。前記ｉ番目の前記第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び前記第１変数群｛ｘ｝の少なくとも一部を含む。前記ｉ番目の前記第１要素関数の変数は、前記ｉ番目の前記第２変数ｙ_ｉと第３変数ｚとの積を含む。前記ｉ番目の前記第２関数は、第２関数群の１つである。前記ｉ番目の前記第３関数は、第３関数群の１つである。前記ｉ番目の前記第１要素関数は、第１要素関数群の１つである。前記第３変数更新は、前記第３変数更新前の前記第３変数ｚに前記ｉ番目の第２要素関数を加えて前記第３変数ｚを更新することを含む。前記第３変数ｚは、第３変数群｛ｚ｝の１つである。前記第２要素関数の変数は、前記第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含む。前記処理部は、少なくとも、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉ、及び、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉの関数の少なくともいずれかを出力する出力する。

実施形態に係る計算装置の例を示す模式図である。 実施形態に係る計算装置の動作を例示するフローチャートである。 実施形態に係る計算装置の動作の一部を例示するフローチャートである。 実施形態に係る計算装置の動作の一部を例示するフローチャートである。 実施形態に係る計算装置の動作の一部を例示するフローチャートである。 実施形態に係る計算装置の動作の一部を例示するフローチャートである。 実施形態に係る計算装置の動作の一部を例示するフローチャートである。 実施形態に係る計算装置の動作の一部を例示するフローチャートである。 実施形態に係る計算装置の動作の一部を例示するフローチャートである。 実施形態に係る計算装置の動作の一部を例示するフローチャートである。 実施形態に係る計算装置の例を示す模式図である。 実施形態に係る計算の特性を例示するグラフである。 計算結果を例示するグラフである。 図１４（ａ）及び図１４（ｂ）は、計算結果を例示するグラフである。

以下に、本発明の各実施の形態について図面を参照しつつ説明する。
本願明細書と各図において、既出の図に関して前述したものと同様の要素には同一の符号を付して詳細な説明は適宜省略する。

（第１実施形態）
図１は、実施形態に係る計算装置の例を示す模式図である。
図１に示すように、本実施形態に係る計算装置１１０は、例えば、処理部２０及び保持部１０を含む。処理部２０は、例えば、ＣＰＵ（Central Processing Unit）などを含む。処理部２０は、例えば電子回路などを含む。保持部１０は、種々のデータを保持可能である。保持部１０は、例えば、記憶部である。保持部１０は、ＲＯＭ（Read Only Memory）及びＲＡＭ（Random Access Memory）の少なくともいずれかを含んでも良い。計算装置１１０は、計算システムでも良い。

この例では、計算装置１１０に、取得部３１が設けられている。取得部３１は、例えば、種々のデータを入手可能である。取得部３１は、例えば、Ｉ／Ｏポートなどを含む。取得部３１は、出力部の機能を有しても良い。取得部３１は、例えば、通信機能を有しても良い。

図１に示すように、計算装置１１０は、操作部３２及び表示部３３などを含んでも良い。操作部３２は、例えば、操作機能を有する装置（例えばキーボード、マウス、タッチ式入力パネル、または音声認識入力装置など）を含んでも良い。表示部３３は、各種のディスプレイを含んでも良い。

計算装置１１０に含まれる複数の要素において、無線及び有線の少なくともいずれかの方法により、互いに通信が可能である。計算装置１１０に含まれる複数の要素が設けられる場所が、互いに異なっても良い。計算装置１１０として、例えば、汎用コンピュータが用いられても良い。計算装置１１０として、例えば、互いに接続された複数のコンピュータが用いられても良い。計算装置１１０の少なくとも一部（例えば、処理部２０及び保持部１０など）として、専用の回路が用いられても良い。計算装置１１０として、例えば、互いに接続された複数の回路が用いられても良い。

計算装置１１０に含まれる複数の要素の例については、後述する。

以下、実施形態に係る計算装置１１０において、実施される動作の例について説明する。

図２は、実施形態に係る計算装置の動作を例示するフローチャートである。
図２に示すように、パラメータを設定する（ステップＳ２０１）。パラメータは、例えば、第１パラメータ群｛Ｊ｝を含む。パラメータは、例えば、第２パラメータ群｛ｈ｝をさらに含んでも良い。これらのパラメータの例については、後述する。

複数の変数を設定する（ステップＳ２０２）。変数は、例えば、第１変数群｛ｘ｝及び第２変数群｛ｙ｝を含む。変数は、例えば、第３変数群｛ｚ｝をさらに含んでも良い。ステップＳ２０２において、変数は適切な値に初期化される。初期化の例については、後述する。

複数の変数の計算（例えば更新）を行う（ステップＳ２１０）。例えば、複数の変数の時間発展が計算される。例えば、第１変数群｛ｘ｝が更新され、第２変数群｛ｙ｝が更新される。第３変数群｛ｚ｝がさらに更新されても良い。これらの計算は、所定の条件（後述する）が満たされるまで繰り返される。ステップＳ２１０は、例えば、サブルーチンである。

サブルーチン（変数の更新）の後に、例えば、関数を算出する（ステップＳ２２０）。例えば、第１変数群｛ｘ｝に含まれる少なくとも１つの要素の関数が算出される。１つの例において、この関数は、第１変数群｛ｘ｝に含まれる少なくとも１つの要素の符号である。

この関数を出力する（ステップＳ２３０）。例えば、１つの例において、第１変数群｛ｘ｝に含まれる少なくとも１つの要素の符号が出力される。ステップＳ２３０において、更新後の第１変数群｛ｘ｝に含まれる少なくとも１つが出力されても良い。この場合、ステップＳ２２０は、省略されても良い。

実施形態においては、上記の複数の変数の計算（ステップＳ２１０）において、例えば、第１変数群｛ｘ｝の更新が、第２変数群｛ｙ｝を用いて行われる。そして、第２変数群｛ｙ｝の更新が、第１変数群｛ｘ｝を用いて行われる。これらの更新が、複数回行われる。１つの例において、複数回の更新の１つにおいて、第１変数群｛ｘ｝の更新の後に、第２変数群｛ｙ｝の更新が行われる。別の１つの例において、例えば、複数回の更新の１つにおいて、第２変数群｛ｙ｝の更新の後に、第１変数群｛ｘ｝の更新が行われる。第３変数群｛ｚ｝の更新が、さらに行われても良い。

計算装置１１０により、最適化問題を高速に計算できる。最適化問題は、例えば、組み合わせ最適化問題（例えば離散最適化問題）である。例えば、大規模なイジング問題を高速に解くことができる。

以下、計算装置１１０で実施される計算の例として、イジング問題について、説明する。

例えば、イジングエネルギーＥ_{Ｉｓｉｎｇ}は、以下の第１式で表される。

上記の第１式において、「Ｎ」はイジングスピンの数である。「ｓ_ｉ」は、ｉ番目のイジングスピンである。例えば、「ｓ_ｉ」＝±１である。「Ｊ」は、例えば、１つの行列である。上記の第１パラメータ群｛Ｊ｝の１つの例が、行列Ｊである。行列Ｊは、実対称行列である。実対称行列において、対角成分（対角要素）が全てゼロである。

上記の第１式に関して、量子分岐マシンの古典モデル（以下、古典分岐マシンと呼ぶ）が提案されている。古典分岐マシンにおいては、運動方程式は、以下の第２〜第４式で与えられる。

第２〜第４式において、「Ｎ」は、例えば、イジングスピンの数に対応する。「Ｄ」は、例えば、「離調」に対応する。「ｃ」は、定数である。「ｐ」は、例えば、「ポンプレート」に対応する。「Ｋ」は、例えば、「カー係数」に対応する。これらの値は、例えば、予め設定されても良い。第２〜第４式において、第２パラメータ群｛ｈ｝は設けられなくても良い。その場合は、第３式および第４式の中の｛ｈ｝の要素を含む項は無視される。

上記の第２〜第４式において、ｐ（ｔ）をゼロから十分大きな値へ増加させた時の「ｘ_ｉ」の最終値の符号「±１」が、最適解（基底状態）のイジングスピン「ｓ_ｉ」となる。「ａ（ｔ）」は、「ｐ（ｔ）」とともに増加するパラメータである。「ａ（ｔ）」は、例えば、以下の第５式で表される。

上記の古典分岐マシンは、第２〜第４式において、「Ｈ」をハミルトニアンとしたハミルトン力学系であると見なすことができる。

一方、シミュレーティッドアニーリングが知られている。この方法では、逐次更新アルゴリズムが採用される。逐次更新アルゴリズムにおいては、複数のスピンが、１つずつ更新される。このような逐次更新アルゴリズムは並列計算に向かない。

これに対して、上記の古典分岐マシンの運動方程式をデジタル計算機で離散解法によって解くことが考えられる。このアルゴリズムは、シミュレーティッドアニーリングとは異なり、並列更新アルゴリズムである。並列更新アルゴリズムにおいては、複数の変数を同時に更新することができる。このため、並列計算による高速化が期待できる。

上記の第２〜第５式を用いるアプローチには、以下の課題があると考えられる。最も計算量が大きい行列Ｊを用いた計算が、第１変数ｘ及び第２変数ｙの両方の更新に必要となる。上記の運動方程式は、数値的に容易には解けないため、例えば、計算量が大きい離散解法（例えば４次のルンゲ・クッタ法など）が必要となる。

これに対して、実施形態の１つの例においては、第２〜第４式に示した連立常微分方程式ではなく、例えば、以下の第６〜第８式に示す連立常微分方程式を用いる。

上記の「ｙ_Ｂ」は、第２変数群｛ｙ｝に含まれる全ての要素（例えば、第１番目の第２変数ｙ_１〜第Ｎ番目の第２変数ｙ_Ｎ）に対応する。「ｆ_ｉ（ｙ_Ｂ，ｔ）」は、ｉ番目の第２関数ｙ_ｉの関数である。

第７式に示すように、第１変数更新は、第１変数更新前のｉ番目の第１変数ｘ_ｉにｉ番目の第１関数を加えてｉ番目の第１変数ｘ_ｉを更新する。ｉ番目の第１関数は、「ｆ_ｉ（ｙ_Ｂ，ｔ）」である。ｉ番目の第１関数は、第１関数群の１つである。

第１計算方法（第１計算装置）において、「ｆ_ｉ（ｙ_Ｂ，ｔ）」として、以下の第９式が用いられる。このとき、第７式は、第１０式で表される。

従って、第１計算方法（第１計算装置）においては、第１変数ｘの更新は、上記の第１０式に基づいて行われる。そして、第２変数ｙの更新は、上記の第８式に基づいて行われる。第１０式に示すように、第１計算方法においては、ｉ番目の第１関数の変数は、パラメータ「Ｄ」と「ｙ_ｉ」の積を含む。パラメータ「Ｄ」は１でも良い。

第１計算方法（第１計算装置）において、処理手順が繰り返される。処理手順は、第１変数更新及び第２変数更新を含む。

実施形態の別の例においては、以下の第１１式〜第１３式が用いられる。

第２計算方法（第２計算装置）においても、処理手順が繰り返される。第２計算方法（第２計算装置）においては、処理手順は、第１変数更新、第２変数更新及び第３変数更新を含む。

第２計算方法（第２計算装置）において、第１変数更新は、第１１式（または第１０式でも良い）に基づいて行われる。例えば、第１変数更新は、第１変数更新前のｉ番目の第１変数ｘ_ｉにｉ番目の第１関数を加えてｉ番目の第１変数ｘ_ｉを更新することを含む。「ｉ」は、ｉは１以上Ｎ以下の整数である。「Ｎ」は、２以上の整数である。ｉ番目の第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つである。ｉ番目の第１関数の変数は、ｉ番目の第２変数ｙ_ｉ含む。ｉ番目の第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つである。ｉ番目の第１関数は、第１関数群の１つである。

第２計算方法（第２計算装置）において、ｉ番目の第１関数は、「ｆ_ｉ（ｙ_Ｂ，ｔ）」である。第２計算方法（第２計算装置）において、ｉ番目の第１関数は、「Ｄｙ_ｉ」でも良い（第１０式）。第２計算方法（第２計算装置）において、ｉ番目の第１関数は、「ｙ_ｉ」でも良い（第１１式）。

第１２式に示すように、第２変数更新は、第２変数更新前のｉ番目の第２変数ｙ_ｉに、ｉ番目の第２関数、ｉ番目の第３関数及びｉ番目の第１要素関数の演算結果を加えて、ｉ番目の第２変数ｙ_ｉを更新することを含む。

ｉ番目の第２関数は、第１２式の右辺の第１項及び第２項に対応する。ｉ番目の第３関数は、第１２式の右辺の第３項に対応する。ｉ番目の第１要素関数は、第１２式の第４項に対応する。

ｉ番目の第２関数の変数は、ｉ番目の第１変数ｘ_ｉを含む。ｉ番目の第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び第１変数群｛ｘ｝の少なくとも一部を含む。ｉ番目の第１要素関数の変数は、ｉ番目の第２変数ｙ_ｉと第３変数ｚとの積を含む。ｉ番目の第２関数は、第２関数群の１つである。ｉ番目の第３関数は、第３関数群の１つである。ｉ番目の第１要素関数は、第１要素関数群の１つである。

第１３式に示すように、第３変数更新は、第３変数更新前の第３変数ｚにｉ番目の第２要素関数を加えて第３変数ｚを更新することを含む。第２要素関数の変数は、第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータＴ（ｔ）と、の差を含む。この２乗平均は、第２変数群｛ｙ｝の上記の少なくとも一部の分散に比例する。第３変数ｚは、第３変数群｛ｚ｝の１つである。

例えば、１〜Ｎの集合が、共通部分を持たない複数の部分集合に分割される。複数の部分集合の数は、「Ｎ_ｓ」である。複数の部分集合のそれぞれに含まれる要素の数は、「Ｎ_ｌ」個である。例えば、上記の第３変数群｛ｚ｝として、ｚ_ｌ（ｌ＝１〜Ｎ_ｓ）が用いられる。第１２式の「ｉ」がｌ番目の部分集合「Ｓ_ｌ」に含まれる場合、第１２式の右辺最終項の第３変数ｚとして、「ｚ_ｌ」が用いられる。「ｚ_ｌ」は、第１３式において第３変数ｚを「ｚ_ｌ」とした式により、更新される。このとき、第１３式の「ｉ」に関する和は、「ｉ」が「Ｓ_ｌ」に含まれる場合に行われる。「Ｎ_ｌ」は、「Ｓ_ｌ」の要素数（「Ｓ_ｌ」に含まれる要素の数）である。第１３式の「ｉ」に関する和は、「Ｎ_ｌ」個の第２変数ｙについての２乗和となる。第１３式において、第１３式の右辺第１項の分母の「Ｎ」は、「Ｎ_ｌ」に置き換えられる。このとき、第１３式の第１項は、複数の部分集合の１つである「Ｓ_ｌ」に含まれる番号に対応する第２変数ｙの２乗平均に対応する。

第１計算方法（第１計算装置）及び第２計算方法（第２計算装置）の例（例えば後述する図３の例）においては、「処理手順」は、ステップＳ１１０、ステップＳ１２０及びステップＳ１３０を含む。「処理手順」は、例えば、ステップＳ１１０、ステップＳ１２０及びステップＳ１３０における「ｉ」に関するループ処理を含む。「処理手順」は、複数回繰り返される。

例えば、第１計算パラメータＴ（ｔ）は、処理手順の繰り返しにおいて、更新される。図３の例では、第１計算パラメータＴ（ｔ）はステップＳ１３０において更新される。例えば、処理手順の後の第１計算パラメータＴ（ｔ）は、その処理手順の前の第１計算パラメータＴ（ｔ）とは異なる。例えば、処理手順の後の第１計算パラメータＴ（ｔ）は、その処理手順の前の第１計算パラメータＴ（ｔ）よりも小さい。

第６〜第１３式において、「Ｎ」は、例えば、イジングスピンの数に対応する。「Ｄ」は、例えば、「離調」に対応する。「ｃ」は、定数である。「ｐ」は、例えば、「ポンプレート」（例えば演算パラメータ）に対応する。「Ｋ」は、例えば、「カー係数」に対応する。これらの値は、例えば、予め設定されても良い。第６〜第１３式において、第２パラメータ群｛ｈ｝は設けられなくても良い。その場合は、第８式の中の｛ｈ｝の要素を含む項は、無視される。第１２式において、第８式と同様に、｛ｈ｝の要素を含む項が設けられても良い。

最も計算量が大きい行列Ｊに関する積和演算は、第２変数ｙの更新にのみ行われ、第１変数ｘの更新では行われない。従って、計算量が削減される。これらの式において、第１変数ｘの時間微分は、第２変数ｙを含む。例えば、第１変数ｘの時間微分は、第１変数ｘを含まない。第２変数ｙの時間微分は、第１変数ｘを含む。第１計算方法においては、例えば、第２変数ｙの時間微分は、第２変数ｙを含まない。ハミルトニアンにおいて、ｘ及びｙは、互いに分離されている。このため、計算量が小さく安定な離散解法が適用可能となる。例えば、シンプレクティック・オイラー法と呼ばれる方法が適用できる。上記の第６〜第１３式においては、「ｘ」の時間微分から「ｐ」が消去されている。

このような方法を用いた場合に、高い性能（例えば、高い精度）が維持できることが分かった。実施形態に係る計算装置では、上記の分離可能なハミルトニアンを持つハミルトン力学系（新しい古典分岐マシン）の運動方程式が、例えば、シンプレクティック・オイラー法で解かれる。実施形態に係る計算装置は、このような新しいアルゴリズムの計算が並列計算によって、できる限り高速に実行されるように構成される。

実施形態において、例えば、「ｐ（ｔ）」をゼロから十分大きな値へ増加させた時の第１変数ｘ_ｉの最終値の符号（「±１」）が、最適解（基底状態）のイジングスピンｓ_ｉとなる。

例えば、第１変数ｘ_ｉ及び第２変数ｙ_ｉは、変数の設定（ステップＳ２０２）において、適切な値に、初期化される。第３変数ｚが設けられる場合、第３変数ｚは、変数の設定（ステップＳ２０２）において、適切な値に、初期化される。例えば、これらの変数は、絶対値が０．１以下の乱数によって、ランダムに初期化される。

以下、ステップＳ２１０（図２参照）のいくつかの例について説明する。

図３は、実施形態に係る計算装置の動作の一部を例示するフローチャートである。
図３は、ステップＳ２１０を例示している。図３に示す複数の式は、第２計算方法に関する１つの例である。第１計算方法においては、第２計算方法に関する例の式から、第６式〜第１０式に応じて、図３に例示する式が適宜変更される。図３に示すように、「ｔ」、「ｐ」「ａ」及び「Ｔ」を初期化する（ステップＳ１０１）。１つの例において、「ｔ」、「ｐ」及び「ａ」は０とされる。「Ｔ」は、「Ｔ０」とされる。「Ｔ０」は、１つの初期値である。図３において、「ｔ_１」は、時刻ｔの最終値に対応する。「ｄｔ」は、時刻ｔの１ステップ当たりの増加量である。「ｄｐ」は、パラメータｐの１ステップ当たりの増加量である。

図３の例において、ステップＳ１０５とステップＳ１０６との間の処理が、１つのループに対応する。「ｔ」が「ｔ_１」よりも小さいときに、以下に説明する一連の処理を含む「処理手順」の繰り返しが行われる（ステップＳ１０５）。例えば、「ｐ」があらかじめ適切に設定された「Ｐ」よりも小さいときに「処理手順」の繰り返しが行われる、としても良い。例えば、「処理手順」は「ｉ」に関するループ処理を含む（ステップＳ３０１ａ及びステップＳ３０１ｂ）。

例えば、パラメータνが０とされる（ステップＳ４０１）。

１つのループ処理（ステップＳ３０１ａ及びステップＳ３０１ｂ）の中において、ｉ番目の第１変数ｘ_ｉの更新が行われる（ステップＳ１１０）。例えば、第１計算方法（第１計算装置）及び第２計算方法（第２計算装置）においては、更新前の第１変数ｘ_ｉに、ｄｔ＊Ｄ＊ｙ_ｉ（第１０式参照）を加えて得られる値を、更新後の第１変数ｘ_ｉとする。ここで、「＊」は、積の記号である。更新前の第１変数ｘ_ｉに、ｄｔ＊ｙ_ｉ（第１１式参照）を加えて得られる値を、更新後の第１変数ｘ_ｉとしても良い。

ｉ番目の第２変数ｙ_ｉの更新が行われる（ステップＳ１２０）。この例では、ステップＳ１２１及びステップＳ１２２とが行われる。

第１計算方法において、ステップＳ１２１は、第１変数群{ｘ}に基づく更新に対応する。一方、第２計算方法においては、ステップＳ１２１は、第１変数群{ｘ}及び第３変数群{ｚ}に基づく更新に対応する。第１計算方法及び第２計算方法において、ステップＳ１２２は、第１パラメータ{Ｊ}及び第１変数群{ｘ}に基づく更新に対応する。ステップＳ１２１とステップＳ１２２との順序は、入れ替えが可能である。

ステップＳ１２１の少なくとも一部と、ステップＳ１２２の少なくとも一部と、が同時に行われても良い。ステップＳ１２１は、例えば第１サブ更新に対応する。ステップＳ１２２は、第２サブ更新に対応する。ステップＳ１２１は、「ｉ」に関してループ処理される（ステップＳ３０１ａ及びステップＳ３０１ｂ）。ステップＳ１２２は、「ｉ」に関してループ処理される（ステップＳ３０２ａ及びステップＳ３０２ｂ）。

第１計算方法において、第１サブ更新では、例えば、更新前の第２変数ｙ_ｉに、ｄｔ＊［（ｐ−Ｄ−ｘ_ｉ＊ｘ_ｉ）＊ｘ_ｉ−ｃ＊ｈ_ｉ＊ａ］を加えて得られる値を、更新後の第２変数ｙ_ｉとする。

第１計算方法において、第２サブ更新では、例えば、更新前の第２変数ｙ_ｉに、ｄｔ＊ｃ＊Σ（Ｊ_ｉ，ｊ＊ｘ_ｊ）を加えて得られる値を更新後の第２変数ｙ_ｉとする。「Σ」は、ｊに関する和を表す。例えば、「ｄｔ＊ｃ＊Ｊ」をＪ行列としてもよい。この場合、「ｄｔ＊ｃ＊」の演算を実際に行わなくても良い。

一方、第２計算方法において、第１サブ更新では、例えば、更新前の第２変数ｙ_ｉに、ｄｔ＊［（ｐ−Ｄ−ｘ_ｉ＊ｘ_ｉ）＊ｘ_ｉ−ｃ＊ｈ_ｉ＊ａ−ｚ＊ｙ_ｉ］を加えて得られる値を、更新後の第２変数ｙ_ｉとする。第２計算方法おける第２サブ更新は、例えば、第１計算方法おける第２サブ更新と同様で良い。

第２計算方法においては、パラメータνは、１つのループ処理（ステップＳ３０１ａ及びステップＳ３０１ｂ）の中において更新される（ステップＳ４０２）。パラメータνの更新においては、例えば、更新前のパラメータνに、ｙ_ｉ＊ｙ_ｉを加えて得られ値を、更新後のパラメータνとする。

例えば、１つのループ処理（ステップＳ３０１ａ及びステップＳ３０１ｂ）が終了すると、第３変数群{ｚ}が更新される（第３変数更新）。

第３変数更新（ステップＳ４０３）は、第３変数更新前の第３変数ｚに、ｄｔ＊ｒ＊（ν／Ｎ−Ｔ）を加えて得られる値を、更新後の第３変数ｚとする。「ｄｔ＊ｒ＊（ν／Ｎ−Ｔ）」は、ｉ番目の第２要素関数の１つの例に対応する。第２要素関数の変数は、第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータＴ（ｔ）と、の差を含む。ステップＳ４０３における「Ｔ」は、第１計算パラメータＴ（ｔ）の１つの例である。

更新のためのパラメータの更新処理を行う（ステップＳ１３０）。すなわち、更新前の「ｔ」に「ｄｔ」を加えて得られる値を更新後の「ｔ」とする。更新前の「ｐ」に「ｄｐ」を加えて得られる値を更新後の「ｐ」とする。「ａ」は、例えば、ｐ^１／２である。更新前の「Ｔ」に「−ｄＴ」を加えて得られる値を更新後の「Ｔ」とする（第１計算パラメータＴ（ｔ）の更新）。「Ｔ」の更新は、第２計算方法において行われ、第１計算方法では行われない。

そして、「ｔ」が「ｔ_１」よりも小さいときに、ステップＳ１０５に戻る（ステップＳ１０６）。例えば、「ｐ」があらかじめ適切に設定された「Ｐ」よりも小さいときに、ステップＳ１０５に戻っても良い。

「ｔ」が「ｔ_１」以上のときに、更新を終了し、図２に示すステップＳ２２０またはステップＳ２３０に進む。

図４は、実施形態に係る計算装置の動作の一部を例示するフローチャートである。
図４は、ステップＳ２１０を例示している。図４に示す例では、１つの繰り返し処理において、ステップＳ１１０の前にステップＳ１２０（図３参照）が実施される。このように、ステップＳ１１０とステップＳ１２０とにおいて、順序は、任意である。

１つの例において、「Ｋ」は１とする。例えば、「Ｎ」、「Ｄ」、「ｃ」、「ｔ_１」、「ｄｔ」及び「ｄｐ」は、予め適切な値に設定することができる。「Ｔ０」は、予め適切な値に設定することができる。

図３及び図４の例では、「ｐ」の更新において、線形増加が適用される。実施形態において、「ｐ」の更新には、任意の増加関数が用いられても良い。実施形態において、上記のように、２種類の更新方法がある。すなわち、１つの更新方法では、第１変数ｘ_ｉの更新の後に、更新された第１変数ｘ_ｉを用いて第２変数ｙ_ｉを更新する。別の更新方法では、第２変数ｙ_ｉの更新の後に、更新された第２変数ｙ_ｉを用いて第１変数ｘ_ｉを更新する。これらの２つの方法が、図３及び図４にそれぞれ対応する。

このように、本実施形態に係る計算装置１１０において、処理部２０（図１参照）は、処理手順（ステップＳ２１０：図２参照）を繰り返す。第１計算方法においては、この処理手順は、例えば、第１変数更新（ステップＳ１１０）及び第２変数更新（ステップＳ１２１及びステップＳ１２２）を含む。第２計算方法においては、この処理手順は、例えば、第１変数更新（ステップＳ１１０）、第２変数更新（ステップＳ１２１及びステップＳ１２２）及び第３変数更新（ステップＳ４０３）を含む。

第１計算方法及び第２計算方法において、第１変数更新は、第１変数更新前のｉ番目（ｉは１以上Ｎ以下の整数であり、Ｎは２以上の整数）の第１変数ｘ_ｉにｉ番目の第１関数を加えてｉ番目の第１変数ｘ_ｉを更新することを含む。ｉ番目の第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つである。ｉ番目の第１関数の変数は、ｉ番目の第２変数ｙ_ｉを含む。ｉ番目の第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つである。ｉ番目の第１関数は、第１関数群の１つである。

第１計算方法において、第２変数更新は、第２変数更新前のｉ番目の第２変数ｙ_ｉにｉ番目の第２関数及びｉ番目の第３関数を加えて、ｉ番目の第２変数ｙ_ｉを更新することを含む。ｉ番目の第２関数の変数は、ｉ番目の第１変数ｘ_ｉを含む。ｉ番目の第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び第１変数群｛ｘ｝の少なくとも一部を含む。ｉ番目の第２関数は、第２関数群の１つである。ｉ番目の第３関数は、第３関数群の１つである。

一方、第２計算方法において、第２変数更新は、第２変数更新前のｉ番目の第２変数ｙ_ｉに、ｉ番目の第２関数、ｉ番目の第３関数及びｉ番目の第１要素関数の演算結果を加えてｉ番目の第２変数ｙ_ｉを更新することを含む。ｉ番目の第２関数の変数は、ｉ番目の第１変数ｘ_ｉを含む。ｉ番目の第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び第１変数群｛ｘ｝の少なくとも一部を含む。ｉ番目の第１要素関数の変数は、ｉ番目の第２変数ｙ_ｉと第３変数ｚとの積を含む。ｉ番目の第２関数は、第２関数群の１つである。ｉ番目の第３関数は、第３関数群の１つである。ｉ番目の第１要素関数は、第１要素関数群の１つである。

第２計算方法において、第３変数更新は、第３変数更新前の第３変数ｚにｉ番目の第２要素関数を加えて第３変数ｚを更新することを含む。第３変数ｚは、第３変数群｛ｚ｝の１つである。第２要素関数の変数は、第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータＴ（ｔ）と、の差を含む。

そして、処理部２０は、上記の処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉ、及び、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉの関数の少なくともいずれかを出力する。例えば、第１変数群｛ｘ｝に含まれる任意の第１変数、及び、処理手順を繰り返した後に得られる任意の第１変数の関数の少なくともいずれかが出力される。例えば、第１変数群｛ｘ｝に含まれる全ての第１変数、及び、処理手順を繰り返した後に得られる全ての第１変数のそれぞれの関数の少なくともいずれかが出力されても良い。

第１計算方法及び第２計算方法において、上記の第１関数（第１関数群）は、第１変数群｛ｘ｝から独立している。第１変数群｛ｘ｝の値を変更しても、第１関数（第１関数群）の結果の値は変化しない。第１計算方法において、上記の第２関数（第２関数群）は、第２変数群｛ｙ｝から独立している。第２変数群｛ｙ｝の値を変更しても、第２関数（第２関数群）の結果の値は変化しない。第１計算方法及び第２計算方法において、上記の第３関数（第３関数群）は、第２変数群｛ｙ｝から独立している。第２変数群｛ｙ｝の値を変更しても、第３関数（第３関数群）の結果の値は変化しない。

第１計算方法（第１計算装置）及び第２計算方法（第２計算装置）においては、ｉ番目の第１関数は、例えば、ｄｔ＊Ｄ＊ｙ_ｉである（例えば、第１０式参照）。ｉ番目の第１関数は、例えば、ｄｔ＊ｙ_ｉでも良い（例えば、第１１式参照）。

第１計算方法（第１計算装置）において、第２関数は、例えば、ｄｔ＊［（ｐ−Ｄ−ｘ_ｉ＊ｘ_ｉ）＊ｘ_ｉ−ｃ＊ｈ_ｉ＊ａ］である。第３関数は、例えば、ｄｔ＊ｃ＊Σ（Ｊ_ｉ，ｊ＊ｘ_ｊ）である。

第２計算方法（第２計算装置）において、第２関数は、例えば、ｄｔ＊［（ｐ−Ｄ−ｘ_ｉ＊ｘ_ｉ）＊ｘ_ｉ−ｃ＊ｈ_ｉ＊ａ］である。第３関数は、例えば、ｄｔ＊ｃ＊Σ（Ｊ_ｉ，ｊ＊ｘ_ｊ）である。第１要素関数は、例えば、ｄｔ＊［−ｚ＊ｙ_ｉ］である。第２計算方法（第２計算装置）において、第２要素関数は、例えば、「ｄｔ＊ｒ＊（ν／Ｎ−Ｔ）」である。

実施形態において、繰り返して行われる上記の処理手順の１つにおいて、第１変数更新、第２変数更新及び第３変数更新の順序は相互に入れ替えが可能である。

実施形態に係る計算装置で行われるアルゴリズムは、例えば、以下を含む。
例えば、行列Ｊ（第１パラメータ群｛Ｊ｝の１つの例）を取得する。または、行列Ｊを計算により定める。行列Ｊは、例えば、イジングモデルのパラメータである。このとき、ベクトルｈ（第２パラメータ群｛ｈ｝の１つの例）をさらに取得しても良い。または、ベクトルｈは、計算により定められても良い。

第１計算方法においては、２種類の変数（第１変数群{ｘ}及び第２変数群{ｙ}）を用いる。一方の変数の更新には、他方の変数の値を用いる。一方の変数の更新において、その一方の変数の値を用いない。一方の変数の更新の後に、更新後のその一方の変数の値を用いて、他方を更新する。

第２計算方法においては、３種類の変数（第１変数群{ｘ}、第２変数群{ｙ}及び第３変数群{ｚ}）を用いる。

上記の第２関数は、例えば、ｉ番目の第１変数ｘ_ｉの非線形関数である第４関数を含む。第４関数は演算パラメータ「ｐ」も含む。複数の種類の変数の更新とともに「ｐ」は変化する。

複数の種類の変数の更新とともに「ｐ」が変化すると、上記の第４関数の実数根の数が変化する。「関数の実数根」とは、関数の値がゼロとなる変数の値（実数）のことである。第２変数更新において第４関数のみを考慮した場合、第４関数の実数根は非線形力学系の固定点（fixed point）に対応する。（ハミルトン力学系では、固定点はハミルトニアンの極値に対応する。）よって、第４関数の実数根の数が変化することは、固定点の数が変化することに対応する。これは、非線形力学系の分岐現象に対応する。実施形態に係る計算装置で用いるアルゴリズムでは、変数の初期値が、初期の１つの安定固定点（stable fixed point）の近傍に設定される。「ｐ」を変化させることで分岐（bifurcation）を起こす。分岐後の複数の安定な固定点（変数の値はこの複数の安定固定点のうちの１つの近傍へと変化する）と、解きたい組合せ最適化の離散変数と、を対応させる。これにより、分岐現象によって、その組合せ最適化問題を解く。例えば、上記の例では、分岐後の安定固定点における各ｘの値が、正負の２値であり、その符号と、イジングスピン（イジング問題の離散変数）と、が、対応付けられている。初期の安定固定点が原点であるため、各ｘの初期値と各ｙの初期値と、を原点の近傍（つまり、絶対値が０．１以下の小さな乱数）の値に設定する。

第４関数は、例えば、ｄｔ＊（ｐ−Ｄ’−ｘｉ＊ｘｉ）＊ｘｉである。「Ｄ’」は、０≦Ｄ’≦Ｄを満たす適切な定数である。初期時刻ではｐ＝０であり、第４関数の根はｘｉ＝０の１つだけである、ｐがＤ’よりも大きくなると、根は３つとなり、そのうちの正負２つの根が、イジングスピンと対応付けられる。第２関数を、例えば、ｄｔ＊［（ｐ−Ｄ−ｘｉ＊ｘｉ）＊ｘｉ−ｃ＊ｈｉ＊ａ］とした場合、第２関数は、ｄｔ＊（ｐ−Ｄ’−ｘｉ＊ｘｉ）＊ｘｉ＋ｄｔ＊［−（Ｄ−Ｄ’）＊ｘｉ−ｃ＊ｈｉ＊ａ］のように、第４関数と１次関数との和で表すことができる。従って、第２関数は、第４関数を含む。

上記の第４関数は、例えば、３次関数である。このような処理により、例えば、ニューラルネットワークで用いられている非線形関数（例えばシグモイド関数）を用いた計算よりも、計算が容易になる。

実施形態においては、時間ステップ（例えば「ｄｔ」）を大きくすると、計算が高速になる。一方、時間ステップを過度に大きくすると、計算が不安定になる。このことを考慮して、計算の一部における時間ステップを大きくし、計算の別の一部では、時間ステップを小さくても良い。例えば、計算量の大きい、行列Ｊ及び第１変数ｘの積和演算を含む更新には、大きい時間ステップを適用しても良い。その他の更新には、小さい時間ステップを適用しても良い。これにより、さらなる高速化が可能である。

以下、このような計算を行う場合の例について説明する。
図５は、実施形態に係る計算装置の動作の一部を例示するフローチャートである。
図５は、ステップＳ２１０を例示している。図５に示す例では、１つの「処理手順」（ステップＳ１０５とステップＳ１０６との間にある処理のまとまり）の中に、小さいループ（ステップＳ１０７ａ〜ステップＳ１０７ｂ）が設けられている。ステップＳ１０７ａにおいて、ループ変数「ｍ」は、１以上Ｍ以下である。Ｍは、例えば、２以上の整数である。小さいループの中において、ステップＳ１１０及びステップＳ１２１がＭ回繰り返される。ステップＳ１１０とステップＳ１２１の順序は、入れ替えが可能である。この後、ステップＳ１２２に進む。

第１計算方法において、図５の例では、ステップＳ１１０及びステップＳ１２１の繰り返しの後に、ステップＳ１２２が実施される。第２計算方法において、図５の例では、ステップＳ１１０、ステップＳ１２１及びステップＳ４０３の繰り返しの後に、ステップＳ１２２が実施される。

図６は、実施形態に係る計算装置の動作の一部を例示するフローチャートである。
図６は、ステップＳ２１０を例示している。図６に示す例でも、１つのループ（ステップＳ１０５〜ステップＳ１０６）の中に、小さいループ（ステップＳ１０７ｃ〜ステップＳ１０７ｄ）が設けられている。ステップＳ１０７ｃにおいて、ループ変数「ｍ」は、１以上Ｍ以下である。ステップＳ１２２が実施された後に、小さいループが実施される。第１計算方法において、小さいループは、ステップＳ１２１及びステップＳ１１０の繰り返しに対応する。第２計算方法において、小さいループは、ステップＳ１２１、ステップＳ１１０及びステップＳ４０３の繰り返しに対応する。小さいループの中において、ステップＳ１２１及びステップＳ１１０、または、ステップＳ１２１、ステップＳ１１０及びステップＳ４０３が、Ｍ回繰り返される。ステップＳ１２１とステップＳ１１０の順序は、入れ替えが可能である。ステップＳ１２１、ステップＳ１１０及びステップＳ４０３の順序は、入れ替えが可能である。

図５の例では、第１変数ｘ_ｉの更新後に、第２変数ｙ_ｉの更新が行われる。図６の例では、第２変数ｙ_ｉの更新後に、第１変数ｘ_ｉの更新が行われる。例えば、行列Ｊに関する積和演算を含まない更新の時間ステップ「ｄｔ’」は、「ｄｔ／Ｍ」とされる。一方、行列Ｊに関する積和演算を含む１回の更新（大きいループ）において、Ｍ回の小さいループ（行列Ｊに関する積和演算を含まない更新）を行う。上記により、例えば、大きなループの時間ステップｄｔを比較的大きな値にすることができる。例えば、高速な計算が可能となる。

このように、実施形態の１つの例では、上記の第２変数更新（ステップＳ１２０）は、第１サブ更新（ステップＳ１２１）及び第２サブ更新（ステップＳ１２２）を含む。

第１サブ更新（ステップＳ１２１）は、第１サブ更新前のｉ番目の第２変数ｙ_ｉに第２関数を加えてｉ番目の第２変数ｙ_ｉを更新することを含む。第２サブ更新（ステップＳ１２２）は、第２サブ更新前のｉ番目の第２変数ｙ_ｉに第３関数を加えてｉ番目の第２変数ｙ_ｉを更新することを含む。この場合も、第２関数は、第２変数群｛ｙ｝から独立している。第３関数は、第２変数群｛ｙ｝から独立している。

第１計算方法において、例えば、第１変数更新及び第１サブ更新を交互にＭ回（Ｍは２以上の整数）行った後に、第２サブ更新を行う。または、第２サブ更新後に、第１変数更新及び第１サブ更新を交互にＭ回行う。第１変数更新と第１サブ更新とを交互に行う順序は、入れ替えが可能である。

第２計算方法において、例えば、第１変数更新、第１サブ更新及び第３変数更新をＭ回（Ｍは２以上の整数）行った後に、第２サブ更新を行う。または、第２サブ更新後に、第１変数更新、第１サブ更新及び第３変数更新をＭ回行う。第１変数更新、第１サブ更新及び第３変数更新の順序は、入れ替えが可能である。

第２計算方法において、第３変数更新（ステップＳ４０３）の少なくとも一部は、第２サブ更新（ステップＳ１２２）の少なくとも一部と同時に実施されても良い。

第３関数は、例えば、第１パラメータ群｛Ｊ｝の上記の少なくとも一部と、第１変数群｛ｘ｝の上記の少なくとも一部と、の積和演算を含む。

１つの例において、上記の処理手順を繰り返した後の第４関数の実数根の数は、２以上である。処理手順を繰り返した後の第４関数の根の１つは正である。処理手順を繰り返した後の第４関数の根の別の１つは負である。処理部２０（図１参照）は、例えば、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉの符号（すなわち±１）を出力する。

既に説明したように、第２関数は、ｉ番目の第２パラメータｈ_ｉを含んでも良い。ｉ番目の第２パラメータｈ_ｉは、第２パラメータ群｛ｈ｝の１つである。

実施形態において、処理部２０は、保持部１０に保持されたデータを読み出し、データを更新し、更新されたデータを保持部１０に保持させる。

例えば、第１変数更新（ステップＳ１１０）は、第１変数更新前のｉ番目の第１変数ｘ_ｉを保持部１０から取得し、第１変数更新後のｉ番目の第１変数ｘ_ｉを保持部１０に保持させることを含む。第２変数更新（ステップＳ１２０）は、第２変数更新前のｉ番目の第２変数ｙ_ｉを保持部１０から取得し、第２変数更新後のｉ番目の第２変数ｙ_ｉを保持部１０に保持させることを含む。第３変数更新（ステップＳ４０３）は、第３変数更新前の第３変数ｚを保持部１０から取得し、第３変数更新後の第３変数ｚを保持部１０に保持させることを含む。

例えば、第１変数更新は、ｉ番目の第２変数ｙ_ｉを保持部１０から取得してｉ番目の第１関数を計算し、ｉ番目の第１変数ｘ_ｉにｉ番目の第１関数を加えてｉ番目の第１変数ｘ_ｉを更新することをさらに含んでも良い。

例えば、第２変数更新は、ｉ番目の第１変数ｘ_ｉを保持部１０から取得して第２関数を計算し、第１パラメータ群｛Ｊ｝の上記の少なくとも一部及び第１変数群｛ｘ｝の上記の少なくとも一部を保持部１０から取得して第３関数を計算し、ｉ番目の第２変数ｙ_ｉに第２関数及び第３関数を加えてｉ番目の第２変数ｙ_ｉを更新することをさらに含んでも良い。

例えば、第２変数更新は、ｉ番目の第１変数ｘ_ｉを保持部１０から取得して第２関数を計算し、第１パラメータ群｛Ｊ｝の上記の少なくとも一部及び第１変数群｛ｘ｝の上記の少なくとも一部を保持部１０から取得して第３関数を計算し、ｉ番目の第２変数ｙ_ｉを保持部１０から取得して第１要素関数を計算することを含んでも良い。例えば、第２変数更新は、ｉ番目の第２変数ｙ_ｉに、第２関数、第３関数及び第１要素関数の演算結果を加えてｉ番目の第２変数ｙ_ｉを更新すること含んでも良い。

実施形態において、例えば、行列Ｊが疎行列の場合、疎行列の圧縮形式を用いても良い。疎行列の圧縮形式には、例えば、ＣＯＯ（coordinate）形式またはＣＳＲ（compressed sparse row）形式などが適用できる。疎行列の圧縮形式を用いることで、例えば、メモリサイズを節約できる。疎行列の圧縮形式を用いることで、例えば、行列Ｊ及び第１変数ｘの積和演算を高速に実施できる。

以下、定数「ｃ」の例について述べる。例えば、離調「Ｄ」は、行列Ｊの最大固有値λｍａｘのｃ倍よりも大きくされる（例えば、非特許文献２参照）。「Ｄ」が大きすぎると、無駄な計算時間が生じる。このため、例えば、「Ｄ」がλｍａｘのｃ倍と実質的に等しいように設定する。この場合、ｃ＝Ｄ／λｍａｘとなる。一方、１つの例において、行列Ｊは、実対称行列である。この場合に、行列Ｊのサイズが十分大きい場合、λｍａｘは、２σ×Ｎ^１／２と実質的に同じになる。この関係は、ランダム行列のウィグナーの半円則に基づく。「σ」は、行列Ｊの非対角成分の標準偏差である。この場合、ｃ＝Ｄ／（２σ×Ｎ^１／２）と設定するとよい。この場合の計算例について、後述する。

本実施形態において、精度を向上させる方法として、上記の非線形関数として用いられる関数を変更することが考えられる。例えば、第１計算方法及び第２計算方法において、上記の第９式〜第１３式に関して説明した関数を適宜用いることができる。さらに、第８式の代わりに、以下の第１４式を用いても良い。

第１４式において、「ｎ」は、２以上の偶数である。このような関数を用いることで、例えば、イジング問題の解の精度を向上することができる。

実施形態に係る計算装置１１０で実施される上記のアルゴリズムは、種々の構成により実施できる。計算装置１１０は、例えば、ＰＣクラスタを含んでも良い。計算装置１１０は、例えば、ＧＰＵ（Graphics Processing Unit）を含んでも良い。計算装置１１０は、例えば専用回路を含んでも良い。専用回路は、例えば、ＦＰＧＡ（field-programmable gate array）、ゲートアレイ、及び、ＡＳＩＣ（application specific integrated circuit）の少なくともいずれかを含んでも良い。計算装置１１０は、例えば、並列デジタル計算装置を含んでも良い。

図７〜図１０は、実施形態に係る計算装置の動作の一部を例示するフローチャートである。
これらの図は、ステップＳ２１０の別の例を示している。これらの図に示すように、図３〜図６の例において、ステップＳ１１０とステップＳ１２１とは、互いに入れ替えられても良い。

図１１は、実施形態に係る計算装置の例を示す模式図である。
図１１に示すように、実施形態に係る計算装置１１１は、複数の回路（第１回路１５Ａ、第２回路１５Ｂ及び第３回路１５Ｃなど）を含む。これらの複数の回路のそれぞれは、例えば、１つのコンピュータである。これらの複数の回路のそれぞれは、例えば、１つの半導体回路でも良い。これらの複数の回路は、互いに通信（例えば、データの送受信）が可能である。計算装置１１１には、さらに制御回路部１５Ｘが設けられている。制御回路部１５Ｘにより、複数の回路における通信が制御される。

例えば、複数の回路のそれぞれには、処理部及び保持部（記憶部）が設けられる。この他、制御部が設けられても良い。

複数の回路（第１回路１５Ａ、第２回路１５Ｂ及び第３回路１５Ｃなど）により、並列計算が行われる。複数の回路の数は、任意である。

例えば、第１回路１５Ａには、第１計算部２０Ａ及び第１保持領域１０Ａが設けられる。この例では、第１回路１５Ａは、第１制御部１６Ａをさらに含む。第２回路１５Ｂには、第２計算部２０Ｂ及び第２保持領域１０Ｂが設けられる。この例では、第２回路１５Ｂは、第２制御部１６Ｂをさらに含む。このような構成は、第３回路１５Ｃにも設けられる。

処理部２０は、上記の複数の計算部（第１計算部２０Ａ及び第２計算部２０Ｂなど）を含む。例えば、第１計算部２０Ａは、第３関数の計算の一部（例えば、ステップＳ３０２ａ、ステップＳ１２２及びステップＳ３０２ｂのｉに関するループの一部）を実施する。第２計算部２０Ｂは、第３関数の計算の別の一部（例えば、ステップＳ３０２ａ、ステップＳ１２２及びステップＳ３０２ｂのｉに関するループの別の一部）を実施する。これらの計算の少なくとも一部は、並列して実施される。例えば、第１計算部２０Ａにおける第３関数の計算の上記の一部の実施の少なくとも一部と、第２計算部２０Ｂにおける第３関数の計算の別の一部の実施の少なくとも一部と、は、同時に行われる。並列計算により、計算が高速化できる。第３関数の計算量は多い。このため、第３関数（例えば、ステップＳ３０２ａ、ステップＳ１２２及びステップＳ３０２ｂのｉに関するループ）の計算を並列化することにより、効果的に高速化が行われる。

並列計算において、例えば、第１計算部２０Ａは、第３関数の計算の一部を実施するのに必要な、第１パラメータ群｛Ｊ｝の一部を、第１保持領域１０Ａに保持させる。このように、第１回路１５Ａのなかで、第３関数の計算の一部に必要な保持と処理が行われる。一方、第２計算部２０Ｂは、第３関数の計算の別の一部を実施するのに必要な、第１パラメータ群｛Ｊ｝の別の一部を、第２保持領域１０Ｂに保持させる。このように、第２回路１５Ｂのなかで、第３関数の計算の別の一部に必要な保持と処理が行われる。

例えば、第１パラメータ群｛Ｊ｝は、第１計算使用部分と、第２計算使用部分と、を含む。第１計算使用部分は、第３関数の計算の一部で用いられる。第２計算使用部分は、第３関数の計算の別の一部で用いられる。第１計算部２０Ａは、上記の第１計算使用部分を第１保持領域１０Ａに保持させる。第２計算部２０Ｂは、上記の第２計算使用部分を第２保持領域１０Ｂに保持させる。

実施形態において、ｉ番目の第１変数ｘ_ｉの第１変数更新及びｉ番目の第２変数ｙ_ｉの第１サブ更新（第２関数による更新）と、ｊ番目（「ｊ」は「ｉ」とは異なる）の第１変数ｘ_ｊの第１変数更新及びｊ番目の第２変数ｙ_ｊの第１サブ更新（第２関数による更新）と、を並列計算しても良い。この場合、処理部２０に、第３計算部２０Ｃ及び第４計算部２０Ｄが設けられても良い。第３計算部２０Ｃは、第１回路１５Ａに設けられる。第４計算部２０Ｄは、第２回路１５Ｂに設けられる。これらの計算部は、機能ブロックである。第３計算部２０Ｃの少なくとも一部で実施される処理は、第１計算部２０Ａの少なくとも一部で実施されても良い。第４計算部２０Ｄの少なくとも一部で実施される処理は、第２計算部２０Ｂの少なくとも一部で実施されても良い。

例えば、第３計算部２０Ｃは、第１変数更新の計算の一部、及び、第２関数の計算（第１サブ更新）の一部を実施する。第４計算部２０Ｄは、第１変数更新の計算の別の一部、及び、第２関数の計算（第１サブ更新）の別の一部を実施する。

既に説明したように、第１変数更新は、第１変数更新前のｉ番目の第１変数ｘ_ｉにｉ番目の第１関数を加えてｉ番目の第１変数ｘ_ｉを更新することを含む。同様に、第１変数更新は、第１変数更新前のｊ番目（ｊはｉとは異なり、１以上Ｎ以下の整数）の第１変数ｘ_ｊにｊ番目の第１関数を加えてｊ番目の第１変数ｘ_ｊを更新することをさらに含む。ここで、ｊ番目の第１変数ｘ_ｊは、第１変数群｛ｘ｝の１つである。ｊ番目の第１関数の変数は、ｊ番目の第２変数ｙ_ｊを含む。ｊ番目の第２変数ｙ_ｊは、第２変数群｛ｙ｝の１つである。

例えば、第１変数更新の計算の一部は、ｉ番目の第１変数ｘ_ｉの更新の計算を含む。第１変数更新の計算の別の一部は、ｊ番目の第１変数ｘ_ｊの更新の計算を含む。

一方、第２関数の計算の一部は、ｉ番目の第１変数ｘ_ｉを変数として含む第２関数の計算を含む。第２関数の計算の別の一部は、ｊ番目の第１変数ｘ_ｊを変数として含む第２関数の計算を含む。

第３計算部２０Ｃにおける、第１変数更新の計算の一部、及び、第２関数の計算（第１サブ更新）の一部の実施の少なくとも一部（例えば、図３〜１０のフローチャートにおけるステップＳ３０１ａとステップＳ３０１ｂとの間のｉループの一部）と、第４計算部２０Ｄにおける第１変数更新の計算の別の一部（例えば、図３〜１０のフローチャートにおけるステップＳ３０１ａとステップＳ３０１ｂとの間のｉループの別の一部）、及び、第２関数の計算（第１サブ更新）の別の一部の実施の少なくとも一部は、同時に行われても良い。このような並列計算により、高速に計算できる。

第２計算方法（第２計算装置）においては、上記の「第１変数更新の計算の一部」で更新される第１変数ｘに対応する第２変数ｙの２乗平均を、第１３式における２乗平均とする。これにより、例えば、「並列計算」において、２乗平均を計算するための、「第１変数更新の計算の一部」と「第１変数更新の計算の別の一部」との間での第２変数ｙの通信が、不要となる。並列計算による高速化が容易になる。部分集合「Ｓ_ｌ」が並列計算が可能になるように定義される。

例えば、処理部２０の１つの部分は、第２変数群{ｙ}の一部の更新と、第２変数群{ｙ}のその一部に関する２乗平均の算出と、を実施する。この場合、第２変数群{ｙ}の上記の一部は、上記の複数の部分集合の１つ（例えば、「Ｓ_ｌ」）に対応する。処理部２０の別の１つの部分は、第２変数群{ｙ}の別の一部の更新と、第２変数群{ｙ}のその別の一部に関する２乗平均の算出と、を実施する。この場合、第２変数群{ｙ}の上記の別の一部は、上記の複数の部分集合の別の１つ（例えば、「Ｓ_ｋ」）に対応する。例えば、処理部２０の上記の１つの部分で算出された２乗平均が、第２変数群{ｙ}の上記の一部の更新に用いられる。例えば、処理部２０の上記の別の１つの部分で算出された２乗平均が、第２変数群{ｙ}の上記の別の一部の更新に用いられる。

処理部２０の一部において、第２変数群{ｙ}の一部の更新と、その一部の２乗平均の算出と、が行われるため、通信（データの授受）が効率的になる。または、通信が不要になる。

例えば、第２変数群{ｙ}の上記の一部の更新と、第２変数群{ｙ}の上記の一部に関する２乗平均（または分散値）の算出と、の少なくとも一部の実施は、第２変数群{ｙ}の上記の別の一部の更新と、第２変数群{ｙ}の上記の別の一部に関する２乗平均(または分散値）の算出と、の少なくとも一部の実施と、同時に行われる。並列の計算により、処理時間が短縮できる。

例えば、処理部２０の１つの部分は、第２変数群{ｙ}の一部の更新（第１の更新）と、第２変数群{ｙ}の一部に関する２乗平均の算出（第１の算出）と、第３変数群｛ｚ｝に含まれる１つの第３変数ｚの更新（第２の更新）と、を実施する。第２の更新は、算出された２乗平均を用いて行われる。処理部２０の別の１つの部分は、第２変数群{ｙ}の別の一部の更新（第３の更新）と、第２変数群{ｙ}の別の一部に関する２乗平均の算出（第２の算出）と、第３変数群｛ｚ｝に含まれる別の１つの第３変数ｚの更新（第４の更新）と、を実施する。第４の更新は、算出された上記の別の一部に関する２乗平均を用いて行われる。

例えば、第１の更新、第１の算出、及び、第２の更新の少なくとも一部の実施は、第３の更新、第２の算出、及び、第４の更新の少なくとも一部の実施と同時に行われても良い。

以下、実施形態に係る計算装置の計算例について説明する。以下の計算例において、計算時間は、パラメータの設定の時間を含まない。計算時間は、パラメータの設定の後に、微分方程式を解くのに要した時間に対応する。

第１計算例では、ＰＣクラスタにより計算が行われる。第１計算例では、変数及びパラメータは、「float」（３２ビット浮動小数点実数）として扱われる。計算コアの数を「Ｑ」とする。「Ｑ」は、Ｎの約数である。Ｌ＝Ｎ／Ｑとする。

ＰＣクラスタで上記のアルゴリズムを並列して計算する際に、ＭＰＩ（Message Passing Interface）が用いられる。ＭＰＩは、分散メモリ型並列計算に対応する。ＭＰＩにおいては、複数の計算コアのそれぞれは、Ｌ個の第１分割変数（ｘ）と、Ｌ個の第２分割変数（ｙ）と、の１つの組み合わせを処理する。

例えば、ｉ番目の計算コアは、{ｘ_ｎ｜ｎ＝（ｉ−１）Ｌ＋１， …、 ｉＬ}、及び、{ｙ_ｎ｜ｎ＝（ｉ−１）Ｌ＋１， …、 ｉＬ}を保持し、これらの更新を行う。

ｉ番目の計算コアは、{ｈ_ｎ｜ｎ＝（ｉ−１）Ｌ＋１， …、 ｉＬ}、及び、{Ｊ_ｍ，ｎ｜ｍ＝（ｉ−１）Ｌ＋１， …、 ｉＬ；ｎ＝１…Ｎ}も保持できる。{ｈ_ｎ｜ｎ＝（ｉ−１）Ｌ＋１， …、 ｉＬ}、及び、{Ｊ_ｍ，ｎ｜ｍ＝（ｉ−１）Ｌ＋１， …、 ｉＬ；ｎ＝１…Ｎ}は、{ｙ_ｎ｜ｎ＝（ｉ−１）Ｌ＋１， …、 ｉＬ}の更新に用いられる。

例えば、{ｙ_ｎ｜ｎ＝（ｉ−１）Ｌ＋１， …、 ｉＬ}のそれぞれの更新には、{ｘ_ｎ｜ｎ＝１， …、 Ｎ}の全成分が用いられる。例えば、Allgather関数により、{ｘ_ｎ｜ｎ＝１， …、 Ｎ}の情報が、全ての計算コアに供給される。すなわち、情報（データ）が共有される。

実施形態において、複数の計算コアの間において、通信が行われる。すなわち、データの送受信が行われる。この通信において、{ｙ_ｎ｜ｎ＝（ｉ−１）Ｌ＋１， …、 ｉＬ}に関する通信、及び、{Ｊ_ｍ，ｎ｜ｍ＝（ｉ−１）Ｌ＋１， …、 ｉＬ；ｎ＝１…Ｎ}に関する通信は不要である。

第２計算方法（第２計算装置）においては、ｉ番目の計算コアは、{ｚ_ｎ｜ｎ＝（ｉ−１）Ｌ＋１， …、 ｉＬ}の更新に必要な第１３式の第２変数ｙの２乗平均を、{ｙ_ｎ｜ｎ＝（ｉ−１）Ｌ＋１， …、 ｉＬ}の２乗平均とする。これにより、第２変数ｙの２乗平均を計算するための、各コア間での第２変数ｙの通信が、不要となる。

例えば、第１変数群{ｘ}に関しての通信を行わず、第１パラメータ群{Ｊ}及び第１変数群{ｘ}の積和演算を分割して並列して実施し、その結果を通信して第２変数群{ｙ}の更新を行う方法が考えられる。この方法においては、第１パラメータ群{Ｊ}及び第１変数群{ｘ}の積和演算を分割して実施する。

以下、Ｎ＝２０００である場合（第１計算例）と、Ｎ＝１０００００である場合（第２計算例）と、の計算の例について説明する。

第１計算例（Ｎ＝２０００の場合）として、「Ｋ_２０００」問題（非特許文献１参照）の計算例について説明する。「Ｋ_２０００」問題は、Ｎ＝２０００の全結合イジングモデルである。行列Ｊの非対角成分は、±１のいずれかである。ベクトルｈの成分は、すべてゼロである。従って、ベクトルｈを含む項の計算は行われない。この場合において、行列Ｊの非対角成分の標準偏差σは、１である。このため、「ｃ」をｃ_０＝Ｄ／（２Ｎ^１／２）に設定する。「Ｋ_２０００」問題における行列Ｊの実際の最大固有値は、８８．８１３３２４である。一方、ランダム行列理論における理論値は、２σＮ^１／２＝８９．４４２７１９であり、これらの値は、互いに非常に近い。

以下に説明する第１計算例では、Ｑ＝２５であり、ｄｐ*（ｔ_１／ｄｔ）＝Ｄ＝２であり、ｔ_１＝５０である。

実施形態に係る第１計算方法（第１計算装置、第８式及び第１０式）おいて、Ｍ＝１で、ｄｔ＝０．２５の場合、計算時間は、７．６ｍｓである。このときに得られたイジングエネルギーの１００回の平均値は、約−６６０８６である。この値を「カット数」（第１５式、第１６式、及び、非特許文献１参照）に換算した値は、３２５２３である。「カット数」が大きいことは、精度が高いことに対応する。カット数Ｎｍｃは、以下の第１６式及び第１７式で表される。

実施形態に係る第１計算方法（第１計算装置、第８式及び第１０式）において、Ｍ＝５で、ｄｔ＝０．５の場合、計算時間は、４．１ｍｓである。このときに得られたイジングエネルギーの１００回の平均値は、約−６６１３７である。この値を「カット数」に換算すると、３２５４９である。Ｍ＝５の場合における「ｄｔ」は、Ｍ＝１の場合の「ｄｔ」の２倍にできる。Ｍ＝５の場合の計算時間は、Ｍ＝１の場合の計算時間の約半分になる。高速な計算が可能である。

一方、非特許文献１によると、コヒーレントイジングマシンにおいて、５ｍｓにおける１００回の「カット数」（非特許文献１参照）の平均値は、３２４５７である。一方、非特許文献１によると、シミュレーティッドアニーリングにおいて、５０ｍｓにおける１００回の「カット数」の平均値は、３２３１４である。このように、実施形態に係る計算により、コヒーレントイジングマシン及びシミュレーティッドアニーリングよりも短時間で、より高精度の解が得られる。

図１２は、実施形態に係る計算の特性を例示するグラフである。
図１２は、計算により得られる「カット数」と、「ｃ」と、の関係の例を示している。図１２の縦軸は、「ｃ／ｃ_０」である。既に説明したように、ｃ_０＝Ｄ／（２Ｎ^１／２）である。図１２の縦軸は、カット数Ｎｍｃである。

図１２に示すように、「ｃ／ｃ_０」が、約１以上、約１．５以下のときに大きなカット数Ｎｍｃが得られる。

以下、第２計算例（Ｎ＝１０００００の場合）について説明する。第２計算例では、行列Ｊの非対角成分、及び、ベクトルｈの成分は、「乱数」により設定される。「乱数」においては、−１〜１の値が一様に設定される。この場合の行列Ｊの非対角成分の標準偏差σは、１／（３^１／２）である。このため、ｃ＝３^１／２Ｄ／（２Ｎ^１／２）と設定する。その他のパラメータに関しては、Ｑ＝１２５０であり、ｄｐ*（ｔ_１／ｄｔ）＝Ｄ＝２であり、ｔ_１＝５０であり、ｄｔ＝０．５であり、Ｍ＝５である。

図１３は、計算結果を例示するグラフである。
図１３には、実施形態に係る第１計算方法（第１計算装置、第８式及び第１０式）での第２計算例の結果と、参考例に係る計算での第２計算例の結果と、が示されている。図１３の横軸は、計算時間Ｔｃ（秒）である。縦軸は、イジングエネルギーの１００回の平均値Ｖａｖｅである。図１３には、実施形態に係る計算結果１１０Ｅと、参考例の計算結果１１９Ｒと、が示されている。これらの計算例において、計算コアの数は、１２５０である。

実施形態に係る第１計算方法（第１計算装置、第８式及び第１０式）の計算結果１１０Ｅにおいて、「ｎ」は、第１４式の非線形関数を用いた場合の第１４式中の「ｎ」の値を示している。計算結果１１０Ｅにおいて、「ｎ」が書かれていない曲線は、第８式が用いられた場合に対応する。

参考例の計算結果１１９Ｒでは、シミュレーティッドアニーリングによる計算が行われる。シミュレーティッドアニーリングでは、スピン反転によるエネルギー変化が、ＭＰＩにより、並列計算される。これらの計算においては、逆温度を線形に増加させている。参考例の計算結果１１９Ｒの複数の曲線において、増加率が互いに異なる。

図１３から分かるように、実施形態に係る第１計算方法（第１計算装置、第８式及び第１０式）の計算結果１１０Ｅにおいて、最終的な平均値Ｖａｖｅは、低い（絶対値が大きい）。これに対して、参考例の計算結果においては、最終的な平均値Ｖａｖｅは十分に低くならない（絶対値が十分に大きくならない）。このように、実施形態においては、高い精度の計算結果が得られる。実施形態によれば、参考例（シミュレーティッドアニーリング）に比べて、同じ精度が得られるのに要する計算時間Ｔｃは、１／１０以下である。実施形態に係る計算は、参考例に比べて、１０倍以上高速である。

以下、本実施形態に係る計算の第３計算例について説明する。第３計算例では、ＧＰＵにより、上記の計算が行われる。この計算において、例えば、変数及びパラメータは、float（３２ビット浮動小数点実数）として扱われる。

この方法では、第１変数群{ｘ}、第２変数群{ｙ}、第１パラメータ群{Ｊ}及び第２パラメータ群{ｈ}は、デバイス変数として定義される。第１パラメータ群{Ｊ}は、行列Ｊである。行列Ｊ及び第１変数ｘの積和演算を用いた第２変数ｙの更新は、行列ベクトル積関数により行われる。第１変数ｘ及び第２変数ｙに関するその他の更新に関して、ｉ番目の成分（ｘ_ｉ，ｙ_ｉ）の更新が１つのスレッドで行われる。

第３計算例では、１つのＧＰＵを用いて、第１計算例と同様の条件で、「Ｋ_２０００」問題が計算される。第３計算例の計算時間は、１４．７ｍｓ、１００回の「カット数」の平均値３２５４９が得られる。第３計算例においても、参考例のシミュレーティッドアニーリングの結果よりも高速に計算できる。

以下、第１計算方法及び第２計算方法（第１計算装置及び第２計算装置）の例について説明する。

図１４（ａ）及び図１４（ｂ）は、計算結果を例示するグラフである。
これらの図は、「Ｇ２２」問題（非特許文献１参照）に関しての計算結果の例を示す。これらの図には、第１計算方法（第１計算装置、第８式及び第１０式）による計算結果Ｍ１１と、第２計算方法（第２計算装置、第１１式〜第１３式）による計算結果Ｍ１２と、が示されている。計算結果Ｍ１１において、第２変数ｙの初期値ｙｉ（０）は、−０．１以上０．１以下の範囲の乱数の値である。計算結果Ｍ１２において、第２変数ｙの初期値ｙｉ（０）は、−０．１以上０．１以下の範囲の乱数の値である。第２計算方法（第２計算装置、第１１式〜第１３式）において、第１計算パラメータＴ（ｔ）は、処理手順（ループ）ごとに、１から０に向けて、単調に減少する。乱数での値を変えて、１００回の計算が行われる。

図１４（ａ）は、１００回の計算における平均値である。図１４（ｂ）は、１００回の計算における最大値である。これらの図の横軸は、処理手順の繰り返し回数ＮＬである。縦軸は、カット数Ｎｍｃである。１つのループにおいて、Ｎ個の第１変数ｘとＮ個の第２変数ｙとが全て更新される。１つのループにおいて、第３変数ｚも更新される。

図１４（ａ）及び図１４（ｂ）に示すように、処理手順の繰り返し回数ＮＬが３０００以上において、第２計算方法（第２計算装置、第１１式〜第１３式）によるカット数Ｎｍｃは、第１計算方法（第１計算装置、第８式及び第１０式）によるカット数Ｎｍｃよりも大きい。例えば、比較的時間をかけて良い解を探索したい場合、第２計算方法（第２計算装置）が有利であると考えられる。

第２計算方法（第２計算装置）におけるｉ番目の第１関数に含まれる第１計算パラメータＴ（ｔ）は、例えば、温度に対応する。第２計算方法（第２計算装置）によれば、例えば、大きな初期運動量（第２変数ｙの初期値ｙｉ（０））でも、良い解に収束する。大きな初期運動量により、より広い範囲の探索が容易になる。

(第２実施形態）
第２実施形態に係る計算装置（計算方法）においては、上記の第２計算方法に関して説明した第１１式〜第１３式の代わりに、以下の第１７式〜第２０式が用いられる。

第１７式に示すように、第３変数ｂが導入される。第３変数ｂは、第１実施形態に関して説明した第３変数群{ｚ}とは異なる。

第２実施形態に係る計算装置は、以下に説明する処理手順を繰り返す処理部２０を含む。処理手順は、第１変数更新、第２変数更新、第３変数更新及び第４変数更新を含む。

第１変数更新は、例えば、第１７式に基づいて行われる。第１変数更新は、例えば、第１変数更新前のｉ番目の第１変数ｘ_ｉにｉ番目の第１関数を加えてｉ番目の第１変数ｘ_ｉを更新することを含む。ｉ番目の第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つである。ｉ番目の第１関数の変数は、ｉ番目の第２変数ｙ_ｉ及び第３変数ｂを含む。ｉ番目の第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つである。ｉ番目の第１関数は、第１関数群の１つである。第１７式に示す例では、第１関数は、ｙ_ｉ／ｂ^２である。

第２変数更新は、第１８式に基づいて行われる。第２変数更新は、例えば、第２変数更新前のｉ番目の第２変数ｙ_ｉに、ｉ番目の第２関数及びｉ番目の第３関数を加えてｉ番目の第２変数ｙ_ｉを更新することを含む。ｉ番目の第２関数の変数は、ｉ番目の第１変数ｘ_ｉを含む。ｉ番目の第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部、及び、第１変数群｛ｘ｝の少なくとも一部を含む。

ｉ番目の第２関数は、例えば、第１８式の「ｘ_ｉ］を含む項に対応する。ｉ番目の第３関数は、第１８式の「ｊ」に関する和の項に対応する。

第３変数更新（第３変数ｂの更新）は、第１９式に基づいて行われる。第３変数更新は、例えば、第３変数更新前の第３変数ｂに第１要素関数を加えて、第３変数ｂを更新することを含む。第３変数ｂは、第３変数群｛ｂ｝の１つである。第１要素関数の変数は、第４変数Ｐ_ｂを含む。第１９式の例では、第１要素関数は、「（Ｐ_ｂ／Ｖ）」に対応する。例えば、「Ｖ」は、自由度（第３変数ｂ）の質量に相当する。

第４変数更新（第４変数Ｐ_ｂの更新）は、例えば、第２０式に基づいて行われる。第４変数更新は、例えば、第４変数更新前の第４変数Ｐ_ｂに第２要素関数を加えて第４変数Ｐ_ｂを更新することを含む。第４変数Ｐ_ｂは、第４変数群｛Ｐ_ｂ｝の１つである。

第２０式に示す例では、第２要素関数の変数は、第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含む。２乗平均は、第２０式の第１項に対応する。第１計算パラメータは、第２０式の第２項に対応する。

処理部２０は、少なくとも、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉ、及び、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉの関数の少なくともいずれかを出力する。

第２実施形態においても、同時更新（並列更新）が可能である。例えば、乱数が不要である。例えば、処理手順の繰り返しは、積和演算の実施により実施できる。最適化問題を高速に計算できる計算装置が提供できる。

第２実施形態に係る計算装置（及び計算方法）に、第１実施形態に関して説明した処理（例えば、並列処理など）が適用できる。

(第３実施形態）
第３実施形態に係る計算装置（計算方法）においては、上記の第２計算方法に関して説明した第１式から第１３式の代わりに、以下の第２１式〜第２６式が用いられる。

第２１式に示すように、第３変数ｂが導入される。第３変数ｂは、第１実施形態に関して説明した第３変数群{ｚ}とは異なる。

第３実施形態に係る計算装置は、以下に説明する処理手順を繰り返す処理部２０を含む。処理手順は、第１変数更新、第２変数更新、第３変数更新及び第４変数更新を含む。

第１変数更新は、例えば、第２１式に基づいて行われる。第１変数更新は、第１変数更新前のｉ番目の第１変数ｘ_ｉにｉ番目の第１関数を加えて、ｉ番目の第１変数ｘ_ｉを更新することを含む。ｉ番目の第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つである。ｉ番目の第１関数の変数は、ｉ番目の第２変数ｙ_ｉ及び第３変数ｂを含む。ｉ番目の第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つである。ｉ番目の第１関数は、第１関数群の１つである。

第２１式の例においては、ｉ番目の第１関数は、ｙ_ｉ／ｂである。

第２変数更新は、第２２式に基づいて行われる。第２変数更新は、第２変数更新前のｉ番目の第２変数ｙ_ｉに、ｉ番目の第２関数及びｉ番目の第３関数の和と、第３変数ｂと、の積を加えて、ｉ番目の第２変数ｙ_ｉを更新することを含む。ｉ番目の第２関数の変数は、ｉ番目の第１変数ｘ_ｉを含む。ｉ番目の第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び第１変数群｛ｘ｝の少なくとも一部を含む。

第２２式の例においては、ｉ番目の第２関数は、例えば、第２２式の「ｘ_ｉ］を含む項に対応する。ｉ番目の第３関数は、第２２式の「ｊ」に関する和の項に対応する。

第３変数更新（第３変数ｂの更新）は、例えば、第２３式に基づいて行われる。第３変数更新は、第３変数更新前の第３変数ｂに第１要素関数を加えて、第３変数ｂを更新することを含む。第３変数ｂは、第３変数群｛ｂ｝の１つである。第１要素関数の変数は、第４変数Ｐ_ｂを含む。

第２３式に示すように、第１要素関数の要素は、第３変数ｂと、第４変数Ｐ_ｂと、の積を含む。第２３式の例では、第１要素関数は、「（ｂ・Ｐ_ｂ）／Ｖ」に対応する。

第４変数更新（第４変数Ｐ_ｂの更新）は、例えば、第２４式に基づいて行われる。第４変数更新は、第４変数更新前の第４変数Ｐ_ｂに第２要素関数を加えて、第４変数Ｐ_ｂを更新することを含む。第４変数Ｐ_ｂは、第４変数群｛Ｐ_ｂ｝の１つである。第２要素関数の変数は、第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含む。

第２４式に示す例では、第２要素関数の変数は、第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含む。この２乗平均は、第２４式の第１項に対応する。第１計算パラメータは、第２４式の第２項及び第３項に対応する。

処理部２０は、少なくとも、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉ、及び、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉの関数の少なくともいずれかを出力する。

第２５式において、「Ｂ」は、第２６式が満たされるように決められる定数である。

第３実施形態においても、同時更新（並列更新）が可能である。例えば、乱数が不要である。例えば、処理手順の繰り返しは、積和演算の実施により実施できる。最適化問題を高速に計算できる計算装置が提供できる。

第３実施形態に係る計算装置（及び計算方法）に、第１実施形態に関して説明した処理（例えば、並列処理など）が適用できる。

（第４実施形態）
第４実施形態は、第１〜第３実施形態に関して説明した計算が可能な回路を含む。

（第５実施形態）
第５実施形態は、計算プログラムに係る。計算プログラムは、コンピュータに、第１〜第３実施形態に関して説明した処理手順を繰り返させる。以下、第１実施形態に関して説明した処理手順を実施する場合について説明する。

この処理手順は、第１変数更新、第２変数更新及び第３変数更新を含む。
第１変数更新は、第１変数更新前のｉ番目（ｉは１以上Ｎ以下の整数であり、Ｎは２以上の整数）の第１変数ｘ_ｉにｉ番目の第１関数を加えてｉ番目の第１変数ｘ_ｉを更新することを含む。ｉ番目の第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つである。ｉ番目の第１関数の変数は、ｉ番目の第２変数ｙ_ｉ含み、ｉ番目の第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つである。ｉ番目の第１関数は、第１関数群の１つである。

第２変数更新は、第２変数更新前のｉ番目の第２変数ｙ_ｉに、ｉ番目の第２関数、ｉ番目の第３関数及びｉ番目の第１要素関数の演算結果を加えてｉ番目の第２変数ｙ_ｉを更新することを含む。ｉ番目の第２関数の変数は、ｉ番目の第１変数ｘ_ｉを含む。ｉ番目の第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び第１変数群｛ｘ｝の少なくとも一部を含む。ｉ番目の第１要素関数の変数は、ｉ番目の第２変数ｙ_ｉと第３変数ｚとの積を含む。ｉ番目の第２関数は、第２関数群の１つである。ｉ番目の第３関数は、第３関数群の１つである。ｉ番目の第１要素関数は、第１要素関数群の１つである。

第３変数更新は、第３変数更新前の第３変数ｚにｉ番目の第２要素関数を加えて第３変数ｚを更新することを含む。第３変数ｚは、第３変数群｛ｚ｝の１つである。第２要素関数の変数は、第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含む。

処理部は、少なくとも、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉ、及び、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉの関数の少なくともいずれかを出力する。

本実施形態に係る計算プログラムは、第２実施形態及び第３実施形態に関して説明した処理が適用できる。

（第６実施形態）
第６実施形態は、コンピュータ読み取り可能な記録媒体である。この記録媒体は、コンピュータに、第１〜第３実施形態に関して説明した処理手順を繰り返させるプログラムが記録されている。以下、第１実施形態に関して説明した処理手順を実施する場合について説明する。

この処理手順は、第１変数更新、第２変数更新及び第３変数更新を含む。
第１変数更新は、第１変数更新前のｉ番目（ｉは１以上Ｎ以下の整数であり、Ｎは２以上の整数）の第１変数ｘ_ｉにｉ番目の第１関数を加えてｉ番目の第１変数ｘ_ｉを更新することを含む。ｉ番目の第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つである。ｉ番目の第１関数の変数は、ｉ番目の第２変数ｙ_ｉ含み、ｉ番目の第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つである。ｉ番目の第１関数は、第１関数群の１つである。

第２変数更新は、第２変数更新前のｉ番目の第２変数ｙ_ｉに、ｉ番目の第２関数、ｉ番目の第３関数及びｉ番目の第１要素関数の演算結果を加えてｉ番目の第２変数ｙ_ｉを更新することを含む。ｉ番目の第２関数の変数は、ｉ番目の第１変数ｘ_ｉを含む。ｉ番目の第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び第１変数群｛ｘ｝の少なくとも一部を含む。ｉ番目の第１要素関数の変数は、ｉ番目の第２変数ｙ_ｉと第３変数ｚとの積を含む。ｉ番目の第２関数は、第２関数群の１つである。ｉ番目の第３関数は、第３関数群の１つである。ｉ番目の第１要素関数は、第１要素関数群の１つである。

第３変数更新は、第３変数更新前の第３変数ｚにｉ番目の第２要素関数を加えて第３変数ｚを更新することを含む。第３変数ｚは、第３変数群｛ｚ｝の１つである。第２要素関数の変数は、第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含む。

処理部は、少なくとも、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉ、及び、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉの関数の少なくともいずれかを出力する。

本実施形態に係る記録媒体には、第２実施形態及び第３実施形態に関して説明した処理が適用できる。

（第７実施形態）
本実施形態は、計算方法に係る。計算方法は、第１〜第３実施形態に関して説明した処理手順を行う。以下、第１実施形態に関して説明した処理手順を実施する場合について説明する。

この処理手順は、第１変数更新、第２変数更新及び第３変数更新を含む。
第１変数更新は、第１変数更新前のｉ番目（ｉは１以上Ｎ以下の整数であり、Ｎは２以上の整数）の第１変数ｘ_ｉにｉ番目の第１関数を加えてｉ番目の第１変数ｘ_ｉを更新することを含む。ｉ番目の第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つである。ｉ番目の第１関数の変数は、ｉ番目の第２変数ｙ_ｉ含み、ｉ番目の第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つである。ｉ番目の第１関数は、第１関数群の１つである。

第２変数更新は、第２変数更新前のｉ番目の第２変数ｙ_ｉに、ｉ番目の第２関数、ｉ番目の第３関数及びｉ番目の第１要素関数の演算結果を加えてｉ番目の第２変数ｙ_ｉを更新することを含む。ｉ番目の第２関数の変数は、ｉ番目の第１変数ｘ_ｉを含む。ｉ番目の第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び第１変数群｛ｘ｝の少なくとも一部を含む。ｉ番目の第１要素関数の変数は、ｉ番目の第２変数ｙ_ｉと第３変数ｚとの積を含む。ｉ番目の第２関数は、第２関数群の１つである。ｉ番目の第３関数は、第３関数群の１つである。ｉ番目の第１要素関数は、第１要素関数群の１つである。

第３変数更新は、第３変数更新前の第３変数ｚにｉ番目の第２要素関数を加えて第３変数ｚを更新することを含む。第３変数ｚは、第３変数群｛ｚ｝の１つである。第２要素関数の変数は、第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含む。

処理部は、少なくとも、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉ、及び、処理手順を繰り返した後に得られるｉ番目の第１変数ｘ_ｉの関数の少なくともいずれかを出力する。

本実施形態に係る記録媒体には、第２実施形態及び第３実施形態に関して説明した処理が適用できる。

上記の種々の情報（データ）の処理（指示）は、例えば、プログラム（ソフトウェア）に基づいて実行される。例えば、コンピュータが、このプログラムを記憶し、このプログラムを読み出すことにより、上記の種々の情報の処理が行われる。

上記の種々の情報の処理は、コンピュータに実行させることのできるプログラムとして、磁気ディスク（フレキシブルディスク及びハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＣＤ−Ｒ、ＣＤ−ＲＷ、ＤＶＤ−ＲＯＭ、ＤＶＤ±Ｒ、ＤＶＤ±ＲＷなど）、半導体メモリ、または、他の記録媒体に記録されても良い。

例えば、記録媒体に記録された情報は、コンピュータ（または組み込みシステム）により読み出されることが可能である。記録媒体において、記録形式（記憶形式）は任意である。例えば、コンピュータは、記録媒体からプログラムを読み出し、このプログラムに基づいてプログラムに記述されている指示をＣＰＵで実行させる。コンピュータにおいて、プログラムの取得（または読み出し）は、ネットワークを通じて行われても良い。

記録媒体からコンピュータ（または組み込みシステム）にインストールされたプログラムに基づいてコンピュータ上で稼働している種々のソフトウェアにおいて、上記の情報の処理の少なくとも一部が実施されても良い。このソフトウェアは、例えば、ＯＳ（オペレーティングシステム）などを含む。このソフトウェアは、例えば、ネットワーク上で動作するミドルウェアなどを含んでも良い。

実施形態における記録媒体は、ＬＡＮまたはインターネットなどにより得られたプログラムをダウンロードして記憶された記録媒体も含まれる。複数の記録媒体に基づいて、上記の処理が行われても良い。

実施形態に係るコンピュータは、１つまたは複数の装置（例えばパーソナルコンピュータなど）を含む。実施形態に係るコンピュータは、ネットワークにより接続された複数の装置を含んでも良い。

実施形態は、例えば、以下の構成（例えば技術案）を含んでも良い。
（構成１）
処理手順を繰り返す処理部を備え、
前記処理手順は、第１変数更新、第２変数更新及び第３変数更新を含み、
前記第１変数更新は、前記第１変数更新前の前記ｉ番目（ｉは１以上Ｎ以下の整数であり、前記Ｎは２以上の整数）の第１変数ｘ_ｉにｉ番目の第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することを含み、前記ｉ番目の前記第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つであり、前記ｉ番目の前記第１関数の変数は、前記ｉ番目の第２変数ｙ_ｉ含み、前記ｉ番目の前記第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つであり、前記ｉ番目の前記第１関数は、第１関数群の１つであり、
前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉに、前記ｉ番目の第２関数、前記ｉ番目の第３関数及び前記ｉ番目の第１要素関数の演算結果を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、前記ｉ番目の前記第２関数の変数は、前記ｉ番目の前記第１変数ｘ_ｉを含み、前記ｉ番目の前記第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び前記第１変数群｛ｘ｝の少なくとも一部を含み、前記ｉ番目の前記第１要素関数の変数は、前記ｉ番目の前記第２変数ｙ_ｉと第３変数ｚとの積を含み、前記ｉ番目の前記第２関数は、第２関数群の１つであり、前記ｉ番目の前記第３関数は、第３関数群の１つであり、前記ｉ番目の前記第１要素関数は、第１要素関数群の１つであり、
前記第３変数更新は、前記第３変数更新前の第３変数ｚに前記ｉ番目の第２要素関数を加えて前記第３変数ｚを更新することを含み、前記第３変数ｚは、第３変数群｛ｚ｝の１つであり、前記第２要素関数の変数は、前記第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含み、
前記処理部は、少なくとも、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉ、及び、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉの関数の少なくともいずれかを出力する、計算装置。

（構成２）
前記処理手順の後の前記第１計算パラメータは、前記処理手順の前の前記第１計算パラメータとは異なる、構成１記載の計算装置。

（構成３）
前記処理手順の後の前記第１計算パラメータは、前記処理手順の前の前記第１計算パラメータよりも小さい、構成１記載の計算装置。

（構成４）
前記第２変数更新は、第１サブ更新及び第２サブ更新を含み、
前記第１サブ更新は、前記第１サブ更新前の前記ｉ番目の前記第２変数ｙ_ｉに前記第２関数を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、
前記第２サブ更新は、前記第２サブ更新前の前記ｉ番目の前記第２変数ｙ_ｉに前記第３関数を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、
前記第１変数更新、前記第１サブ更新及び前記第３変数更新をＭ回（Ｍは２以上の整数）行った後、前記第２サブ更新を行う、または、前記第２サブ更新後、前記第１変数更新、前記第１サブ更新及び前記第３変数更新をＭ回行う、構成１〜３のいずれか１つに記載の計算装置。

（構成５）
前記第３変数更新の少なくとも一部は、前記第２サブ更新の少なくとも一部と同時に実施される、構成４記載の計算装置。

（構成６）
前記処理部の１つの部分は、前記第２変数群{ｙ}の前記一部の更新と、前記第２変数群{ｙ}の前記一部に関する前記２乗平均の算出と、を実施し、
前記処理部の別の１つの部分は、前記第２変数群{ｙ}の別の一部の更新と、第２変数群{ｙ}の別の一部に関する２乗平均の算出と、を実施し、
前記第２変数群{ｙ}の前記一部の前記更新と、前記第２変数群{ｙ}の前記一部に関する前記２乗平均の前記算出と、の少なくとも一部の実施は、前記第２変数群{ｙ}の前記別の一部の前記更新と、前記第２変数群{ｙ}の前記別の一部に関する前記２乗平均の前記算出と、の少なくとも一部の実施と同時に行われる、構成１〜５のいずれか１つに記載の計算装置。

（構成７）
前記処理部の１つの部分は、第１の更新、第１の算出、及び、第２の更新を実施し、前記第１の更新は、前記第２変数群{ｙ}の前記一部の更新を含み、前記第１の算出は、前記第２変数群{ｙ}の一部に関する２乗平均の算出を含み、前記第２の更新は、前記第３変数群｛ｚ｝に含まれる１つの前記第３変数ｚの前記２乗平均を用いた更新を含み、
前記処理部の別の１つの部分は、第３の更新、第２の算出、及び、第４の更新を実施し、前記第３の更新は、前記第２変数群{ｙ}の別の一部の更新を含み、前記第２の算出は、第２変数群{ｙ}の前記別の一部に関する前記２乗平均の算出を含み、前記第４の更新は、前記第３変数群｛ｚ｝に含まれる別の１つの前記第３変数ｚの前記２乗平均を用いた更新を含み、
前記第１の更新、前記第１の算出、及び、前記第２の更新の少なくとも一部の実施は、前記第３の更新、前記第２の算出、及び、前記第４の更新の少なくとも一部の実施と同時に行われる、構成１〜５のいずれか１つに記載の計算装置。

（構成８）
前記ｉ番目の前記第２関数は、前記ｉ番目の前記第１変数ｘ_ｉを変数とする前記ｉ番目の第４関数を含み、
前記ｉ番目の第４関数は、前記第１変数ｘ_ｉの非線形関数を含み、
前記ｉ番目の第４関数は、演算パラメータｐを含み、
前記演算パラメータｐは、前記処理手順を繰り返すと変化し、
前記処理手順を繰り返した後の前記ｉ番目の前記第４関数の実数根の数は、前記処理手順を繰り返す前の前記ｉ番目の前記第４関数の実数根の数とは異なる、構成１〜７のいずれか１つに記載の計算装置。

（構成９）
前記第１変数更新は、前記第１変数更新前の前記ｉ番目の前記第１変数ｘ_ｉを保持部から取得し、前記第１変数更新後の前記ｉ番目の前記第１変数ｘ_ｉを前記保持部に保持させることを含み、
前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉを前記保持部から取得し、前記第２変数更新後の前記ｉ番目の前記第２変数ｙ_ｉを前記保持部に保持させることを含む、構成１〜７のいずれか１つに記載の計算装置。

（構成１０）
前記第３変数更新は、前記第３変数更新前の前記第３変数ｚを前記保持部から取得し、前記第３変数更新後の前記第３変数ｚを前記保持部に保持させることを含む、構成９に記載の計算装置。

（構成１１）
前記第１変数更新は、前記第２変数群{ｙ}の少なくとも一部を前記保持部から取得して前記ｉ番目の前記第１関数を計算し、前記ｉ番目の前記第１変数ｘ_ｉに前記ｉ番目の前記第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することをさらに含み、
前記第２変数更新は、前記ｉ番目の前記第１変数ｘ_ｉを前記保持部から取得して前記ｉ番目の前記第２関数を計算し、前記第１パラメータ群｛Ｊ｝の前記少なくとも一部及び前記第１変数群｛ｘ｝の前記少なくとも一部を前記保持部から取得して前記ｉ番目の前記第３関数を計算し、前記ｉ番目の前記第２変数ｙ_ｉに前記ｉ番目の前記第２関数及び前記ｉ番目の前記第３関数を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することをさらに含む、を含む、構成９または１０記載の計算装置。

（構成１２）
前記処理部は、第１計算部と、第２計算部と、を含み、
前記第１計算部は、前記第３関数の計算の一部を実施し、
前記第２計算部は、前記第３関数の前記計算の別の一部を実施し、
前記第１計算部における前記第３関数の前記計算の前記一部の前記実施の少なくとも一部と、前記第２計算部における前記第３関数の前記計算の前記別の一部の前記実施の少なくとも一部は、同時に行われる、構成９〜１１のいずれか１つに記載の計算装置。

（構成１３）
前記保持部は、第１保持領域及び第２保持領域を含み、
前記第１パラメータ群｛Ｊ｝は、第１計算使用部分と、第２計算使用部分と、を含み、前記第１計算使用部分は、前記第３関数の前記計算の前記一部で用いられ、前記第２計算使用部分は、前記第３関数の前記計算の前記別の一部で用いられ、
前記第１計算部は、前記第１計算使用部分を前記第１保持領域に保持させ、
前記第２計算部は、前記第２計算使用部分を前記第２保持領域に保持させる、構成９〜１２のいずれか１つに記載の計算装置。

（構成１４）
処理手順を繰り返す処理部を備え、
前記処理手順は、第１変数更新、第２変数更新、第３変数更新及び第４変数更新を含み、
前記第１変数更新は、前記第１変数更新前のｉ番目（ｉは１以上Ｎ以下の整数であり、前記Ｎは２以上の整数）の第１変数ｘ_ｉに前記ｉ番目の第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することを含み、前記ｉ番目の前記第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つであり、前記ｉ番目の前記第１関数の変数は、前記ｉ番目の第２変数ｙ_ｉ及び第３変数ｂを含み、前記ｉ番目の前記第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つであり、前記ｉ番目の前記第１関数は、第１関数群の１つであり、
前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉに前記ｉ番目の第２関数及び前記ｉ番目の第３関数を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、前記ｉ番目の前記第２関数の変数は、前記ｉ番目の前記第１変数ｘ_ｉを含み、前記ｉ番目の前記第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び前記第１変数群｛ｘ｝の少なくとも一部を含み、
前記第３変数更新は、前記第３変数更新前の第３変数ｂに第１要素関数を加えて前記第３変数ｂを更新することを含み、前記第１要素関数の変数は、第４変数Ｐ_ｂを含み、
前記第４変数更新（Ｐ_ｂの更新）は、前記第４変数更新前の第４変数Ｐ_ｂに第２要素関数を加えて前記第４変数Ｐ_ｂを更新することを含み、前記第２要素関数の変数は、前記第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含み、
前記処理部は、少なくとも、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉ、及び、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉの関数の少なくともいずれかを出力する、計算装置。

（構成１５）
前記ｉ番目の前記第１関数は、ｙ_ｉ／ｂ^２である、構成１４記載の計算装置。

（構成１６）
処理手順を繰り返す処理部を備え、
前記処理手順は、第１変数更新、第２変数更新、第３変数更新及び第４変数更新を含み、
前記第１変数更新は、前記第１変数更新前のｉ番目（ｉは１以上Ｎ以下の整数であり、前記Ｎは２以上の整数）の第１変数ｘ_ｉに前記ｉ番目の第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することを含み、前記ｉ番目の前記第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つであり、前記ｉ番目の前記第１関数の変数は、前記ｉ番目の第２変数ｙ_ｉ及び第３変数ｂを含み、前記ｉ番目の前記第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つであり、前記ｉ番目の前記第１関数は、第１関数群の１つであり、
前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉに、前記ｉ番目の第２関数及び前記ｉ番目の第３関数の和と、前記第３変数ｂと、の積を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、前記ｉ番目の前記第２関数の変数は、前記ｉ番目の前記第１変数ｘ_ｉを含み、前記ｉ番目の前記第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び前記第１変数群｛ｘ｝の少なくとも一部を含み、
前記第３変数更新は、前記第３変数更新前の第３変数ｂに第１要素関数を加えて前記第３変数ｂを更新することを含み、前記第１要素関数の変数は、第４変数Ｐ_ｂを含み、
前記第４変数更新は、前記第４変数更新前の第４変数Ｐ_ｂに第２要素関数を加えて前記第４変数Ｐ_ｂを更新することを含み、前記第２要素関数の変数は、前記第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含み、
前記処理部は、少なくとも、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉ、及び、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉの関数の少なくともいずれかを出力する、計算装置。

（構成１７）
前記第１要素関数の前記要素は、前記第３変数ｂと、前記第４変数Ｐ_ｂと、の積を含む構成１６記載の計算装置。

（構成１８）
コンピュータに、処理手順を繰り返させる計算プログラムであって、
前記処理手順は、第１変数更新、第２変数更新及び第３変数更新を含み、
前記第１変数更新は、前記第１変数更新前の前記ｉ番目（ｉは１以上Ｎ以下の整数であり、前記Ｎは２以上の整数）の第１変数ｘ_ｉに前記ｉ番目の第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することを含み、前記ｉ番目の前記第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つであり、前記ｉ番目の前記第１関数の変数は、前記ｉ番目の第２変数ｙ_ｉ含み、前記ｉ番目の前記第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つであり、前記ｉ番目の前記第１関数は、第１関数群の１つであり、
前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉに、前記ｉ番目の第２関数、前記ｉ番目の第３関数及び前記ｉ番目の第１要素関数の演算結果を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、前記ｉ番目の前記第２関数の変数は、前記ｉ番目の前記第１変数ｘ_ｉを含み、前記ｉ番目の前記第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び前記第１変数群｛ｘ｝の少なくとも一部を含み、前記ｉ番目の前記第１要素関数の変数は、前記ｉ番目の前記第２変数ｙ_ｉと第３変数ｚとの積を含み、前記ｉ番目の前記第２関数は、第２関数群の１つであり、前記ｉ番目の前記第３関数は、第３関数群の１つであり、前記ｉ番目の前記第１要素関数は、第１要素関数群の１つであり、
前記第３変数更新は、前記第３変数更新前の第３変数ｚに前記ｉ番目の第２要素関数を加えて前記ｉ番目の前記第３変数ｚを更新することを含み、前記ｉ番目の前記第３変数ｚは、第３変数群｛ｚ｝の１つであり、前記第２要素関数の変数は、前記第２変数群｛ｙ｝の少なくとも一部の２乗平均値と、第１計算パラメータと、の差を含み、
前記処理部は、少なくとも、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉ、及び、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉの関数の少なくともいずれかを出力する、計算プログラム。

（構成１９）
コンピュータに、処理手順を繰り返させる計算プログラムを記録したコンピュータ読み取り可能な記録媒体であって、
前記処理手順は、第１変数更新、第２変数更新及び第３変数更新を含み、
前記第１変数更新は、前記第１変数更新前の前記ｉ番目（ｉは１以上Ｎ以下の整数であり、前記Ｎは２以上の整数）の第１変数ｘ_ｉに前記ｉ番目の第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することを含み、前記ｉ番目の前記第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つであり、前記ｉ番目の前記第１関数の変数は、前記ｉ番目の第２変数ｙ_ｉ含み、前記ｉ番目の前記第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つであり、前記ｉ番目の前記第１関数は、第１関数群の１つであり、
前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉに、前記ｉ番目の第２関数、前記ｉ番目の第３関数及び前記ｉ番目の第１要素関数の演算結果を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、前記ｉ番目の前記第２関数の変数は、前記ｉ番目の前記第１変数ｘ_ｉを含み、前記ｉ番目の前記第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び前記第１変数群｛ｘ｝の少なくとも一部を含み、前記ｉ番目の前記第１要素関数の変数は、前記ｉ番目の前記第２変数ｙ_ｉと第３変数ｚとの積を含み、前記ｉ番目の前記第２関数は、第２関数群の１つであり、前記ｉ番目の前記第３関数は、第３関数群の１つであり、前記ｉ番目の前記第１要素関数は、第１要素関数群の１つであり、
前記第３変数更新は、前記第３変数更新前の第３変数ｚに前記ｉ番目の第２要素関数を加えて前記ｉ番目の前記第３変数ｚを更新することを含み、前記ｉ番目の前記第３変数ｚは、第３変数群｛ｚ｝の１つであり、前記第２要素関数の変数は、前記第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含み、
前記処理部は、少なくとも、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉ、及び、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉの関数の少なくともいずれかを出力する、記録媒体。

（構成２０）
処理手順を繰り返し、
前記処理手順は、第１変数更新、第２変数更新及び第３変数更新を含み、
前記第１変数更新は、前記第１変数更新前の前記ｉ番目（ｉは１以上Ｎ以下の整数であり、前記Ｎは２以上の整数）の第１変数ｘ_ｉに前記ｉ番目の第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することを含み、前記ｉ番目の前記第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つであり、前記ｉ番目の前記第１関数の変数は、前記ｉ番目の第２変数ｙ_ｉ含み、前記ｉ番目の前記第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つであり、前記ｉ番目の前記第１関数は、第１関数群の１つであり、
前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉに、前記ｉ番目の第２関数、前記ｉ番目の第３関数及び前記ｉ番目の第１要素関数の演算結果を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、前記ｉ番目の前記第２関数の変数は、前記ｉ番目の前記第１変数ｘ_ｉを含み、前記ｉ番目の前記第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び前記第１変数群｛ｘ｝の少なくとも一部を含み、前記ｉ番目の前記第１要素関数の変数は、前記ｉ番目の前記第２変数ｙ_ｉと第３変数ｚとの積を含み、前記ｉ番目の前記第２関数は、第２関数群の１つであり、前記ｉ番目の前記第３関数は、第３関数群の１つであり、前記ｉ番目の前記第１要素関数は、第１要素関数群の１つであり、
前記第３変数更新は、前記第３変数更新前の第３変数ｚに前記ｉ番目の第２要素関数を加えて前記ｉ番目の前記第３変数ｚを更新することを含み、前記ｉ番目の前記第３変数ｚは、第３変数群｛ｚ｝の１つであり、前記第２要素関数の変数は、前記第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含み、
前記処理部は、少なくとも、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉ、及び、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉの関数の少なくともいずれかを出力する、計算方法。

実施形態によれば、最適化問題を高速に計算できる計算装置、計算プログラム、記録媒体及び計算方法が提供できる。

以上、例を参照しつつ、本発明の実施の形態について説明した。しかし、本発明は、これらの例に限定されるものではない。例えば、計算装置に含まれる処理部、取得部及び保持部などの各要素の具体的な構成に関しては、当業者が公知の範囲から適宜選択することにより本発明を同様に実施し、同様の効果を得ることができる限り、本発明の範囲に包含される。

各例のいずれか２つ以上の要素を技術的に可能な範囲で組み合わせたものも、本発明の要旨を包含する限り本発明の範囲に含まれる。

本発明の実施の形態として上述した計算装置、計算プログラム、記録媒体及び計算方法を基にして、当業者が適宜設計変更して実施し得る全ての計算装置、計算プログラム、記録媒体及び計算方法も、本発明の要旨を包含する限り、本発明の範囲に属する。

本発明の思想の範疇において、当業者であれば、各種の変更例及び修正例に想到し得るものであり、それら変更例及び修正例についても本発明の範囲に属するものと了解される。

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。

１０…保持部、 １０Ａ、１０Ｂ…第１、第２保持領域、 １５Ａ〜１５Ｃ…第１〜第３回路、 １５Ｘ…制御回路部、 １６Ａ、１６Ｂ…第１、第２制御部、 １６Ｘ…回路、 ２０…処理部、 ２０Ａ〜２０Ｄ…第１〜第４計算部、 ３１…取得部、 ３２…操作部、 ３３…表示部、 １１０、１１１…計算装置、 １１０Ｅ、１１９Ｒ…計算結果、Ｍ１１、Ｍ１２…計算結果、 ＮＬ…処理手順の繰り返し回数、 Ｎｍｃ…カット数、 Ｔ…計算パラメータ、 Ｔｃ…計算時間、 Ｖａｖｅ…平均値

Claims (11)

1. 処理手順を繰り返す処理部を備え、
   前記処理手順は、第１変数更新、第２変数更新及び第３変数更新を含み、
   前記第１変数更新は、前記第１変数更新前の前記ｉ番目（ｉは１以上Ｎ以下の整数であり、前記Ｎは２以上の整数）の第１変数ｘ_ｉにｉ番目の第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することを含み、前記ｉ番目の前記第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つであり、前記ｉ番目の前記第１関数の変数は、前記ｉ番目の第２変数ｙ_ｉ含み、前記ｉ番目の前記第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つであり、前記ｉ番目の前記第１関数は、第１関数群の１つであり、
   前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉに、前記ｉ番目の第２関数、前記ｉ番目の第３関数及び前記ｉ番目の第１要素関数の演算結果を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、前記ｉ番目の前記第２関数の変数は、前記ｉ番目の前記第１変数ｘ_ｉを含み、前記ｉ番目の前記第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び前記第１変数群｛ｘ｝の少なくとも一部を含み、前記ｉ番目の前記第１要素関数の変数は、前記ｉ番目の前記第２変数ｙ_ｉと第３変数ｚとの積を含み、前記ｉ番目の前記第２関数は、第２関数群の１つであり、前記ｉ番目の前記第３関数は、第３関数群の１つであり、前記ｉ番目の前記第１要素関数は、第１要素関数群の１つであり、
   前記第３変数更新は、前記第３変数更新前の第３変数ｚに前記ｉ番目の第２要素関数を加えて前記第３変数ｚを更新することを含み、前記第３変数ｚは、第３変数群｛ｚ｝の１つであり、前記第２要素関数の変数は、前記第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含み、
   前記処理部は、少なくとも、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉ、及び、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉの関数の少なくともいずれかを出力する、計算装置。
2. 前記処理手順の後の前記第１計算パラメータは、前記処理手順の前の前記第１計算パラメータとは異なる、請求項１記載の計算装置。
3. 前記第２変数更新は、第１サブ更新及び第２サブ更新を含み、
   前記第１サブ更新は、前記第１サブ更新前の前記ｉ番目の前記第２変数ｙ_ｉに前記第２関数を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、
   前記第２サブ更新は、前記第２サブ更新前の前記ｉ番目の前記第２変数ｙ_ｉに前記第３関数を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、
   前記第１変数更新、前記第１サブ更新及び前記第３変数更新をＭ回（Ｍは２以上の整数）行った後、前記第２サブ更新を行う、または、前記第２サブ更新後、前記第１変数更新、前記第１サブ更新及び前記第３変数更新をＭ回行う、請求項１または２に記載の計算装置。
4. 前記処理部の１つの部分は、第１の更新、第１の算出、及び、第２の更新を実施し、前記第１の更新は、前記第２変数群{ｙ}の前記一部の更新を含み、前記第１の算出は、前記第２変数群{ｙ}の一部に関する２乗平均の算出を含み、前記第２の更新は、前記第３変数群｛ｚ｝に含まれる１つの前記第３変数ｚの前記２乗平均を用いた更新を含み、
   前記処理部の別の１つの部分は、第３の更新、第２の算出、及び、第４の更新を実施し、前記第３の更新は、前記第２変数群{ｙ}の別の一部の更新を含み、前記第２の算出は、第２変数群{ｙ}の前記別の一部に関する前記２乗平均の算出を含み、前記第４の更新は、前記第３変数群｛ｚ｝に含まれる別の１つの前記第３変数ｚの前記２乗平均を用いた更新を含み、
   前記第１の更新、前記第１の算出、及び、前記第２の更新の少なくとも一部の実施は、前記第３の更新、前記第２の算出、及び、前記第４の更新の少なくとも一部の実施と同時に行われる、請求項１〜３のいずれか１つに記載の計算装置。
5. 前記ｉ番目の前記第２関数は、前記ｉ番目の前記第１変数ｘ_ｉを変数とする前記ｉ番目の第４関数を含み、
   前記ｉ番目の第４関数は、前記第１変数ｘ_ｉの非線形関数を含み、
   前記ｉ番目の第４関数は、演算パラメータｐを含み、
   前記演算パラメータｐは、前記処理手順を繰り返すと変化し、
   前記処理手順を繰り返した後の前記ｉ番目の前記第４関数の実数根の数は、前記処理手順を繰り返す前の前記ｉ番目の前記第４関数の実数根の数とは異なる、請求項１〜４のいずれか１つに記載の計算装置。
6. 前記第１変数更新は、前記第１変数更新前の前記ｉ番目の前記第１変数ｘ_ｉを保持部から取得し、前記第１変数更新後の前記ｉ番目の前記第１変数ｘ_ｉを前記保持部に保持させることを含み、
   前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉを前記保持部から取得し、前記第２変数更新後の前記ｉ番目の前記第２変数ｙ_ｉを前記保持部に保持させることを含み、
   前記第３変数更新は、前記第３変数更新前の前記第３変数ｚを前記保持部から取得し、前記第３変数更新後の前記第３変数ｚを前記保持部に保持させることを含む、請求項１〜５のいずれか１つに記載の計算装置。
7. 前記処理部は、第１計算部と、第２計算部と、を含み、
   前記第１計算部は、前記第３関数の計算の一部を実施し、
   前記第２計算部は、前記第３関数の前記計算の別の一部を実施し、
   前記第１計算部における前記第３関数の前記計算の前記一部の前記実施の少なくとも一部と、前記第２計算部における前記第３関数の前記計算の前記別の一部の前記実施の少なくとも一部は、同時に行われる、請求項１〜６のいずれか１つに記載の計算装置。
8. 前記保持部は、第１保持領域及び第２保持領域を含み、
   前記第１パラメータ群｛Ｊ｝は、第１計算使用部分と、第２計算使用部分と、を含み、前記第１計算使用部分は、前記第３関数の前記計算の前記一部で用いられ、前記第２計算使用部分は、前記第３関数の前記計算の前記別の一部で用いられ、
   前記第１計算部は、前記第１計算使用部分を前記第１保持領域に保持させ、
   前記第２計算部は、前記第２計算使用部分を前記第２保持領域に保持させる、請求項６記載の計算装置。
9. コンピュータに、処理手順を繰り返させる計算プログラムであって、
   前記処理手順は、第１変数更新、第２変数更新及び第３変数更新を含み、
   前記第１変数更新は、前記第１変数更新前の前記ｉ番目（ｉは１以上Ｎ以下の整数であり、前記Ｎは２以上の整数）の第１変数ｘ_ｉに前記ｉ番目の第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することを含み、前記ｉ番目の前記第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つであり、前記ｉ番目の前記第１関数の変数は、前記ｉ番目の第２変数ｙ_ｉ含み、前記ｉ番目の前記第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つであり、前記ｉ番目の前記第１関数は、第１関数群の１つであり、
   前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉに、前記ｉ番目の第２関数、前記ｉ番目の第３関数及び前記ｉ番目の第１要素関数の演算結果を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、前記ｉ番目の前記第２関数の変数は、前記ｉ番目の前記第１変数ｘ_ｉを含み、前記ｉ番目の前記第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び前記第１変数群｛ｘ｝の少なくとも一部を含み、前記ｉ番目の前記第１要素関数の変数は、前記ｉ番目の前記第２変数ｙ_ｉと第３変数ｚとの積を含み、前記ｉ番目の前記第２関数は、第２関数群の１つであり、前記ｉ番目の前記第３関数は、第３関数群の１つであり、前記ｉ番目の前記第１要素関数は、第１要素関数群の１つであり、
   前記第３変数更新は、前記第３変数更新前の第３変数ｚに前記ｉ番目の第２要素関数を加えて前記ｉ番目の前記第３変数ｚを更新することを含み、前記ｉ番目の前記第３変数ｚは、第３変数群｛ｚ｝の１つであり、前記第２要素関数の変数は、前記第２変数群｛ｙ｝の少なくとも一部の２乗平均値と、第１計算パラメータと、の差を含み、
   前記処理部は、少なくとも、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉ、及び、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉの関数の少なくともいずれかを出力する、計算プログラム。
10. コンピュータに、処理手順を繰り返させる計算プログラムを記録したコンピュータ読み取り可能な記録媒体であって、
    前記処理手順は、第１変数更新、第２変数更新及び第３変数更新を含み、
    前記第１変数更新は、前記第１変数更新前の前記ｉ番目（ｉは１以上Ｎ以下の整数であり、前記Ｎは２以上の整数）の第１変数ｘ_ｉに前記ｉ番目の第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することを含み、前記ｉ番目の前記第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つであり、前記ｉ番目の前記第１関数の変数は、前記ｉ番目の第２変数ｙ_ｉ含み、前記ｉ番目の前記第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つであり、前記ｉ番目の前記第１関数は、第１関数群の１つであり、
    前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉに、前記ｉ番目の第２関数、前記ｉ番目の第３関数及び前記ｉ番目の第１要素関数の演算結果を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、前記ｉ番目の前記第２関数の変数は、前記ｉ番目の前記第１変数ｘ_ｉを含み、前記ｉ番目の前記第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び前記第１変数群｛ｘ｝の少なくとも一部を含み、前記ｉ番目の前記第１要素関数の変数は、前記ｉ番目の前記第２変数ｙ_ｉと第３変数ｚとの積を含み、前記ｉ番目の前記第２関数は、第２関数群の１つであり、前記ｉ番目の前記第３関数は、第３関数群の１つであり、前記ｉ番目の前記第１要素関数は、第１要素関数群の１つであり、
    前記第３変数更新は、前記第３変数更新前の第３変数ｚに前記ｉ番目の第２要素関数を加えて前記ｉ番目の前記第３変数ｚを更新することを含み、前記ｉ番目の前記第３変数ｚは、第３変数群｛ｚ｝の１つであり、前記第２要素関数の変数は、前記第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含み、
    前記処理部は、少なくとも、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉ、及び、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉの関数の少なくともいずれかを出力する、記録媒体。
11. 処理手順を繰り返し、
    前記処理手順は、第１変数更新、第２変数更新及び第３変数更新を含み、
    前記第１変数更新は、前記第１変数更新前の前記ｉ番目（ｉは１以上Ｎ以下の整数であり、前記Ｎは２以上の整数）の第１変数ｘ_ｉに前記ｉ番目の第１関数を加えて前記ｉ番目の前記第１変数ｘ_ｉを更新することを含み、前記ｉ番目の前記第１変数ｘ_ｉは、第１変数群｛ｘ｝の１つであり、前記ｉ番目の前記第１関数の変数は、前記ｉ番目の第２変数ｙ_ｉ含み、前記ｉ番目の前記第２変数ｙ_ｉは、第２変数群｛ｙ｝の１つであり、前記ｉ番目の前記第１関数は、第１関数群の１つであり、
    前記第２変数更新は、前記第２変数更新前の前記ｉ番目の前記第２変数ｙ_ｉに、前記ｉ番目の第２関数、前記ｉ番目の第３関数及び前記ｉ番目の第１要素関数の演算結果を加えて前記ｉ番目の前記第２変数ｙ_ｉを更新することを含み、前記ｉ番目の前記第２関数の変数は、前記ｉ番目の前記第１変数ｘ_ｉを含み、前記ｉ番目の前記第３関数の変数は、第１パラメータ群｛Ｊ｝の少なくとも一部及び前記第１変数群｛ｘ｝の少なくとも一部を含み、前記ｉ番目の前記第１要素関数の変数は、前記ｉ番目の前記第２変数ｙ_ｉと第３変数ｚとの積を含み、前記ｉ番目の前記第２関数は、第２関数群の１つであり、前記ｉ番目の前記第３関数は、第３関数群の１つであり、前記ｉ番目の前記第１要素関数は、第１要素関数群の１つであり、
    前記第３変数更新は、前記第３変数更新前の第３変数ｚに前記ｉ番目の第２要素関数を加えて前記ｉ番目の前記第３変数ｚを更新することを含み、前記ｉ番目の前記第３変数ｚは、第３変数群｛ｚ｝の１つであり、前記第２要素関数の変数は、前記第２変数群｛ｙ｝の少なくとも一部の２乗平均と、第１計算パラメータと、の差を含み、
    前記処理部は、少なくとも、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉ、及び、前記処理手順を前記繰り返した後に得られる前記ｉ番目の前記第１変数ｘ_ｉの関数の少なくともいずれかを出力する、計算方法。

Images