def main():
if clean_dir:
clean_outdir(outdir)
x0 = generate_vf_landmarks(magnification, ecc0, ecc1, n_prototypes, vf='full', filename=full_filename('vf', outdir))
b0 = generate_V2_boundary_retinotopy(boundary_len, map_h, branch='vd', filename=full_filename('boundary', outdir), tilt=10.0)
y = initialize_map(map_h, map_w, x0, filename=full_filename('init', outdir), random_seed = random_seed)
optimize(x0, b0, b1, b2, map_h, map_w, y,
k0 = k0,
kr = kr,
eta0 = eta0,
outdir = outdir,
save_interval = save_interval,
n = n_iterations)
def main(rseed, outputdir):
if clean_dir:
clean_outdir(outputdir)
x0 = generate_vf_landmarks(magnification,
ecc0,
ecc1,
n_prototypes,
vf='full',
filename=full_filename('vf', outputdir))
b0 = generate_V2_boundary_retinotopy(boundary_len,
map_h,
branch='vd',
filename=full_filename(
'boundary', outputdir),
tilt=10.0)
y = initialize_map(map_h,
map_w,
x0,
filename=full_filename('init', outdir),
random_seed=rseed)
# loop through b1 and b2
b1s = [0.03 * 1.6**i for i in range(0, 8)]
b2s = [0.03 * 1.6**i for i in range(0, 8)]
print("b1:", b1s)
print("b2:", b2s)
for b1 in b1s:
for b2 in b2s:
print("b1 = %10.4f \t b2 = %10.4f" % (b1, b2))
optimize(
x0,
b0,
b1,
b2,
map_h,
map_w,
y,
k0=k0,
kr=kr,
eta0=eta0,
outdir=outputdir,
save_interval=0, # disable saving intermediate maps
n=n_iterations,
report_iterations=False,
log_iterations=False)