Почему численное решение совпадает с аналитическим решением на языке Си? - PullRequest
Новая
       8

Почему численное решение совпадает с аналитическим решением на языке Си?

0 голосов
/ 12 сентября 2018

Я кодировал задачу cfd 1 Dimmension, но мое численное решение приближается к аналитическому решению (до 6 знаков после запятой).

Я использую метод TDMA для численного решения и для аналитического решения. Я прямо подставляю значение x в функцию T (x).

Аналитическое решение T (x) оказываетсяT(x) = -(x^2)/2 +11/21(x);

E.г.Тогда 4 точки сетки;

  x0  = 0.000000,    x1  = 0.333333 ,   x2  = 0.666666 ,   x3  = 0.999999 .
T(x0) = 0.000000 , T(x1) = 0.119048 , T(x2) = 0.126984 , T(x3) = 0.023810.

А для численного решения я использовал технику TDMA, см. Код ниже.

Введите n = 4 для результатов. 

#include<stdio.h>

void temp_matrix(int n, double *a, double *b, double *c, double *d, double *T);

int main() {
  int Bi = 20.0;
  int n;
  printf("%s ", "Enter the Number of total Grid Points");
  scanf("%d", &n);
  float t = (n - 1);
  double dx = 1.0 / t;
  int i;
  printf("\n");

  double q; // analytical solution below
  double z[n];
  for (i = 0; i <= n - 1; i++) {
    q = (dx) * i;
    z[i] = -(q * q) / 2 + q * (11.0 / 21);
    printf("\nT analytical %lf ", z[i]);
  }

  double b[n - 1];
  b[n - 2] = -2.0 * Bi * dx - 2.0;
  for (i = 0; i <= n - 3; i++) {
    b[i] = -2.0;
  }

  double a[n - 1];
  a[n - 2] = 2.0;
  a[0] = 0;
  for (i = 1; i < n - 2; i++) {
    a[i] = 1.0;
  }

  double c[n - 1];
  for (i = 0; i <= n - 2; i++) {
    c[i] = 1.0;
  }

  double d[n - 1];
  for (i = 0; i <= n - 2; i++) {
    d[i] = -(dx * dx);
  }

  double T[n];
  temp_matrix(n, a, b, c, d, T);

  return 0;
}

void temp_matrix(int n, double *a, double *b, double *c, double *d, double *T) {
  int i;
  double beta[n - 1];
  double gama[n - 1];
  beta[0] = b[0];
  gama[0] = d[0] / beta[0];
  for (i = 1; i <= n - 2; i++) {
    beta[i] = b[i] - a[i] * (c[i - 1] / beta[i - 1]);
    gama[i] = (d[i] - a[i] * gama[i - 1]) / beta[i];
  }
  int loop;
  for (loop = 0; loop < n - 1; loop++)
    for (loop = 0; loop < n - 1; loop++)

      T[0] = 0;
  T[n - 1] = gama[n - 2];

  for (i = n - 2; i >= 1; i--) {
    T[i] = gama[i - 1] - (c[i - 1] * (T[i + 1])) / beta[i - 1];
  }
  printf("\n");
  for (i = 0; i < n; i++) {
    printf("\nT numerical %lf", T[i]);

  }
}

1 Ответ

0 голосов
/ 12 сентября 2018

Почему численное решение совпадает с аналитическим решением на языке Си?

Они отличаются примерно на 3 бита.

Печатайте с достаточной точностью, чтобы увидеть разницу.

Используя приведенное ниже, мы видим разницу в последнем шестнадцатеричном значении x620 против x619 T[3]. Это только 1 часть в разнице 10 15 . 

#include<float.h>
printf("T analytical %.*e\t%a\n", DBL_DECIMAL_DIG - 1, z[i], z[i]);
printf("T numerical  %.*e\t%a\n", DBL_DECIMAL_DIG - 1, T[i], T[i]);

C позволяет выполнять double math при long double math при FLT_EVAL_METHOD == 2, а затем те же аналитические / численные результаты. Ваши результаты могут отличаться от моих из-за этого, а также из-за других тонких нюансов FP. 

printf("FLT_EVAL_METHOD %d\n", FLT_EVAL_METHOD);

выход 

T analytical 0.0000000000000000e+00 0x0p+0
T analytical 1.1904761904761907e-01 0x1.e79e79e79e7ap-4
T analytical 1.2698412698412700e-01 0x1.0410410410411p-3
T analytical 2.3809523809523836e-02 0x1.861861861862p-6

T numerical  0.0000000000000000e+00 0x0p+0
T numerical  1.1904761904761904e-01 0x1.e79e79e79e79ep-4
T numerical  1.2698412698412698e-01 0x1.041041041041p-3
T numerical  2.3809523809523812e-02 0x1.8618618618619p-6

FLT_EVAL_METHOD 0
...