How did it come to be that stable evaluation of log(1 + x), called log1p in C99 and logp1 in IEEE-754, are apparently still not available in #Fortran as of #Fortran2023?