4G以内的质数表生成程序仅主程序MATLAB代码.pdf

function cp() time_slice=clock; % To record the time that process started. % Author : 孙乙方 % Concise Prime % Code Aim: Fast % Time Spent on This Program: 56 Hours % Time Limit: 60 Hours % Techniques % In IF syntax, 0 is faster than ~. % When using counters, FOR loop is much more faster than WHILE loop. % To Matlab Default Type DOUBLE (which is a dynamic file space consumption % one), within range [-2^53,2^53], the best resolution is 1, on X32 and % X64 machines. % To FOR loop, the looping variable is DOUBLE type in INT32 range on X32 machines and % in INT64 range on X64 machines. % When x and y are both DOUBLE type, mod(x,y)=x-floor(x/y)*y, and the % latter one is faster. % Theorems % Sieve of Eratosthenes % Between n and 2*n, there must be one prime number. (n=1) % Technologies % 1-SuperConcise Technology, to make the code consise, in order to reduce % errors and improve efficiency. % 2-PreRetraction Technology, to retract the table to 2 and odd numbers % before the computation begins to reduce time of process. % 3-TypeDouble Technology, to use the DOUBLE datatype of Matlab to % improve efficiency. % 4-ConstantVM Technology, to use virtual memory file constantly to % reduce code complexity. % 5-QuickFile Technology, to use FWRITE function instead of SAVE % function to improve file operation speed. % 6-AvatarSwitch Technology, to use file as a switch instead of copying % and deleting files. % 7-SmartCounter Technology, to use smart counter variable to count the % amount of new prime numbers on variable b side. % 8-FileOnDemand Technology, to read data from file on demand to save % time. % Load Progress load(INIT,vm_amount,ite_finished,direction,number_for_check,pos_for_check,data_tail,data_len,top); % Variables G41=4294967295; % 4G-1 M12 % 12M % Prime Table Bou