def multi_denss(niter, **kwargs):
try:
time.sleep(1)
# Processing keyword args for compatibility with RAW GUI
kwargs['path'] = '.'
kwargs['output'] = kwargs['output'] + '_' + str(niter)
np.random.seed(niter + int(time.time()))
kwargs['seed'] = np.random.randint(2**31 - 1)
kwargs['quiet'] = True
if niter <= superargs.nmaps - 1:
sys.stdout.write("\r Running denss job: %i / %i " %
(niter + 1, superargs.nmaps))
sys.stdout.flush()
fname = kwargs['output'] + '.log'
logger = logging.getLogger(kwargs['output'])
logger.setLevel(logging.INFO)
fh = logging.FileHandler(fname)
formatter = logging.Formatter('%(asctime)s - %(message)s')
fh.setFormatter(formatter)
logger.addHandler(fh)
kwargs['my_logger'] = logger
logger.info('BEGIN')
logger.info('Script name: %s', sys.argv[0])
logger.info('DENSS Version: %s', __version__)
logger.info('Data filename: %s', superargs.file)
logger.info('Output prefix: %s', kwargs['output'])
logger.info('Mode: %s', superargs.mode)
result = saxs.denss(**kwargs)
logger.info('END')
return result
except KeyboardInterrupt:
print("KeyboardInterrupt")
pass
qdata, Idata, sigqdata, qbinsc, Imean, chis, rg, supportV, rho, side = saxs.denss(
q=args.q,
I=args.I,
sigq=args.sigq,
dmax=args.dmax,
ne=args.ne,
voxel=args.voxel,
oversampling=args.oversampling,
rho_start=args.rho_start,
add_noise=args.add_noise,
limit_dmax=args.limit_dmax,
limit_dmax_steps=args.limit_dmax_steps,
recenter=args.recenter,
recenter_steps=args.recenter_steps,
recenter_mode=args.recenter_mode,
positivity=args.positivity,
extrapolate=args.extrapolate,
output=args.output,
steps=args.steps,
ncs=args.ncs,
ncs_steps=args.ncs_steps,
ncs_axis=args.ncs_axis,
seed=args.seed,
shrinkwrap=args.shrinkwrap,
shrinkwrap_old_method=args.shrinkwrap_old_method,
shrinkwrap_sigma_start=args.shrinkwrap_sigma_start,
shrinkwrap_sigma_end=args.shrinkwrap_sigma_end,
shrinkwrap_sigma_decay=args.shrinkwrap_sigma_decay,
shrinkwrap_threshold_fraction=args.shrinkwrap_threshold_fraction,
shrinkwrap_iter=args.shrinkwrap_iter,
shrinkwrap_minstep=args.shrinkwrap_minstep,
chi_end_fraction=args.chi_end_fraction,
write_xplor_format=args.write_xplor_format,
write_freq=args.write_freq,
enforce_connectivity=args.enforce_connectivity,
enforce_connectivity_steps=args.enforce_connectivity_steps,
cutout=args.cutout,
quiet=args.quiet,
DENSS_GPU=args.DENSS_GPU,
my_logger=my_logger)
qdata, Idata, sigqdata, qbinsc, Imean, chis, rg, supportV, rho, side = saxs.denss(
q=q,
I=I,
sigq=sigq,
dmax=args.dmax,
ne=args.ne,
voxel=args.voxel,
oversampling=args.oversampling,
limit_dmax=args.limit_dmax,
limit_dmax_steps=args.limit_dmax_steps,
recenter=args.recenter,
recenter_steps=args.recenter_steps,
recenter_mode=args.recenter_mode,
positivity=args.positivity,
flatten_low_density=args.flatten_low_density,
minimum_density=args.minimum_density,
maximum_density=args.maximum_density,
rho_start=args.rho_start,
extrapolate=args.extrapolate,
output=args.output,
steps=args.steps,
ncs=args.ncs,
ncs_steps=args.ncs_steps,
ncs_axis=args.ncs_axis,
seed=args.seed,
shrinkwrap=args.shrinkwrap,
shrinkwrap_sigma_start=args.shrinkwrap_sigma_start,
shrinkwrap_sigma_end=args.shrinkwrap_sigma_end,
shrinkwrap_sigma_decay=args.shrinkwrap_sigma_decay,
shrinkwrap_threshold_fraction=args.shrinkwrap_threshold_fraction,
shrinkwrap_iter=args.shrinkwrap_iter,
shrinkwrap_minstep=args.shrinkwrap_minstep,
chi_end_fraction=args.chi_end_fraction,
write_xplor_format=args.write_xplor_format,
write_freq=args.write_freq,
enforce_connectivity=args.enforce_connectivity,
enforce_connectivity_steps=args.enforce_connectivity_steps,
cutout=args.cutout)