Hybrid QAOA

.gitignore

@@ -329,4 +329,3 @@ __pycache__/
 *.btm.cs
 *.odx.cs
 *.xsd.cs
-

Build/build.ps1

@@ -40,6 +40,9 @@ Build-One 'publish' '../MachineLearning.sln'
 Write-Host "##[info]Build Jupyter magic library"
 Build-One 'publish' '../Magic.sln'
 
+Write-Host "##[info]Build QAOA library"
+Build-One 'publish' '../QAOA.sln'
+
 if (-not $all_ok) {
     throw "At least one test failed execution. Check the logs."
 }

Build/manifest.ps1

@@ -9,6 +9,7 @@
         "Microsoft.Quantum.Chemistry",
         "Microsoft.Quantum.Numerics",
         "Microsoft.Quantum.MachineLearning"
+        "Microsoft.Quantum.QAOA"
     );
     Assemblies = @(
         ".\Standard\src\bin\$Env:BUILD_CONFIGURATION\netstandard2.1\Microsoft.Quantum.Standard.dll",
@@ -17,5 +18,6 @@
         ".\Chemistry\src\DataModel\bin\$Env:BUILD_CONFIGURATION\netstandard2.1\Microsoft.Quantum.Chemistry.DataModel.dll",
         ".\Chemistry\src\Runtime\bin\$Env:BUILD_CONFIGURATION\netstandard2.1\Microsoft.Quantum.Chemistry.Runtime.dll",
         ".\Chemistry\src\Tools\bin\$Env:BUILD_CONFIGURATION\netcoreapp3.1\qdk-chem.dll"
+        ".\QAOA\src\bin\$Env:BUILD_CONFIGURATION\netcoreapp3.1\QAOA.dll"
     ) | ForEach-Object { Get-Item (Join-Path $PSScriptRoot ".." $_) };
 } | Write-Output;

Build/pack.ps1

@@ -72,6 +72,12 @@ if ($Env:ENABLE_PYTHON -eq "false") {
     Pack-Wheel '../Python/qsharp'
 }
 
+Write-Host "##[info]Pack QAOA library"
+Pack-One '../QAOA/src/QAOA.csproj'
+
+Write-Host "##[info]Pack QAOA library"
+Pack-One '../QAOA/src/QAOA.csproj'
+
 if (-not $all_ok) {
     throw "At least one test failed execution. Check the logs."
 }

Build/test.ps1

@@ -40,6 +40,15 @@ Test-One '../Chemistry/tests/JupyterTests/JupyterTests.csproj'
 Write-Host "##[info]Testing Numerics/tests/NumericsTests.csproj"
 Test-One '../Numerics/tests/NumericsTests.csproj'
 
+Write-Host "##[info]Testing QAOA/tests/QaoaTests/QaoaTests.csproj"
+Test-One '../QAOA/tests//QaoaTests/QaoaTests.csproj'
+
+Write-Host "##[info]Testing QAOA/tests/HybridQaoaTests/HybridQaoaTests.csproj"
+Test-One '../QAOA/tests/HybridQaoaTests/HybridQaoaTests.csproj'
+
+Write-Host "##[info]Testing QAOA/tests/JupyterTests/JupyterTests.csproj"
+Test-One '../QAOA/tests/JupyterTests/JupyterTests.csproj'
+
 if (-not $all_ok) {
     throw "At least one test failed execution. Check the logs."
 }

QAOA.sln

