public void Create()
{
var values = new int[] { 0, 1, 2 };
var spins = new Spins(values.Select(v => new Spin(v)));
var initInteraction = new Dense <double>(1);
var ising = OpenJij.MakeTransverseIsing <Dense <double> >(spins, initInteraction, 0.1, 2);
this.DisposeAndCheckDisposedState(initInteraction);
this.DisposeAndCheckDisposedState(ising);
}
private static void Main()
{
//generate dense graph with size N=5
const uint n = 500;
using var dense = new Dense <double>(n);
//set interactions
for (uint i = 0; i < n; i++)
{
for (uint j = 0; j < n; j++)
{
dense.J[i, j] = (i == j) ? 0 : -1;
}
}
//set local fields
for (uint i = 0; i < n; i++)
{
dense.H[i] = -1;
}
//generate random engine (mersenne twister)
using var randRngine = new StdMt19937(0x1234);
//create classical Ising system
using var system = OpenJij.MakeClassicalIsing(dense.GenSpin(randRngine), dense);
//generate schedule list
//from beta=0.1 to beta=100, 10 samples, 10 Monte Carlo step for each tempearature
using var scheduleList = OpenJij.MakeClassicalScheduleList(0.1, 100, 10, 200);
//do annealing (updater: SingleSpinFlip)
var stopWatch = new Stopwatch();
stopWatch.Start();
Algorithm <SingleSpinFlip <ClassicalIsing <Dense <double> > > > .Run(system, randRngine, scheduleList);
stopWatch.Stop();
double ticks = stopWatch.ElapsedTicks;
double time = (ticks / Stopwatch.Frequency) * 1000;
Console.WriteLine($"time {time}[ms]");
//show spins
Console.Write("The result spins are [");
foreach (var elem in Result.GetSolution(system))
{
Console.Write($"{elem} ");
}
Console.WriteLine("]");
}
public void MakeClassicalScheduleList()
{
var scheduleList = OpenJij.MakeClassicalScheduleList(0.1, 100, 10, 200);
this.DisposeAndCheckDisposedState(scheduleList);
}