- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我有一个程序需要从文件中读取数据,但我不想给出文件名并希望将数据放入我的命令文件中,例如 nbody < input ???
我提供的代码(称为 nbody.c)是 N 体模拟的简化版本,它说明了上述理论。可以使用以下命令编译、链接代码并创建名为 nbody 的可执行文件。 gcc –lm nbody.c –o nbody
为了运行可执行文件,提供了一个名为 input 的示例数据输入文件作为起点。可以使用以下命令执行串行代码,该命令需要几秒钟才能运行完成。 nbody <输入
我的程序是;
#include <stdio.h>
#include <math.h>
int m = 1000;
void simulate(int m, int n, double delt, double eps, double x[], double y[], double z[],
int iter, double xn[], double yn[], double zn[]);
int main(int argc, char *argv[])
{
double x[m], y[m], z[m], xn[m], yn[m], zn[m];
double delt, eps, resultant;
int i, n, iter;
char ch, file_name[25];
FILE *fp, *ofp;
//file_name=argv[0];
//char *mode = "r";
printf("Enter the name of file you wish to see\n"); //instead of this, i wanted it reads automatically, like taking data from the file and calculate it.
gets(file_name);
fp = fopen(file_name, "r");
rewind(fp);
printf("The contents of %s file are : \n", file_name);
fscanf(fp,"%d %lf %lf", &n, &delt, &eps);
printf("n = %d, delta t = %lf and tolerance = %lf \n",n, delt, eps);
fclose(fp);
//while(ch=fgets(fp)!=EOF){
//printf("enter n, delta t and tolerance\n");
//scanf("%d %lf %lf", &n, &delt, &eps);
//printf("%c",ch)
//printf("enter n, delta t and tolerance\n");
//scanf("%d %lf %lf", &n, &delt, &eps);
simulate(m, n, delt, eps, x, y, z, iter, xn, yn, zn);
printf("n = %d, delta t = %lf and tolerance = %lf \n",n, delt, eps);
for (i=1; i <= n; i++) {
resultant = sqrt(x[i]*x[i]+y[i]*y[i]+z[i]*z[i]);
printf("%-5.7lf %-5.7lf %-5.7lf %-5.7lf\n", x[i], y[i], z[i], resultant);
}
return 0;
}
void simulate(int m, int n, double delt, double eps, double x[], double y[], double z[],
int iter, double xn[], double yn[], double zn[])
{
/*
Author C.Ierotheou
Aug 2009
Copyright University of Greenwich
routine to crudely simulate body-body interactions (inefficient method)
*/
double zero = 0.0, one = 1.0, twopi = 6.2831853071795864769252866;
int i, itest, j;
double aux,auy,delt2,delt3,eps2,fx,fy,fz,ftx,fty,ftz,step;
// return if n<5
if (n < 5) return;
//initialization
iter = 0;
delt2 = 0.5*delt*delt;
delt3 = delt2*delt2;
eps2 = eps*eps;
//initial distribution of points
auy = zero;
itest=n+1;
j = n+2;
fx = (double)(n);
x[1] = zero;
y[1] = zero;
z[1] = one;
for (i=2; i <= n; i++) {
step = (double)(i+i-j)/fx;
auy = fmod(auy+3.6/sqrt(itest*(one-step*step)),twopi);
aux = sin(auy);
x[i] = aux*step;
y[i] = aux*sin(acos(step));
z[i] = cos(auy);
}
// main iteration loop
Mainloop:
iter=iter+1;
itest=0;
for (i=1; i <= n; i++) {
// total sum of force vectors
fx=zero;
fy=zero;
fz=zero;
for (j=1; j<=n ; j++) {
if (j != i) {
aux=pow(x[i]-x[j],2.0)+pow(y[i]-y[j],2.0)+pow(z[i]-z[j],2.0);
aux=aux*sqrt(aux);
fx=fx+(x[i]-x[j])/aux;
fy=fy+(y[i]-y[j])/aux;
fz=fz+(z[i]-z[j])/aux;
}
}
// tangential component of force
aux=x[i]*fx+y[i]*fy+z[i]*fz;
ftx=fx-x[i]*aux;
fty=fy-y[i]*aux;
ftz=fz-z[i]*aux;
aux=ftx*ftx+fty*fty+ftz*ftz;
if (aux > eps2) {
itest=1;
aux=sqrt(one-aux*delt3);
xn[i]=x[i]*aux+ftx*delt2;
yn[i]=y[i]*aux+fty*delt2;
zn[i]=z[i]*aux+ftz*delt2;
}
}
for (i=1; i <= n; i++) {
x[i]=xn[i];
y[i]=yn[i];
z[i]=zn[i];
}
if (itest == 1) goto Mainloop;
}
最佳答案
你问:
In C, how to read data from file using x
您不会在 C 中执行此操作。这种功能由您的 shell/控制台/运行时环境提供。
如果您有一个程序x
,可以从stdin
读取其输入,并且您有一个文件y
,其中包含适用于x
,您使用:
x < y
来自外壳。
更新
我想我知道你的问题可能出在哪里。你有:
char ch, file_name[25];
FILE *fp, *ofp;
//file_name=argv[0];
//char *mode = "r";
printf("Enter the name of file you wish to see\n"); //instead of this, i wanted it reads automatically, like taking data from the file and calculate it.
gets(file_name);
fp = fopen(file_name, "r");
您希望从命令行传递输入文件的名称,而不是使用 gets
读取输入文件的名称。我希望我在这方面走在正确的道路上。
在这种情况下,您需要使用:
char* file_name = argv[1]; // argv[0] is the program name
// argv[1] is the first argument to the program
然后,您可以使用:
x y
其中x
是您的程序,y
包含x
所需的输入数据。
关于c - 在C中,如何使用x<y(y是文件名)从文件中读取数据,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/24457742/
问题故障解决记录 -- Java RMI Connection refused to host: x.x.x.x .... 在学习JavaRMI时,我遇到了以下情况 问题原因:可
我正在玩 Rank-N-type 并尝试输入 x x .但我发现这两个函数可以以相同的方式输入,这很不直观。 f :: (forall a b. a -> b) -> c f x = x x g ::
这个问题已经有答案了: How do you compare two version Strings in Java? (31 个回答) 已关闭 8 年前。 有谁知道如何在Java中比较两个版本字符串
这个问题已经有答案了: How do the post increment (i++) and pre increment (++i) operators work in Java? (14 个回答)
下面是带有 -n 和 -r 选项的 netstat 命令的输出,其中目标字段显示压缩地址 (127.1/16)。我想知道 netstat 命令是否有任何方法或选项可以显示整个目标 IP (127.1.
我知道要证明 : (¬ ∀ x, p x) → (∃ x, ¬ p x) 证明是: theorem : (¬ ∀ x, p x) → (∃ x, ¬ p x) := begin intro n
x * x 如何通过将其存储在“auto 变量”中来更改?我认为它应该仍然是相同的,并且我的测试表明类型、大小和值显然都是相同的。 但即使 x * x == (xx = x * x) 也是错误的。什么
假设,我们这样表达: someIQueryable.Where(x => x.SomeBoolProperty) someIQueryable.Where(x => !x.SomeBoolProper
我有一个字符串 1234X5678 我使用这个正则表达式来匹配模式 .X|..X|X. 我得到了 34X 问题是为什么我没有得到 4X 或 X5? 为什么正则表达式选择执行第二种模式? 最佳答案 这里
我的一个 friend 在面试时遇到了这个问题 找到使该函数返回真值的 x 值 function f(x) { return (x++ !== x) && (x++ === x); } 面试官
这个问题在这里已经有了答案: 10年前关闭。 Possible Duplicate: Isn't it easier to work with foo when it is represented b
我是 android 的新手,我一直在练习开发一个针对 2.2 版本的应用程序,我需要帮助了解如何将我的应用程序扩展到其他版本,即 1.x、2.3.x、3 .x 和 4.x.x,以及一些针对屏幕分辨率
为什么案例 1 给我们 :error: TypeError: x is undefined on line... //case 1 var x; x.push(x); console.log(x);
代码优先: # CASE 01 def test1(x): x += x print x l = [100] test1(l) print l CASE01 输出: [100, 100
我正在努力温习我的大计算。如果我有将所有项目移至 'i' 2 个空格右侧的函数,我有一个如下所示的公式: (n -1) + (n - 2) + (n - 3) ... (n - n) 第一次迭代我必须
给定 IP 字符串(如 x.x.x.x/x),我如何或将如何计算 IP 的范围最常见的情况可能是 198.162.1.1/24但可以是任何东西,因为法律允许的任何东西。 我要带198.162.1.1/
在我作为初学者努力编写干净的 Javascript 代码时,我最近阅读了 this article当我偶然发现这一段时,关于 JavaScript 中的命名空间: The code at the ve
我正在编写一个脚本,我希望避免污染 DOM 的其余部分,它将是一个用于收集一些基本访问者分析数据的第 3 方脚本。 我通常使用以下内容创建一个伪“命名空间”: var x = x || {}; 我正在
我尝试运行我的test_container_services.py套件,但遇到了以下问题: docker.errors.APIError:500服务器错误:内部服务器错误(“ b'{” message
是否存在这两个 if 语句会产生不同结果的情况? if(x as X != null) { // Do something } if(x is X) { // Do something } 编
我是一名优秀的程序员,十分优秀!