@@ -0,0 +1,53 @@
+
+Microsoft Visual Studio Solution File, Format Version 12.00
+# Visual Studio Version 16
+VisualStudioVersion = 16.0.29920.165
+MinimumVisualStudioVersion = 10.0.40219.1
+Project("{9A19103F-16F7-4668-BE54-9A1E7A4F7556}") = "QAOA", "QAOA\src\QAOA.csproj", "{159D83BF-56A2-41B7-951B-35292F5A9F23}"
+EndProject
+Project("{9A19103F-16F7-4668-BE54-9A1E7A4F7556}") = "HybridQaoaTests", "QAOA\tests\HybridQaoaTests\HybridQaoaTests.csproj", "{1623F038-210A-4BA2-8D4E-9AF611610829}"
+EndProject
+Project("{9A19103F-16F7-4668-BE54-9A1E7A4F7556}") = "JupyterTests", "QAOA\tests\JupyterTests\JupyterTests.csproj", "{3C8B4AF6-5495-4F23-B7B2-634E5A667837}"
+EndProject
+Project("{9A19103F-16F7-4668-BE54-9A1E7A4F7556}") = "QaoaTests", "QAOA\tests\QaoaTests\QaoaTests.csproj", "{208C0012-D62F-4F29-8738-394A62BBE842}"
+EndProject
+Project("{2150E333-8FDC-42A3-9474-1A3956D46DE8}") = "src", "src", "{AD73B607-0EF5-4C11-8832-22D584939780}"
+EndProject
+Project("{2150E333-8FDC-42A3-9474-1A3956D46DE8}") = "Tests", "Tests", "{80082123-4392-4F0E-8F87-26183C620693}"
+EndProject
+Global
+	GlobalSection(SolutionConfigurationPlatforms) = preSolution
+		Debug|Any CPU = Debug|Any CPU
+		Release|Any CPU = Release|Any CPU
+	EndGlobalSection
+	GlobalSection(ProjectConfigurationPlatforms) = postSolution
+		{159D83BF-56A2-41B7-951B-35292F5A9F23}.Debug|Any CPU.ActiveCfg = Debug|Any CPU
+		{159D83BF-56A2-41B7-951B-35292F5A9F23}.Debug|Any CPU.Build.0 = Debug|Any CPU
+		{159D83BF-56A2-41B7-951B-35292F5A9F23}.Release|Any CPU.ActiveCfg = Release|Any CPU
+		{159D83BF-56A2-41B7-951B-35292F5A9F23}.Release|Any CPU.Build.0 = Release|Any CPU
+		{1623F038-210A-4BA2-8D4E-9AF611610829}.Debug|Any CPU.ActiveCfg = Debug|Any CPU
+		{1623F038-210A-4BA2-8D4E-9AF611610829}.Debug|Any CPU.Build.0 = Debug|Any CPU
+		{1623F038-210A-4BA2-8D4E-9AF611610829}.Release|Any CPU.ActiveCfg = Release|Any CPU
+		{1623F038-210A-4BA2-8D4E-9AF611610829}.Release|Any CPU.Build.0 = Release|Any CPU
+		{3C8B4AF6-5495-4F23-B7B2-634E5A667837}.Debug|Any CPU.ActiveCfg = Debug|Any CPU
+		{3C8B4AF6-5495-4F23-B7B2-634E5A667837}.Debug|Any CPU.Build.0 = Debug|Any CPU
+		{3C8B4AF6-5495-4F23-B7B2-634E5A667837}.Release|Any CPU.ActiveCfg = Release|Any CPU
+		{3C8B4AF6-5495-4F23-B7B2-634E5A667837}.Release|Any CPU.Build.0 = Release|Any CPU
+		{208C0012-D62F-4F29-8738-394A62BBE842}.Debug|Any CPU.ActiveCfg = Debug|Any CPU
+		{208C0012-D62F-4F29-8738-394A62BBE842}.Debug|Any CPU.Build.0 = Debug|Any CPU
+		{208C0012-D62F-4F29-8738-394A62BBE842}.Release|Any CPU.ActiveCfg = Release|Any CPU
+		{208C0012-D62F-4F29-8738-394A62BBE842}.Release|Any CPU.Build.0 = Release|Any CPU
+	EndGlobalSection
+	GlobalSection(SolutionProperties) = preSolution
+		HideSolutionNode = FALSE
+	EndGlobalSection
+	GlobalSection(NestedProjects) = preSolution
+		{159D83BF-56A2-41B7-951B-35292F5A9F23} = {AD73B607-0EF5-4C11-8832-22D584939780}
+		{1623F038-210A-4BA2-8D4E-9AF611610829} = {80082123-4392-4F0E-8F87-26183C620693}
+		{3C8B4AF6-5495-4F23-B7B2-634E5A667837} = {80082123-4392-4F0E-8F87-26183C620693}
+		{208C0012-D62F-4F29-8738-394A62BBE842} = {80082123-4392-4F0E-8F87-26183C620693}
+	EndGlobalSection
+	GlobalSection(ExtensibilityGlobals) = postSolution
+		SolutionGuid = {09481F99-E1AC-44CE-A892-7331820661DF}
+	EndGlobalSection
+EndGlobal

QAOA/src/Examples.cs

