windows - Java高精度運(yùn)算問題求助
問題描述
公司項(xiàng)目里面需要做大量的高精度運(yùn)算,剛開始用double類型運(yùn)算,后來發(fā)現(xiàn)用double類型運(yùn)算后有些值的精確度超過了理想范圍,就用了BigDecimal來計(jì)算,現(xiàn)在的問題是BigDecimal的運(yùn)算效率比double慢幾十倍,數(shù)據(jù)量大的話,慢的要死。請(qǐng)問有沒有好的解決方案?這個(gè)問題急需解決。
//相關(guān)性系數(shù)計(jì)算public BigDecimal getRelativityTool_bydim(RelativityTool u) {BigDecimal sim = new BigDecimal('0'); //最后的皮爾遜相關(guān)度系數(shù)BigDecimal common_items_len = new BigDecimal(this.rating_map_list.size()); //操作數(shù)的個(gè)數(shù)BigDecimal this_sum = new BigDecimal('0'); //第一個(gè)相關(guān)數(shù)的和BigDecimal u_sum = new BigDecimal('0'); //第二個(gè)相關(guān)數(shù)的和BigDecimal this_sum_sq = new BigDecimal('0'); //第一個(gè)相關(guān)數(shù)的平方和BigDecimal u_sum_sq = new BigDecimal('0'); //第二個(gè)相關(guān)數(shù)的平方和BigDecimal p_sum = new BigDecimal('0'); //兩個(gè)相關(guān)數(shù)乘積的和for (int i = 0; i < this.rating_map_list.size(); i++) { BigDecimal this_grade = this.rating_map_list.get(i); BigDecimal u_grade = u.rating_map_list.get(i); //評(píng)分求和 //平方和 //乘積和 this_sum = this_sum.add(this_grade); u_sum = u_sum.add(u_grade); this_sum_sq = this_sum_sq.add(this_grade.pow(2)); u_sum_sq = u_sum_sq.add(u_grade.pow(2)); p_sum = p_sum.add(this_grade.multiply(u_grade));}BigDecimal num = common_items_len.multiply(p_sum).subtract(this_sum.multiply(u_sum));BigDecimal den = sqrt(common_items_len.multiply(this_sum_sq).subtract(this_sum.pow(2)).multiply(common_items_len.multiply(u_sum_sq).subtract(u_sum.pow(2))));if (den.compareTo(new BigDecimal('0')) == 0) { sim = new BigDecimal('1');} else { sim = num.pide(den,5, BigDecimal.ROUND_HALF_UP);}return sim; } //大數(shù)字開方 public static BigDecimal sqrt(BigDecimal x) {BigDecimal n1 = BigDecimal.ONE;BigDecimal ans = BigDecimal.ZERO;while ((n1.multiply(n1).subtract(x)).abs().compareTo(BigDecimal.valueOf(0.001)) == 1) { BigDecimal s1 = x.pide(n1, 2000, BigDecimal.ROUND_HALF_UP); BigDecimal s2 = n1.add(s1); n1 = s2.pide(BigDecimal.valueOf(2), 2000, BigDecimal.ROUND_HALF_UP);}ans = n1;BigDecimal rt = new BigDecimal(ans.toString().split('.')[0]);return rt; }
問題解答
回答1:除了使用C或者C++來做高精度運(yùn)算之外,好像沒有什么辦法可以同時(shí)兼顧性能和精度了。
回答2:大學(xué)計(jì)算機(jī)專業(yè)有門課程叫“計(jì)算方法”,專門探討如何在精度有限的計(jì)算過程中保持誤差最小化。樓主有興趣的話可以找下相關(guān)教材。
回答3:后來發(fā)現(xiàn)用double類型運(yùn)算后有些值的精確度超過了理想范圍是超過還是滿足不了?
這里有一段計(jì)算平方根的代碼,我從stackoverflow上找到的,在我自己的機(jī)子測試要比你上面那個(gè)快十倍左右。所以: 一則你可以通過改進(jìn)的算法來提高性能,其二,最好的辦法找一些已有的library來直接用:例如這個(gè)上面列的
public void test_sqrt() { BigDecimal x = BigDecimal.valueOf(Long.MAX_VALUE); BigDecimal x0 = BigDecimal.ZERO; BigDecimal x2 = new BigDecimal(2); BigDecimal x1 = new BigDecimal(Math.sqrt(x.doubleValue())); while (!x0.equals(x1)) {x0 = x1;x1 = x.pide(x0, 2000, BigDecimal.ROUND_HALF_UP);x1 = x1.add(x0);x1 = x1.pide(x2, 2000, BigDecimal.ROUND_HALF_UP); } assertEquals(3037000499L, x1.longValue());}
相關(guān)文章:
1. 如何解決Centos下Docker服務(wù)啟動(dòng)無響應(yīng),且輸入docker命令無響應(yīng)?2. 我在centos容器里安裝docker,也就是在容器里安裝容器,報(bào)錯(cuò)了?3. javascript - 微信h5發(fā)送圖文信息,部分設(shè)備點(diǎn)擊“發(fā)送”按鈕時(shí)沒反應(yīng),問題較難重現(xiàn),如何能找到可能存在問題的點(diǎn)?4. javascript - 音樂播放器-圖片旋轉(zhuǎn)5. android - 為 AppBarLayout 設(shè)置的背景圖片 TransitionDrawable 為什么只在第一次打開的時(shí)候有效?6. javascript - js中遞歸與for循環(huán)同時(shí)發(fā)生的時(shí)候,代碼的執(zhí)行順序是怎樣的?7. MySQL timestamp的默認(rèn)值怎么設(shè)置?8. android - 使用百度sdk調(diào)用SDKInitializer.initialize(this)時(shí)報(bào)錯(cuò)?9. objective-c - IOS 分享到微信 提示 應(yīng)用消息數(shù)據(jù)錯(cuò)誤10. docker 17.03 怎么配置 registry mirror ?
