มีวิธีทางโปรแกรมในการรับ double ที่ใกล้เคียงกับ 1.0 มากที่สุด แต่จริงๆแล้วไม่ใช่ 1.0 หรือไม่?
วิธีแฮ็ควิธีหนึ่งในการทำเช่นนี้คือการ memcpy สองเท่าให้เป็นจำนวนเต็มขนาดเท่ากันแล้วลบออก วิธีการทำงานของรูปแบบจุดลอยตัวของ IEEE754 สิ่งนี้จะจบลงด้วยการลดเลขชี้กำลังทีละหนึ่งในขณะที่เปลี่ยนส่วนเศษส่วนจากศูนย์ทั้งหมด (1.000000000000) เป็นค่าทั้งหมด (1.111111111111) อย่างไรก็ตามมีเครื่องจักรที่เก็บจำนวนเต็ม little-endian ในขณะที่ floating-point ถูกจัดเก็บ big-endian ดังนั้นจึงไม่สามารถใช้ได้เสมอไป
nextafter()เป็นวิธีเดียวที่เหมาะสมในการบรรลุสิ่งที่เขาต้องการ
1.0000...ไบนารีกำลังลดลง0.111111....และเพื่อทำให้เป็นปกติคุณต้องเลื่อนไปทางซ้าย: 1.11111...ซึ่งคุณต้องลดเลขชี้กำลัง แล้วคุณอยู่ห่างจาก 1.0 2 ulp ดังนั้นไม่การลบค่าหนึ่งออกจากค่าอินทิกรัลไม่ได้ให้สิ่งที่ถามที่นี่
คำตอบ:
ใน C และ C ++ ค่าต่อไปนี้ให้ค่าที่ใกล้เคียงที่สุดกับ 1.0:
#include <limits.h>
double closest_to_1 = 1.0 - DBL_EPSILON/FLT_RADIX;
อย่างไรก็ตามโปรดทราบว่าใน C ++ เวอร์ชันที่ใหม่กว่าlimits.hจะถูกเลิกใช้งานในรูปแบบclimits. แต่ถ้าคุณใช้รหัสเฉพาะ C ++ อยู่แล้วคุณสามารถใช้
#include <limits>
typedef std::numeric_limits<double> lim_dbl;
double closest_to_1 = 1.0 - lim_dbl::epsilon()/lim_dbl::radix;
และตามที่ Jarod42 เขียนในคำตอบของเขาตั้งแต่ C99 หรือ C ++ 11 คุณยังสามารถใช้nextafter:
#include <math.h>
double closest_to_1 = nextafter(1.0, 0.0);
แน่นอนใน C ++ คุณสามารถ (และสำหรับ C ++ รุ่นที่ใหม่กว่า) ควรรวมcmathและใช้std::nextafterแทน
std::ceil(std::nextafter(1., std::numeric_limits<double>::max()))นอกจากนี้ยังเป็นวิธีที่ดีเพื่อเพิ่มคู่เป็นจำนวนเต็มซึ่งแสดงต่อไปนี้:
nextafter , และนี่คือวิธีการดำเนินการคิดถึงมันในกรณีที่คนอื่นต้องการที่จะเห็นว่ามันทำ โดยทั่วไป: บิดบิตดิบ
ใน C คุณสามารถใช้สิ่งนี้:
#include <float.h>
...
double value = 1.0+DBL_EPSILON;
DBL_EPSILON คือความแตกต่างระหว่าง 1 และค่าต่ำสุดที่มากกว่า 1 ที่สามารถแสดงได้
คุณจะต้องพิมพ์เป็นตัวเลขหลายหลักจึงจะเห็นค่าที่แท้จริง
บนแพลตฟอร์มของฉันให้printf("%.16lf",1.0+DBL_EPSILON)1.0000000000000002
1.0ถามเฉพาะเกี่ยวกับ BTW ยังให้ค่าที่ใกล้เคียงที่สุดมากกว่า 1 และไม่ใช่ค่าที่ใกล้เคียงที่สุดกับ 1 (ซึ่งอาจน้อยกว่า 1) ดังนั้นฉันยอมรับว่านี่เป็นคำตอบบางส่วน แต่ฉันคิดว่ามันสามารถมีส่วนร่วมได้
ใน C ++ คุณสามารถใช้สิ่งนี้ได้
1 + std::numeric_limits<double>::epsilon()