你需要一个任意长度的字符串作为答案吗?或者您需要任意精度或任意指数大小?
或者……双精度浮点答案(返回值)就足够了;鉴于我们“只”使用了一些“任意”大小的对数?
双精度浮点数具有11位有符号指数:因此,如果您的大数字字符串的长度≤1022位≈307十进制数字(所以字符串长度为306个字符,包括小数点),您是安全的!更准确地说,如果得到的十进制指数的绝对值≤307,则应该是安全的.你需要更大的指数吗? (换句话说,我想:你是在使用现实世界的数字还是理论/纯数学?)
为什么不使用一些字符串处理,以及一些简单的浮点日志算法?对于任何真实世界的数字,这应该非常快……
function bclog10($n){
//←Might need to implement some validation logic here!
$pos=strpos($n,'.');
if($pos===false){
$dec_frac='.'.substr($n,0,15);$pos=strlen($n);
}else{ $dec_frac='.'.substr(substr($n,0,$pos).substr($n,$pos+1),0,15);
}
return log10((float)$dec_frac)+(float)$pos;
}
您可以使用一些众所周知的对数算法转换基数:
function bclogn($n,$base=M_E){//$base should be float: default is e
return bclog10($n)*log(10)/log($base);
}
我已经测试了这些功能,它们适合我,我提供的例子;提供与Windows 10计算器完全相同的答案,达到PHP使用的双精度算术的极限.
如果你实际上需要超过15位精度和超过307的十进制指数,你可能能够实现自己的“BigFloat”类对象,并以某种方式使用标准的内置浮点函数构建其方法分而治之的方法!那么也许,我们可以将其作为任意精度浮点对数算法的基础,通过将其与上述函数/技术相结合.您可能需要考虑在math.stackexchange.com与人们协商,以了解更多关于这是否可行的方法.
主要编辑:第二次尝试……
function bclog10($n){//By Matthew Slyman @aaabit.com
$m=array();// ↓ Validation, matching/processing regex…
preg_match('/^(-)?0*([1-9][0-9]*)?(\.(0*))?([1-9][0-9]*)?([Ee](-)?0*([1-9][0-9]*))?$/',$n,$m);
if(!isset($m[1])){throw new \Exception('Argument: not decimal number string!');}
$sgn=$m[1];if($sgn==='-'){throw new \Exception('Cannot compute: log(
$abs=$m[2];$pos=strlen($abs);
if(isset($m[4])){$fre=$m[4];}else{$fre='';}$neg=strlen($fre);
if(isset($m[5])){$frc=$m[5];}else{$frc='';}
if(isset($m[7])){$esgn=$m[7]==='-'?-1:1;}else{$esgn=1;}
if(isset($m[8])){$eexp=$m[8];}else{$eexp=0;}
if($pos===0){
$dec_frac='.'.substr($frc,0,15);$pos=-1*$neg;
}else{ $dec_frac='.'.substr($abs.$fre.$frc,0,15);
}
return log10((float)$dec_frac)+(float)$pos+($esgn*$eexp);
}