@@ -0,0 +1,72 @@
+
+
+namespace Quantum.QAOA
+{
+    using Microsoft.Quantum.QAOA.QaoaHybrid;
+    using System;
+
+
+    class Examples
+    {
+        static void Main(string[] args)
+        {
+            //PARAMETERS
+            int numberOfIterations = 50;
+            int p = 3;
+            int numberOfRandomStartingPoints = 1;
+
+            //EXAMPLES
+
+            //Quantum Santa (http://quantumalgorithmzoo.org/traveling_santa/)
+            double[] dtx = { 0.619193, 0.742566, 0.060035, -1.568955, 0.045490 };
+            double[] dtz = { 3.182203, -1.139045, 0.221082, 0.537753, -0.417222 };
+            double[] segmentCosts = { 4.70, 9.09, 9.03, 5.70, 8.02, 1.71 };
+            double[] dh = { 4 * 20 - 0.5 * 4.7, 4 * 20 - 0.5 * 9.09, 4 * 20 - 0.5 * 9.03, 4 * 20 - 0.5 * 5.70, 4 * 20 - 0.5 * 8.02, 4 * 20 - 0.5 * 1.71 };
+            double[] dJ = { 40.0,40.0,20.0,40.0,40.0,40.0,
+                            40.0,40.0,40.0,20.0,40.0,40.0,
+                            40.0,40.0,40.0,40.0,40.0,40.0,
+                            40.0,40.0,40.0,40.0,40.0,40.0,
+                            40.0,40.0,40.0,40.0,40.0,20.0,
+                            40.0,40.0,40.0,40.0,40.0,40.0};
+            ProblemInstance quantumSanta = new ProblemInstance(dh, dJ);
+
+
+            //MaxCut (medium.com/mdr-inc/qaoa-maxcut-using-blueqat-aaf33038f46e)
+            dh = new Double[] { 0,0,0,0,0};
+            dJ = new Double[]{ 0,1,0,1,0,
+                               0,0,1,0,0,
+                               0,0,0,1,1,
+                               0,0,0,0,1,
+                               0,0,0,0,0};
+            ProblemInstance maxCut1 = new ProblemInstance(dh, dJ);
+
+
+            //Rigetti MaxCut unit tests
+            dh = new Double[]{-0.5,0,-1,0.5};
+            dJ = new Double[]{0,1,2,0,
+                              0,0,0.5,0,
+                              0,0,0,2.5,
+                              0,0,0,0};
+            ProblemInstance maxCut2 = new ProblemInstance(dh, dJ);
+
+
+            dh = new Double[] { 0.8, -0.5 };
+            dJ = new Double[]{ 0, -1,
+                               0, 0};
+            ProblemInstance maxCut3 = new ProblemInstance(dh, dJ);
+
+            dh = new Double[] {0, 0 };
+            dJ = new Double[]{ 0, 1,
+                               0, 0};
+            ProblemInstance maxCut4 = new ProblemInstance(dh, dJ);
+
+            //END EXAMPLES
+
+            HybridQaoa cop = new HybridQaoa(numberOfIterations, p, maxCut4, numberOfRandomStartingPoints, true);
+
+            OptimalSolution res = cop.RunOptimization();
+            Console.WriteLine(res.optimalVector);
+
+            }
+    }
+}

QAOA/src/Jupyter/HybridQaoaMagic.cs

