- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我有一个 const int 类型的变量,但它所依赖的参数是 double 类型的。当我尝试将其从“double”转换为“const int”时,它无法正常工作。比如N应该是991,就输入990。我试过好几种方法,只有一种有效,但我不确定这种方法是否一直有效。以下是我尝试过的一些方法:
第一种方法:
const int N = (Ls-1)/dx + 1;
第二种方法:
const int N = static_cast<const int>((Ls-1)/dx) + 1;
第三种方法:
double Z = (Ls-1)/dx + 1;
const int N = Z;
第四种方法(唯一可行的方法):
double Z = (Ls-1)/dx;
const int N = Z + 1;
请注意 dx 是一个值,这样 (Ls-1)/dx 的余数将始终为零(即它始终是一个整数值)。无论如何可以解释为什么其他方法不起作用,以便我可以更好地理解类型转换?
编辑:根据要求,我正在上传整个代码以展示一切是如何工作的:
#include <iostream>
#include <math.h>
#include <stdio.h>
#include <fstream>
#include <cmath>
#include <algorithm>
#define pi 3.14159265
using namespace std;
//Define Fluid Properties
double rho_L = 998; //Liquid Density
double rho_LG = 828.9; //Liquid-Gas Density Ratio
double mu_L = 0.000798; //Liquid Viscosity
double mu_LG = 40.24; //Liquid-Gas Viscosity Ratio
double sigma = 0.0712; //Surface Tension
double nu_G = (mu_L/mu_LG)/(rho_L/rho_LG);
//Define Injector Properties
double Uinj = 56.7; //Injection Velocity
double Dinj = 0.0998; //Injector Diameter
double theta = 15.0*pi/180.0; //Spray Cone Angle
double L = 500.0*Dinj; //Atomization Length
double Ls = L/Dinj; //Normalized Atomization Length
//Define Solver Parameters
double K = 5294; //Viscous Dissipation Coefficient
double Eps = pow(10,-5); //Residual Error
double dx = 0.0001; //Step Size
double Ui = 10; //Initial Guess
//const int Z = static_cast<const int>((Ls-1)/dx + 1) + 1;
const int N = (Ls-1)/dx + 1;//Z;
double deriv (double U, double X, double delta, double m)
{
double dudx;
dudx = -(1.0/delta)*(1.0/U)*(U - sqrt(1.0 - U)/sqrt(m*X*X))*(U - sqrt(1.0 - U)/sqrt(m*X*X));
return (dudx);
}
int main()
{
//Declare Variables
int max_step;
double ERR;
int step;
double DEN;
double SMD;
double m;
double Ug;
double Re;
double Cd;
double delta;
double K1;
double K2;
double K3;
double K4;
//Allocate Memory From Heap
double *U = new double [N];
double *X = new double [N];
//Initialize Vectors and Variables
DEN = 0.5*rho_L - (4.0/3.0)*K*(mu_L)/(Uinj*Dinj*Dinj)*L;
m = 4.0/rho_LG*tan(theta)*tan(theta);
for (int i = 0; i < N; i++)
{
X[i] = 1.0 + dx*i;
}
U[0] = 1.0;
max_step = 1;
ERR = 1;
step = 0;
while(abs(ERR) > Eps && step < max_step)
{
//Calculate Ug
Ug = sqrt(1.0 - (Ui/Uinj))/sqrt(m*Ls*Ls)*Uinj;
//Calculate SMD
SMD = 6.0*sigma/(DEN*(Uinj*Uinj - Ui*Ui));
//Calculate Re # and Drag Coefficient
Re = abs(Ui-Ug)*SMD/nu_G;
if(Re <= 0.01)
{
Cd = (0.1875) + (24.0/Re);
}
else if(Re > 0.01 && Re <= 260.0)
{
Cd = (24.0/Re)*(1.0 + 0.1315*pow(Re,(0.32 - 0.05*log10(Re))));
}
else
{
Cd = (24.0/Re)*(1.0 + 0.1935*pow(Re,0.6305));
}
//Determine New U
delta = (4.0/3.0)*(1.0/Cd)*(rho_LG)*(SMD/Dinj);
//RK4
for (int i = 0; i < N-1; i++)
{
K1 = deriv(U[i],X[i],delta,m);
K2 = deriv(U[i]+0.5*dx*K1,X[i]+0.5*dx,delta,m);
K3 = deriv(U[i]+0.5*dx*K2,X[i]+0.5*dx,delta,m);
K4 = deriv(U[i]+dx*K3,X[i+1],delta,m);
U[i+1] = U[i] + dx/6.0*(K1 + 2.0*K2 + 2.0*K3 + K4);
//if(i >= 0 && i <= 3)
//cout << i << " " << K1 << " " << K2 << " " << K3 << " " << K4 << " " << U[i] << endl;
}
ERR = abs(U[N-1]*Uinj - Ui)/Ui;
Ui = U[N-1]*Uinj;
step = step + 1;
}
SMD = 6.0*sigma/(DEN*(Uinj*Uinj - Ui*Ui));
cout << "U = " << Ui << endl;
cout << "SMD = " << SMD << endl;
cout << "DEN = " << DEN << endl;
cout << "Ug = " << Ug << endl;
cout << "m = " << m << endl;
cout << "delta = " << delta << endl;
cout << "Re = " << Re << endl;
cout << "Cd = " << Cd << endl;
cout << "U* = " << U[N-1] << endl;
cout << "Error = " << ERR << endl;
cout << "step = " << step << endl;
//Output Data Into Text File
ofstream outputdata("result-500-15.txt");
for (int i = 0; i < N; i++)
{
outputdata << X[i] << " " << U[i] << '\n';
}
outputdata.close();
delete [] U;
delete [] X;
return 0;
}
最佳答案
您的猜测是正确的:0.1 没有有限的二进制表达式。这是一个相当复杂的问题,并且有许多极端情况,一般不会通过添加评论中提到的 0.01 来解决。 (这在很大程度上取决于您期望的值等)
你的问题表明商总是应该是一个整数。在这种情况下,保持正确结果的正确方法是不使用任何 double
开头(对于 Ls
、dx
、Z
).要么使用小数类型(C++ 中没有内置,使用你自己的或库),要么使用任意精度的小数类型(同样,使用像 gmp
这样的库 - 如果你知道所有的数字,这是明智的有一个终止的十进制表达式),或者,最简单的:如果 Ls
和 dx
都保证小数点后最多有 n
位,两者都乘以 10^n
并使用整数类型。
好的,您的代码与我的预期大不相同。在这种情况下,我认为正确的做法是固定步数 N
并据此计算 dx
,而不是反过来计算:
const int N = 10000;
double dx = (Ls-1.0)/(double)(N-1);
如果您想从 dx 的值开始并选择 N 以便计算出 dx 的值,请在程序启动时询问用户:
#include <cmath>
double dxestim;
cout << "dx should be close to: ";
cin >> dxestim;
cout << "Candidate values for N: " << endl;
int N1 = (int) floor((Ls-1)/dx + 1.0);
int N2 = (int) ceil((Ls-1)/dx + 1.0);
cout << N1 << " gives dx = " << (Ls-1.0)/(double)(N1-1) << endl;
cout << N2 << " gives dx = " << (Ls-1.0)/(double)(N2-1) << endl;
cout << "Please choose N: ";
cin >> N;
...
关于从 double 到 const int 的 C++ 类型转换不能正常工作,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/14783234/
自从我 faced an issue由于背景图片对于不同分辨率的内容来说太短,我尝试将背景分成 3 部分并自动拉伸(stretch)中间部分以相应地填充顶部和底部图像之间的空间。不幸的是我没能在 CS
我从去年开始就在我的程序中运行这个函数(Linux 和 Windows)。 现在我需要实现一个新功能,我的新构建不再运行。 我还有其他使用 POST 的 CUrl 函数,结果是一样的:没问题,但我的
在评估函数应用方面,Haskell 是只支持普通降阶还是也支持应用降阶?我是否认为正常顺序是 Haskell 惰性的原因? 最佳答案 GHC 运行时不使用术语缩减策略,因为那会非常低效。事实上,GHC
怎么来的multi使用多处理池对多个“进程”上的数据进行分段和处理的函数比仅调用 map 慢(8 秒)。功能(6 秒)? from multiprocessing import Pool import
假设我正在渲染一个 3d GL_TRIANGLE。该对象需要 3 个顶点才能定义:A、B、C。我将此类数据放入缓冲区并通过 glVertexAttribPointer 将其绑定(bind)到着色器。
我有一个字体的三个文件,普通的,粗体的和浅色的。由于 font-weight:light 不存在,我该如何在 font-face 上设置 light 呢? 顺便问一下,font-weight:ligh
我是 C 的新手,我似乎无法弄清楚什么似乎是一个非常简单的指针问题。我的程序将行号添加到文件中。它逐行读入文件,然后在每行的开头添加一个行号。它在每个文件上都可以正常工作,如下所示: soccer@s
我有以下代码,我不确定为什么当它命中 Myclass 的析构函数时我会收到堆损坏检测错误。我相信我正在正确地释放内存?? #include #include using namespace std
有什么方法可以将“正常”数学符号解释为逆波兰符号 (RPN)..? 例如1) 2 + 3*4 - 1 = 234*+1-2) 5 (4-8) = 548- 你可以假设遵循 BODMAS 规则并且必须首
http://www.ergotopia.de/ergonomie-shop/ergonomische-kissen/orthopaedisches-sitzkissen的手机页面应该看起来像右边(检
我正在 Phonegap/Cordova 中构建一个应用程序。应用目前相当简单,但确实需要网络状态和地理定位插件才能工作。 到目前为止,我已经在 Android 上开发了该应用程序(目前它仅由一些基本
我一整天都在做这个,但没有运气 我设法在一行 TfidfVectorizer 中消除了问题 这是我的工作代码 from sklearn.feature_extraction.text import C
也许有人看到一个错误,问题是当我按btn2 (button 2)和btn3 (button 3)应用程序crashes时,但操作仍然有效,即video正在运行并且PDF打开,而button 1正常工作
我正在开发一个应用程序。它的第一页是登录屏幕。成功登录后,我想将用户带到选项卡式 Activity 。我怎样才能在安卓中做到这一点?谢谢 最佳答案 在 Android 中,启动 Activity 是通
我不确定我在这里做错了什么。 :normal! I### 当我对一个单词执行此命令时,我想要的最终结果是: ### word 但是我得到了这个: ###word 最佳答案 Vim 的 :normal是
我必须将 2 个静态矩阵发送到分配动态矩阵的函数,将矩阵 1 乘以矩阵 2,并返回新矩阵的地址。请注意,COMM 很常见。 我尝试删除 free_matrix 行,它工作正常。 void main()
我在我的一个项目中使用 Gnome libglib 并遇到了一个奇怪的错误。我可以输入 GList 的元素数量看起来仅限于 45 个。在第 45 个元素处,它给出了此错误 40 counter 41
我正在尝试获取“顶级”HWND 的尺寸。即,我想要 Firefox/Windows 资源管理器等的主 HWND 的当前尺寸。窗口。如果窗口最小化, GetWindowRect() 将不起作用。 Get
相同的标题:什么是索引 - 正常 - 全文 - 唯一? 最佳答案 普通索引用于通过仅包含行数据的切片或散列来加速操作。 全文索引向数据库的全文搜索 (FTS) 引擎指示它应该将数据存档在给定字段中,以
我正在使用 EnumParser来自 here它在 VC++ 中编译得很好,但是使用 gcc 我有这样的错误: ./Terminator.o: In function `EnumParser::Enu
我是一名优秀的程序员,十分优秀!