Ich habe ein Problem mit einer Planetensimulation. Ich will den Flächeninhalt bestimmen, der von der verbindungslinie zwischen planet und Sonne überstrichen wird. wie man sich denken kann, kommen riesige Zahlen raus (die Erde ist 150 Mio. km von der Sonne entfernt). Ich berechne diese Fläche, indem ich sie in viele Dreiecke aufteile, die einzeln mit dieser Funktion berechnet werden:

long double Flaeche(long double xa, long double ya, long double xb, long double yb, long double xc, long double yc)
   {long double a = sqrtl(powl((xb-xc), 2.0) + powl((yb-yc), 2.0));
    long double b = sqrtl(powl((xc-xa), 2.0) + powl((yc-ya), 2.0));
    long double c = sqrtl(powl((xa-xb), 2.0) + powl((ya-yb), 2.0));
    long double s = (a+b+c)/2.0;
    long double Flaecheninhalt = (sqrtl((s*(s-a)*(s-b)*(s-c))));
    if (Flaecheninhalt<0){
       Flaecheninhalt = Flaecheninhalt * (-1.0);}
    return Flaecheninhalt;
    }

Die Teilflächen werden dann addiert.

Dabei kommen aber manchmal Zahlen raus, die kleiner als 0 sind. Die if-Abfrage

if (Flaecheninhalt<0){
       Flaecheninhalt = Flaecheninhalt * (-1.0);}

bemerkt dies aber nicht.

Läuft der long double über?

könnte dies negativ werden?

long double i = 9.21369875421325665478;
    char zwischenspeicher[1000];
    FILE *fl;
    fl = fopen("Flächeninhalte.txt", "w");
    sprintf(zwischenspeicher, "%g\n", i);
    fprintf(fl, zwischenspeicher);
    fclose(fl);

In der Datei stand jedoch -3.96111e-215 ...

Kann mir jemand helfen?

ofstream f1("Flächeninhalte.txt");
f1<<i<<endl;

(btw, du kannst die fprintf()-Ausgaben genauso formatieren wie bei sprintf() ;))

s^4 + (s^3)*a + (s^2)*b*c + s*a*b*c - (s^3)*b - (s^3)*c - (s^2)*a*b -(s^2)*a*c

was zeigt das bei dem Minusteil das C mit höher ptonzierten S multipliziert wird, daraus folgt ein sehr großen c bei sehr kleinem a und b wird s durch die formel s=(a+b+c)/2; kleiner c und das um einiges. Daher kommen evtl falsch eingelesene wete rein oder der gleichen, weil bei einem dreieck (a,b,c) sollte das normal nicht auftreten. jeddoch lieben wir alle das programmieren weil dinge passieren die so nicht passieren sollten oder^^? Grundprinzipiel kann die formel negative werte ausgeben, ausser ich habe falsch aufgelöst und normal sollte sie aber bei seitenlänmgen a,b,c eines dreiecks immer positiv sein.

long double i = 9.21369875421325665478;
ofstream f1("Flächeninhalte.txt"); 
f1<<i<<endl;

@oelsni
Nun verrat mir bitte mal wie bei einen Dreieck eine Seite sehr groß und die beiden anderen sehr klein sein sollen?

krümmung des raums aufgrund von masse-singularitäten

vielleicht sind die berechnungen ja allesamt richtig. nur unser verstand ist zu klein, um dies zu begreifen

#include <iostream>
#include <sstream>
#include <iomanip>
using namespace std;

int main()
{
    long double i = 9.21369875421325665478;
    stringstream f1;
    f1 << std::setprecision(15) << i << endl; //mehr als 15 Stellen sind auch bei long double nicht drin
    cout << f1.str() << endl;;
    return 0;
}