New Perlin run with twice the map sizes.

scripts/spatial_maps/simplex_input_tuning_correlation/correlation_length_distance_perlin.py

@@ -19,7 +19,7 @@ def filter_run_names_by_par_dict(traj, par_dict):
     return run_name_list
 
 
-def correlation_length_fit_dict(traj, load=True):
+def correlation_length_fit_dict(traj, map_type='perlin', load=True):
     def phi_diff(d_phi):
         return (d_phi + np.pi) % (2 * np.pi) - np.pi
 
@@ -106,8 +106,8 @@ def correlation_length_fit_dict(traj, load=True):
 
             fit_correlation_lengths_dict[corr_len] = np.mean(fit_corr_lens_different_seeds)
 
-        np.save(DATA_FOLDER + 'fit_correlation_lengths_dict.npy', fit_correlation_lengths_dict)
+        np.save(DATA_FOLDER + map_type + '_fit_correlation_lengths_dict.npy', fit_correlation_lengths_dict)
         return fit_correlation_lengths_dict
     else:
-        fit_correlation_lengths_dict = np.load(DATA_FOLDER + 'fit_correlation_lengths_dict.npy')
+        fit_correlation_lengths_dict = np.load(DATA_FOLDER + map_type + '_fit_correlation_lengths_dict.npy')
         return fit_correlation_lengths_dict.item()

scripts/spatial_network/orientation_map/paper_figures_orientation_map.py

@@ -1269,7 +1269,7 @@ def plot_exc_and_inh_hdi_over_fit_corr_len(traj, plot_run_names):
     # linestyles = ['solid', 'dashed', 'solid']
     # markers = [verts, '', 'o']
 
-    corr_len_fit_dict = correlation_length_fit_dict(traj, load=True)
+    corr_len_fit_dict = correlation_length_fit_dict(traj, map_type='pinwheel', load=True)
 
     width = 2 * panel_size
     height = 1.2 * panel_size
@@ -1548,10 +1548,10 @@ if __name__ == "__main__":
     map_seed = 1
     exemplary_head_direction = 0
 
-    # corr_len_fit_dict = correlation_length_fit_dict(traj, load=False)
-    # plt.plot(corr_len_fit_dict.keys(),corr_len_fit_dict.values())
-    # plt.show()
-    # abbrechen
+    corr_len_fit_dict = correlation_length_fit_dict(traj, map_type='pinwheel', load=True)
+    plt.plot(corr_len_fit_dict.keys(),corr_len_fit_dict.values())
+    plt.show()
+    abbrechen
 
 
     print("## Map specifications")
@@ -1603,7 +1603,7 @@ if __name__ == "__main__":
         selected_inhibitory_neuron,
         cut_off_dist=100.)
     #
-    number_of_suggestions = 2
+    number_of_suggestions = 0
     representative_excitatory_neuron_indices = get_neurons_with_given_hdi(hdi_means["polar_exc"], hdi_means[
         "circular_exc"],
                                                                           number_of_suggestions, plot_run_names,

scripts/spatial_network/orientation_map/run_orientation_map.py

@@ -265,7 +265,7 @@ def main():
     traj.f_add_parameter("simulation.dt", 0.1, comment="Network simulation time step in ms")
     traj.f_add_parameter("simulation.duration", 1000, comment="Network simulation duration in ms")
 
-    correlation_length_range = np.linspace(1.0, 800.0, 12, endpoint=True).tolist()
+    correlation_length_range = np.linspace(1.0, 650.0, 24, endpoint=True).tolist()
     # correlation_length_range = [200.0]
     seed_range = range(10)
     # seed_range = [1]