JP7007585B2 - 最適化装置、最適化装置の制御方法及び最適化装置の制御プログラム - Google Patents
最適化装置、最適化装置の制御方法及び最適化装置の制御プログラム Download PDFInfo
- Publication number
- JP7007585B2 JP7007585B2 JP2018049467A JP2018049467A JP7007585B2 JP 7007585 B2 JP7007585 B2 JP 7007585B2 JP 2018049467 A JP2018049467 A JP 2018049467A JP 2018049467 A JP2018049467 A JP 2018049467A JP 7007585 B2 JP7007585 B2 JP 7007585B2
- Authority
- JP
- Japan
- Prior art keywords
- unit
- operating conditions
- arithmetic
- optimization device
- ising model
- Prior art date
- Legal status (The legal status 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 status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/047—Probabilistic or stochastic networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/01—Dynamic search techniques; Heuristics; Dynamic trees; Branch-and-bound
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
-
- 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
- G06F7/58—Random or pseudo-random number generators
- G06F7/588—Random number generators, i.e. based on natural stochastic processes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Algebra (AREA)
- Probability & Statistics with Applications (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
1つの側面では、本発明は、問題に応じた適切な動作条件を自動的に設定可能な最適化装置、最適化装置の制御方法及び最適化装置の制御プログラムを提供することを目的とする。
(第1の実施の形態)
以下に示す最適化装置は、イジングモデルの基底状態(イジング型のエネルギー関数が最小値となるときの各ニューロンの値)を探索することで、巡回セールスマン問題などの最適化問題を計算する。イジング型のエネルギー関数E(x)は、たとえば、以下の式(1)で定義される。
また、i番目のニューロンのスピン反転(値の変化)に伴うエネルギー変化ΔEiは、以下の式(2)で表される。
最適化装置10は、操作部11、管理部12、演算部13、記憶部14を有する。
操作部11には、計算対象の問題が入力される。操作部11には、たとえば、図示しない入力デバイスが接続され、ユーザによる入力デバイスの操作によって、問題が入力される。
管理部12は、たとえば、CPU(Central Processing Unit)やDSP(Digital Signal Processor)などの演算処理装置としてのプロセッサである。ただし、管理部12は、ASIC(Application Specific Integrated Circuit)やFPGA(Field Programmable Gate Array)などの特定用途の電子回路を含んでもよい。なお、複数のプロセッサの集合を「マルチプロセッサ」または単に「プロセッサ」ということがある。
演算部13は、入力された初期動作条件または全体動作条件を用いて、イジングモデルの基底状態を探索する。演算部13は、確率的探索法により、イジングモデルの基底状態を探索する。
計算対象の問題が操作部11に入力されると、管理部12は、その問題をイジングモデルに変換し、初期動作条件とともに演算部13に入力する。なお、管理部12は、初期動作条件を予め演算部13に入力しておいてもよい。
演算部13は、全体動作条件を用いて、上記と同様にイジングモデルの基底状態を探索する。管理部12は、全体動作条件を用いた演算部13における探索結果に基づいて、その全体動作条件をさらに変更し、変更した全体動作条件を演算部13に入力して、演算部13に探索させてもよい。
第2の実施の形態の最適化装置は、交換モンテカルロ法(レプリカ交換法または拡張アンサンブル法などとも呼ばれる)を実現する。
たとえば、図2では、複数(4つ)のレプリカが用意されており、各レプリカの温度はT1~T4(T4>T3>T2>T1)であるとする。交換モンテカルロ法では、温度が隣接して設定されているレプリカ同士で、所定の交換確率にしたがって、ノードの状態(または温度)が交換される。交換確率はmin(1,R)とすることができ、Rは以下の式(4)で表される。
たとえば、温度がT1のレプリカと温度がT2のレプリカ間では、温度がT2のレプリカのエネルギーが、温度がT1のレプリカのエネルギーよりも大きいときには、確率Rでノードの状態の交換が行われる。
最適化装置20は、操作部21と、管理部22と、演算部23a1,23a2,・・・,23anを含む演算部23、記憶部24と、を有する。
管理部22は、イジングモデル変換部22a、交換制御部22bを有する。
イジングモデル変換部22aは、問題に応じて、たとえば、式(1)で定義されるイジング型のエネルギー関数E(x)の重み係数Wijとバイアスbiを求めることで、問題のイジングモデルへの変換を行う。
図4は、第2の実施の形態の最適化装置の動作の一例の流れを示すフローチャートである。
そして、管理部22は、初期動作条件または全体動作条件を変更することで得られた動作条件を演算部23a1~23anに与えて、再びステップS13の探索を実行させる。
(第3の実施の形態)
図5は、第3の実施の形態の最適化装置の一例を示す図である。図5において、図3に示した第2の実施の形態の最適化装置20と同じ要素については同一符号が付されている。
演算部32に含まれる演算部32a1~32anのそれぞれは、第2の実施の形態の最適化装置20の演算部23a1~23anのそれぞれと同様の機能を有するとともに、上記アニーリング条件に基づいて、シミュレーテッド・アニーリングを行う。
図6は、第3の実施の形態の最適化装置の動作の一例の流れを示すフローチャートである。
図7は、アニーリング条件の変更例を示す図である。
図8は、第4の実施の形態の最適化装置の一例を示す図である。図8において、図3に示した第2の実施の形態の最適化装置20と同じ要素については同一符号が付されている。
第2管理部41bは、演算部42aへの初期動作条件の入力、及び探索結果に応じた動作条件の変更(初期動作条件から全体動作条件への変更など)を行う。つまり、第2管理部41bは、第2の実施の形態の最適化装置20の交換制御部22b、または、第3の実施の形態の最適化装置30のアニーリング条件設定部31aと同様の、動作条件設定(または変更)機能を有する。
上記のような最適化装置40においても、第2及び第3の実施の形態の最適化装置20,30と同様の効果が得られるとともに、演算装置42内のハードウェアにより動作条件を変更することが可能になるため、変更にかかる時間を短縮でき計算時間を短縮できる。
ところで、上記各実施の形態において操作部11,21には、問題の代りに、その問題を変換したイジングモデルがユーザによって入力されてもよい。その場合、管理部12,22,31,41は問題のイジングモデルへの変換を行わなくてもよい。ただ、管理部12,22,31,41が問題のイジングモデルへの変換を行うことで、ユーザは問題を入力するだけでよいので利便性が高くなる。
図9は、コンピュータのハードウェア例を示す図である。
コンピュータ50は、CPU51、RAM52、HDD53、画像信号処理部54、入力信号処理部55、媒体リーダ56及び通信インタフェース57を有する。上記ユニットは、バスに接続されている。
プログラムは、コンピュータ読み取り可能な記録媒体(たとえば、記録媒体56a)に記録しておくことができる。記録媒体として、たとえば、磁気ディスク、光ディスク、光磁気ディスク、半導体メモリなどを使用できる。磁気ディスクには、FDおよびHDDが含まれる。光ディスクには、CD、CD-R(Recordable)/RW(Rewritable)、DVDおよびDVD-R/RWが含まれる。プログラムは、可搬型の記録媒体に記録されて配布されることがある。その場合、可搬型の記録媒体から他の記録媒体(たとえば、HDD53)にプログラムをコピーして実行してもよい。
11 操作部
12 管理部
13 演算部
14 記憶部
Claims (6)
- 問題が入力される操作部と、
イジングモデルの基底状態を探索する演算部と、
前記操作部から入力された前記問題を前記イジングモデルに変換し、変換した前記イジングモデルと初期動作条件とを前記演算部に入力し、前記初期動作条件を用いて前記演算部に前記基底状態を探索させた結果に基づいて、前記初期動作条件を変更した全体動作条件を用いて前記基底状態を前記演算部に探索させる、管理部と、
を有する最適化装置。 - 前記演算部は、複数の第1の演算部を有し、
前記管理部は、前記複数の第1の演算部のそれぞれに、互いに異なる前記初期動作条件を入力し、前記複数の第1の演算部のそれぞれにおいて前記基底状態を探索させた結果に基づいた交換確率にしたがって、前記複数の第1の演算部のうち、第2の演算部と第3の演算部に入力する前記初期動作条件を交換する、
請求項1に記載の最適化装置。 - 前記演算部は、複数の第1の演算部を有し、
前記管理部は、前記複数の第1の演算部のそれぞれに、互いに異なる前記初期動作条件を入力し、前記複数の第1の演算部のそれぞれにおいて前記基底状態を探索させた結果に基づいて、前記複数の第1の演算部のそれぞれに、前記初期動作条件を変更した、互いに異なる前記全体動作条件を入力する、
請求項1に記載の最適化装置。 - 前記管理部は、前記問題を前記イジングモデルに変換する第1の管理部と、前記演算部が搭載された演算装置に実装され、前記初期動作条件を前記全体動作条件に変更する第2の管理部とを有する、
請求項1乃至3の何れか一項に記載の最適化装置。 - 問題が入力される操作部と、イジングモデルの基底状態を探索する演算部と、前記演算部を管理する管理部とを有する最適化装置の制御方法において、
前記管理部が、
前記操作部から入力された前記問題を前記イジングモデルに変換し、
変換した前記イジングモデルと初期動作条件とを前記演算部に入力し、
前記初期動作条件を用いて前記演算部に前記基底状態を探索させた結果に基づいて、前記初期動作条件を変更した全体動作条件を用いて前記基底状態を前記演算部に探索させる、
最適化装置の制御方法。 - 問題が入力される操作部と、イジングモデルの基底状態を探索する演算部と、前記演算部を管理する管理部とを有する最適化装置の制御プログラムにおいて、
前記管理部に、
前記操作部から入力された前記問題を前記イジングモデルに変換し、
変換した前記イジングモデルと初期動作条件とを前記演算部に入力し、
前記初期動作条件を用いて前記演算部に前記基底状態を探索させた結果に基づいて、前記初期動作条件を変更した全体動作条件を用いて前記基底状態を前記演算部に探索させる、
最適化装置の制御プログラム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018049467A JP7007585B2 (ja) | 2018-03-16 | 2018-03-16 | 最適化装置、最適化装置の制御方法及び最適化装置の制御プログラム |
US16/284,101 US11599073B2 (en) | 2018-03-16 | 2019-02-25 | Optimization apparatus and control method for optimization apparatus using ising models |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018049467A JP7007585B2 (ja) | 2018-03-16 | 2018-03-16 | 最適化装置、最適化装置の制御方法及び最適化装置の制御プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2019160169A JP2019160169A (ja) | 2019-09-19 |
JP7007585B2 true JP7007585B2 (ja) | 2022-01-24 |
Family
ID=67905541
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2018049467A Active JP7007585B2 (ja) | 2018-03-16 | 2018-03-16 | 最適化装置、最適化装置の制御方法及び最適化装置の制御プログラム |
Country Status (2)
Country | Link |
---|---|
US (1) | US11599073B2 (ja) |
JP (1) | JP7007585B2 (ja) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7248907B2 (ja) | 2019-08-14 | 2023-03-30 | 富士通株式会社 | 最適化装置および最適化装置の制御方法 |
JP7339539B2 (ja) * | 2020-01-15 | 2023-09-06 | 富士通株式会社 | 最適化装置、最適化装置の温度設定方法及び最適化装置の温度設定プログラム |
JP7417074B2 (ja) | 2020-02-19 | 2024-01-18 | 富士通株式会社 | 最適化装置、最適化方法及び最適化装置の制御プログラム |
JP7417078B2 (ja) | 2020-02-27 | 2024-01-18 | 富士通株式会社 | 情報処理装置、情報処理方法およびプログラム |
JP2021149808A (ja) * | 2020-03-23 | 2021-09-27 | 富士通株式会社 | Cpu状態表示方法及びcpu状態表示プログラム |
JP2021182338A (ja) * | 2020-05-20 | 2021-11-25 | 富士通株式会社 | 最適化装置、最適化プログラムおよび最適化方法 |
JP2021184148A (ja) | 2020-05-21 | 2021-12-02 | 富士通株式会社 | 最適化装置、最適化方法、および最適化プログラム |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016103282A (ja) | 2015-12-25 | 2016-06-02 | 株式会社日立製作所 | 情報処理システム及び管理装置 |
JP2018005541A (ja) | 2016-07-01 | 2018-01-11 | 富士通株式会社 | 情報処理装置、イジング装置及び情報処理装置の制御方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030121028A1 (en) | 2001-12-22 | 2003-06-26 | Michael Coury | Quantum computing integrated development environment |
WO2014192153A1 (ja) * | 2013-05-31 | 2014-12-04 | 株式会社日立製作所 | 半導体装置 |
CN107077642B (zh) * | 2014-08-22 | 2021-04-06 | D-波系统公司 | 可用于量子计算的用于求解问题的系统和方法 |
JP5865457B1 (ja) * | 2014-08-29 | 2016-02-17 | 株式会社日立製作所 | 情報処理システム及び管理装置 |
JP6659957B2 (ja) | 2016-06-06 | 2020-03-04 | 富士通株式会社 | 情報処理装置、イジング装置及び情報処理装置の制御方法 |
-
2018
- 2018-03-16 JP JP2018049467A patent/JP7007585B2/ja active Active
-
2019
- 2019-02-25 US US16/284,101 patent/US11599073B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016103282A (ja) | 2015-12-25 | 2016-06-02 | 株式会社日立製作所 | 情報処理システム及び管理装置 |
JP2018005541A (ja) | 2016-07-01 | 2018-01-11 | 富士通株式会社 | 情報処理装置、イジング装置及び情報処理装置の制御方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2019160169A (ja) | 2019-09-19 |
US11599073B2 (en) | 2023-03-07 |
US20190286077A1 (en) | 2019-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7007585B2 (ja) | 最適化装置、最適化装置の制御方法及び最適化装置の制御プログラム | |
JP7071638B2 (ja) | 最適化装置、最適化装置の制御方法及び最適化装置の制御プログラム | |
JP2020046718A (ja) | 最適化装置、最適化装置の制御方法及び最適化装置の制御プログラム | |
JP7093009B2 (ja) | 最適化装置、最適化装置の制御方法及び最適化装置の制御プログラム | |
JP7100254B2 (ja) | 最適化システム、最適化システムの制御方法及び最適化システムの制御プログラム | |
JP7068575B2 (ja) | 最適化システム、最適化装置及び最適化システムの制御方法 | |
JP7053995B2 (ja) | 最適化装置及び最適化装置の制御方法 | |
JP7248907B2 (ja) | 最適化装置および最適化装置の制御方法 | |
US20210065087A1 (en) | Information processing apparatus, combination optimization method, and computer-readable recording medium recording combination optimization program | |
JP7417074B2 (ja) | 最適化装置、最適化方法及び最適化装置の制御プログラム | |
JP2020205049A (ja) | 最適化装置、最適化装置の制御方法及び最適化装置の制御プログラム | |
JP2021168096A (ja) | サンプリング装置、サンプリング方法及びサンプリングプログラム | |
CA3188066A1 (en) | Methods and systems for hyperparameter tuning and benchmarking | |
JP7219402B2 (ja) | 最適化装置、最適化装置の制御方法及び最適化装置の制御プログラム | |
JP7181454B2 (ja) | 最適化装置、最適化装置の制御方法および最適化装置の制御プログラム | |
JP6997404B2 (ja) | 最適化装置、最適化装置の制御方法および最適化装置の制御プログラム | |
JP2020187453A (ja) | 最適化装置および最適化装置の制御方法 | |
JP2020173579A (ja) | 最適化システムおよび最適化システムの制御方法 | |
JP2022094510A (ja) | 最適化プログラム、最適化方法及び情報処理装置 | |
JP2023061477A (ja) | プログラム、データ処理方法及びデータ処理装置 | |
JP2022072685A (ja) | 評価関数生成プログラム、評価関数生成方法、最適化方法及び最適化装置 | |
JP2021131723A (ja) | 情報処理方法、情報処理装置及びプログラム | |
US20240111833A1 (en) | Data processing apparatus and data processing method | |
US20220335487A1 (en) | Non-transitory computer-readable recording medium, data processing method, and data processing apparatus | |
JP2022161128A (ja) | プログラム、データ処理方法及びデータ処理装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A80 | Written request to apply exceptions to lack of novelty of invention |
Free format text: JAPANESE INTERMEDIATE CODE: A80 Effective date: 20180329 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20201210 |
|
RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20201217 |
|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20201217 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20211126 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20211207 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20211220 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7007585 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |