Skip to content

Trouble with high sweep counts #4

@zmorrell

Description

@zmorrell

Hello,
Great package, I am very excited to work with it. I have been having some trouble using higher sweep counts (> ~2000). I have distilled the problem to a minimal example. I get a segmentation fault when I try to execute the following code:

from pysa.sa import Solver
import numpy as np

float_type = 'float32'
n_vars = 2 
problem = 2 * np.random.random((n_vars, n_vars)) - 1 
problem = (problem + problem.T) / 2 

solver = Solver(problem=problem, problem_type='ising', float_type=float_type)

solver.metropolis_update(num_sweeps = 2048)

I have done some initial testing and I can run up to 2033 sweeps without segmentation fault, though I don't know if this number is machine specific.

Is there something that I am doing wrong, or is there some limit on the number of sweeps?

Thank you in advance!

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions