int main() { system("cls"); do { // requests prime number from the user. printf("Enter a prime number: "); scanf("%d",&ptest); // checks if it is prime before assigning it to the 'p' variable. flag=prime(ptest); }while(flag==0); if(flag==1) p=ptest; do { // requests another prime number from the user. printf("\nEnter a second prime number: "); scanf("%d",&qtest); // checks if it is prime before assigning it to the 'q' variable. flag=prime(qtest); }while(flag==0); if(flag==1) q=qtest; // gets the message from the user. printf("\nEnter message: "); fflush(stdin); scanf("%[^\n]%*c",Msg); len=strlen(Msg); for(i=0;i!=len;i++) m[i]=Msg[i]; // calculates the modulus and phi. n=p*q; phi=(p-1)*(q-1); // then e and d. calcE(); printf("\nPossible values of e and d: "); for(i=0;i<j-1;i++) printf("\n%ld\t%ld",e[i],d[i]); encrypt(); decrypt(); getch(); return(0); }
int main() { Particle particle; std::ofstream ofs; ofs.open( "kep" ); for (int i=0;i<1000000;i++) { double x = particle.qx; double y = particle.qy; // for (int s=0;i<sympI;i++) { symp2(&particle, dT); /* particle.px+=calc(x, y) *symp_u[i]* dT; particle.py+=calc(y, x) *symp_u[i]* dT; particle.qx+=particle.px/M*symp_k[i]* dT; particle.qy+=particle.py/M*symp_k[i]* dT; */ // } // std::cout << "px = " << particle.px << ", py = " << particle.py << std::endl; // std::cout << "qx = " << particle.qx << ", px = " << particle.px << std::endl; if (i%500==0) { ofs << particle.qx << " " << particle.qy << " " << 0.0 << " " << calcE(particle) << " " << a_length(particle) << std::endl; } } ofs.close(); }
//Update static void update(void) { calcJD(); calcE(); field_defaultIncidence(Ez); calcMB(); calcH(); }
static inline void update(void) { calcJD(); calcE(); calcMB(); calcH(); field_defaultIncidence(Hz); }
static inline void update(void) { calcJD(); calcE(); scatteredWave(Ey, EPS_EY); calcMB(); calcH(); }
void reportElement_toFile(FILE* f,FILE* reso, bool J2, double theta, Particle pa) { double q = hypot(hypot(pa.qx, pa.qy), pa.qz); double p = hypot(hypot(pa.px, pa.py), pa.pz); double a = a_length(q, p); double Energy = 0; if (J2) { Energy = calcE_J2(q, p, pa); } else { Energy = calcE(q, p); } double h = h_moment(pa); double e = e_anomaly(a, h); double I = I_anomaly(h, pa); double OMEGA = OMEGA_anomaly(I, h, pa); fprintf(f, "%.4f %.4f %.4f %.4f %.4f %.4f %.4f %.4f %.4f\n", pa.qx, pa.qy, pa.qz, Energy, a, h, e, I, theta); fprintf(reso, "%f %f\n", a*(pa.qx/q*cos(theta)+pa.qy/q*sin(theta)), a*(pa.qy/q*cos(theta)-pa.qx/q*sin(theta))); }
//Update static void update(void) { // _CalcJD(); calcJD(); calcE(); // MPI_Barrier(MPI_COMM_WORLD); //ăăăȘă? scatteredWave(Ez, EPS_EZ); //planeWave(Ez, EPS_EZ); //Connection_SendRecvE(); Connection_ISend_IRecvE(); // _CalcMB(); calcMB(); calcH(); Connection_ISend_IRecvH(); //Connection_SendRecvH(); ntff(); }
static void update(){ calcE(); field_defaultIncidence(Ezx); field_defaultIncidence(Ezy); calcH(); }