origin = [0,0,0] N = 2 n_bins = 30 #Bins for the histogram plot time1 = timeit.default_timer() #Timing The Loop for i in range(iterations_of_numbers): vector_magnitude_list = [] for i in range(iterations_of_postitions): #This loop generates the spheres magnitude, vector = Obs.RandSphere(N) vector_magnitude_list.append(magnitude) #print(str(vector_magnitude_list)) N += 2 #How many the number of flagella increases by flagella_number_list.append(N) mean_length_list.append(Obs.mean_value(vector_magnitude_list, n_bins)) #Now we append the empty lists from before #Finally, write to a file for a,b in zip(mean_length_list,flagella_number_list): f.write("%s,%s\n" % (a,b)) f.close() time2 = timeit.default_timer() print("Time taken for code to run... in seconds :- " + str(time2 - time1)) #Finish Timing The Loop Obs.mean_length_plotter() #Finally, send it to be plotted Obs.log_length_plotter(N)
mean_mag_torque_list.append( scipy.linalg.norm( obs.torque(final_force_list, initial_force_list, N, torque_list))) initial_force_list = [] final_force_list = [] mean_force = sum(mean_mag_force_list) / (iterations_of_positions) mean_torque = sum(mean_mag_torque_list) / (iterations_of_positions) mean_force_list.append(mean_force) mean_torque_list.append(mean_torque) number_list.append(N) N += 1 #Writing to a file for a, b in zip(mean_torque_list, number_list): g.write("%s,%s\n" % (a, b)) g.close() for a, b in zip(mean_force_list, number_list): f.write("%s,%s\n" % (a, b)) f.close() #Plotting it all obs.mean_length_plotter() obs.mean_torque_plotter() #obs.log_length_plotter(N)