הרבה זמן לא נגעתי בתלת מימד, עכשיו חזרתי לזה קצת ואני רוצה מחדש (באופן משופר) להכין סיבובים תלת מימדיים של ווקטורים בעזרת מטריצות.

בכל מקרה, כמה שעות עכשיו חירבשתי את כל הקוד שלי בשביל למצוא מה לעזאזל גורם לו לא לפעול, וסוף סוף, מצאתי את הבעיה.
נראה שcos מחזיר ערכים בלתי הגיוניים.
cos על 90 מעלות (ברדיאנים כמובן) מחזיר לי 4.37114e-008-, ןגם הצלחתי להחזיר 8 נקודה משהו.
איך בדיוק זה הגיוני? הקוד הזה פעל לי פעם. אין לי מושג מה הסיפור של C++, אבל זה לא הבאג (או בעיה) היחיד שמצאתי לאחרונה עם הספרייה הבסיסית...

קוד PHP:

 #include <iostream>
using namespace std;

#include <math.h>

#define DEG2RAD 0.01745329251994329576923690768488
#define RAD2DEG 57.2957795130823208767981548141052

#define Deg2Rad(a)((float)(a * DEG2RAD))
#define Rad2Deg(a)((float)(a * RAD2DEG))

int main(int argc, char *argv[])
{
    
    float a = Deg2Rad(90.0f);

    cout<<cos(a)<<"\n";
    cout<<sin(a)<<"\n";
    
    system("pause");
    return EXIT_SUCCESS;
} 


א. אתה יכול לעבוד ב double. ברגע שהעדפת float איבדת מידה של דיוק.
ב. כידוע גם double אינו מדויק.
ג. פונקציות sin ו cos הן "שטוחות" (נגזרת קרובה לאפס) במקסימום/מינימום שלהן ו"תלולות" באזור האפס. לכן טעות החישוב (בשל אי הדיוק של נקודה צפה) מורגשת יותר סביב אפס. שגיאה של E-16 (אם תעבור ל double) היא קטנה למדי.

לא דרוש לי דיוק מירבי, אבל זה בכל מקרה לא מסביר איך פונקצייה טריגונומטרית מחזירה ערכים שיותר גדולים מ 1.
שיניתי את הקוד לdouble סתם לראות איזה באג הפעם יהיה, זה הפך את המספר ל 6 נקודה משהו במקום 4 נקודה משהו.
מה הסיפור של זה?

איפה ראית ערך גדול מ 1?
אני מניח שאתה מודע לזה ש 4.37114e-008- הוא למעשה 0.0000000437114-.

הראל

לא, לא הייתי מודע לזה.
עדיין לא ברור לי איך זה לא שווה ל0, אבל נו שיהיה.

תודה.

תשתמש בספריות מתמטיקה ייעודיות כגון GMP שמציעות דיוק בלתי מוגבל במהירות גבוהה
http://gmplib.org/

ציטוט:

GMP is a free library for arbitrary precision arithmetic, operating on signed integers, rational numbers, and floating point numbers. There is no practical limit to the precision except the ones implied by the available memory in the machine GMP runs on. GMP has a rich set of functions, and the functions have a regular interface. The main target applications for GMP are cryptography applications and research, Internet security applications, algebra systems, computational algebra research, etc. GMP is carefully designed to be as fast as possible, both for small operands and for huge operands. The speed is achieved by using fullwords as the basic arithmetic type, by using fast algorithms, with highly optimised assembly code for the most common inner loops for a lot of CPUs, and by a general emphasis on speed.

_________________________________________________

אזהרה: משרד הבריאות קובע כי העישון מזיק לבריאות !
תראו, אפילו החייזר נהיה ירוק מזה

מנהל פורום מדע בדיוני ופנטסיה ופורום מדע טכנולוגיה וטבע בפרש