@@ -0,0 +1,141 @@
+namespace QAOA.Jupyter
+{
+    using System;
+    using System.Threading.Tasks;
+    using Microsoft.Jupyter.Core;
+    using Microsoft.Quantum.QAOA.QaoaHybrid;
+    using Newtonsoft.Json;
+
+    public class HybridQaoaRunMagic : MagicSymbol
+    {
+
+        public HybridQaoaRunMagic()
+        {
+            this.Name = $"%qaoa.hybridqaoa.run";
+            this.Documentation = new Documentation() { Summary = "Runs a hybrid QAOA algorithm with a classical optimizer that chooses input angles. QAOA parameters are provided as a json" };
+            this.Kind = SymbolKind.Magic;
+            this.Execute = this.Run;
+        }
+
+        /// <summary>
+        /// List of arguments
+        /// </summary>
+        public class Arguments
+        {
+            /// <summary>
+            /// Number of iterations in the fidelity sampling.
+            /// </summary>
+            [JsonProperty(PropertyName = "number_of_iterations")]
+            public int NumberOfIterations { get; set; }
+
+            /// <summary>
+            /// Depth of a QAOA circuit.
+            /// </summary>
+            [JsonProperty(PropertyName = "p")]
+            public int p { get; set; }
+
+            /// <summary>
+            /// Description of a combinatorial problem to be solved.
+            /// </summary>
+            [JsonProperty(PropertyName = "problem_instance")]
+            public ProblemInstance ProblemInstance { get; set; }
+
+            /// <summary>
+            /// Number of random starting points in the angles parameters spaces.
+            /// </summary>
+            [JsonProperty(PropertyName = "number_of_random_starting_points")]
+            public int NumberOfRandomStartingPoints { get; set; } = 1;
+
+            /// <summary>
+            /// Flag whether optimization should be logged into a file.
+            /// </summary>
+            [JsonProperty(PropertyName = "should_log")]
+            public Boolean ShouldLog { get; set; } = false;
+
+            /// <summary>
+            /// Initial beta angles.
+            /// </summary>
+            [JsonProperty(PropertyName = "initial_beta")]
+            public Double[] InitialBeta { get; set; } = null;
+
+            /// <summary>
+            /// Initial gamma angles.
+            /// </summary>
+            [JsonProperty(PropertyName = "initial_gamma")]
+            public Double[] InitialGamma { get; set; } = null;
+        }
+
+        /// <summary>
+        /// Runs a hybrid QAOA.
+        /// </summary>
+        public async Task<ExecutionResult> Run(string input, IChannel channel)
+        {
+            if (string.IsNullOrWhiteSpace(input))
+            {
+                channel.Stderr("Please provide correct arguments");
+                return ExecuteStatus.Error.ToExecutionResult();
+            }
+
+            var args = JsonConvert.DeserializeObject<Arguments>(input);
+
+            HybridQaoa hybridQaoa = new HybridQaoa(args.NumberOfIterations, args.p, args.ProblemInstance, args.NumberOfRandomStartingPoints, args.ShouldLog, args.InitialBeta, args.InitialGamma);
+
+            return hybridQaoa.RunOptimization().ToExecutionResult();
+        }
+    }
+
+        public class HybridQaoaProblemInstanceMagic : MagicSymbol
+        {
+            public HybridQaoaProblemInstanceMagic()
+            {
+                this.Name = $"%qaoa.hybridqaoa.create.problem.instance";
+                this.Documentation = new Documentation() { Summary = "Prepares a problem instance objects that serves as one of arguments to %qaoa.hybridqaoa.run." };
+                this.Kind = SymbolKind.Magic;
+                this.Execute = this.Run;
+            }
+
+            /// <summary>
+            /// List of arguments
+            /// </summary>
+            public class Arguments
+            {
+                /// <summary>
+                /// Coefficents for one-local Hamiltonian terms.
+                /// </summary>
+                [JsonProperty(PropertyName = "one_local_hamiltonian_coefficients")]
+                public Double[] OneLocalHamiltonianCoefficients { get; set; }
+
+                /// <summary>
+                /// Coefficents for one-local Hamiltonian terms.
+                /// </summary>
+                [JsonProperty(PropertyName = "two_local_hamiltonian_coefficients")]
+                public Double[] TwoLocalHamiltonianCoefficients { get; set; }
+
+                /// <summary>
+                /// Depth of a QAOA circuit.
+                /// </summary>
+                [JsonProperty(PropertyName = "problem_size_in_bits")]
+                public int ProblemSizeInBits { get; set; }
+
+            }
+
+            /// <summary>
+            /// Prepares a ProblemInstance object for a hybrid QAOA.
+            /// </summary>
+            public async Task<ExecutionResult> Run(string input, IChannel channel)
+            {
+                if (string.IsNullOrWhiteSpace(input))
+                {
+                    channel.Stderr("Please provide correct arguments");
+                    return ExecuteStatus.Error.ToExecutionResult();
+                }
+
+                var args = JsonConvert.DeserializeObject<Arguments>(input);
+
+                ProblemInstance problemInstance = new ProblemInstance(args.OneLocalHamiltonianCoefficients, args.TwoLocalHamiltonianCoefficients);
+
+                return problemInstance.ToExecutionResult();
+            }
+
+        }
+}

QAOA/src/QAOA.csproj

@@ -0,0 +1,14 @@
+<Project Sdk="Microsoft.Quantum.Sdk/0.12.20072031">
+
+  <PropertyGroup>
+    <OutputType>Exe</OutputType>
+    <TargetFramework>netcoreapp3.1</TargetFramework>
+  </PropertyGroup>
+
+  <ItemGroup>
+    <PackageReference Include="Accord.Math" Version="3.8.0" />
+	<PackageReference Include="Microsoft.Jupyter.Core" Version="1.3.52077" />
+	<PackageReference Include="Newtonsoft.Json" Version="12.0.3" />
+  </ItemGroup>
+
+</Project>