| 
									
										
										
										
											2015-10-04 05:28:11 +02:00
										 |  |  | ///////////////////////////////////////////////////////////////////////////////////
 | 
					
						
							|  |  |  | // Copyright (C) 2015 Edouard Griffiths, F4EXB                                   //
 | 
					
						
							|  |  |  | //                                                                               //
 | 
					
						
							|  |  |  | // This program is free software; you can redistribute it and/or modify          //
 | 
					
						
							|  |  |  | // it under the terms of the GNU General Public License as published by          //
 | 
					
						
							|  |  |  | // the Free Software Foundation as version 3 of the License, or                  //
 | 
					
						
							| 
									
										
										
										
											2019-04-11 14:32:15 +02:00
										 |  |  | // (at your option) any later version.                                           //
 | 
					
						
							| 
									
										
										
										
											2015-10-04 05:28:11 +02:00
										 |  |  | //                                                                               //
 | 
					
						
							|  |  |  | // This program is distributed in the hope that it will be useful,               //
 | 
					
						
							|  |  |  | // but WITHOUT ANY WARRANTY; without even the implied warranty of                //
 | 
					
						
							|  |  |  | // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the                  //
 | 
					
						
							|  |  |  | // GNU General Public License V3 for more details.                               //
 | 
					
						
							|  |  |  | //                                                                               //
 | 
					
						
							|  |  |  | // You should have received a copy of the GNU General Public License             //
 | 
					
						
							|  |  |  | // along with this program. If not, see <http://www.gnu.org/licenses/>.          //
 | 
					
						
							|  |  |  | ///////////////////////////////////////////////////////////////////////////////////
 | 
					
						
							|  |  |  | 
 | 
					
						
							|  |  |  | #include "util/db.h"
 | 
					
						
							|  |  |  | #include <cmath>
 | 
					
						
							|  |  |  | 
 | 
					
						
							| 
									
										
										
										
											2017-08-17 17:17:25 +02:00
										 |  |  | double CalcDb::dbPower(double magsq, double floor) | 
					
						
							| 
									
										
										
										
											2015-10-04 05:28:11 +02:00
										 |  |  | { | 
					
						
							| 
									
										
										
										
											2020-11-04 20:05:11 +01:00
										 |  |  |     if (floor <= 0.0) { | 
					
						
							|  |  |  | 		return -100.0; | 
					
						
							|  |  |  | 	} | 
					
						
							| 
									
										
										
										
											2017-08-17 17:17:25 +02:00
										 |  |  | 
 | 
					
						
							| 
									
										
										
										
											2020-11-04 20:05:11 +01:00
										 |  |  | 	if (magsq > floor) { | 
					
						
							| 
									
										
										
										
											2017-08-17 17:17:25 +02:00
										 |  |  | 	    return 10.0 * log10(magsq); | 
					
						
							| 
									
										
										
										
											2020-11-04 20:05:11 +01:00
										 |  |  | 	} else { | 
					
						
							| 
									
										
										
										
											2017-08-17 17:17:25 +02:00
										 |  |  | 		return 10.0 * log10(floor); | 
					
						
							| 
									
										
										
										
											2015-10-04 05:28:11 +02:00
										 |  |  | 	} | 
					
						
							|  |  |  | } | 
					
						
							| 
									
										
										
										
											2017-07-25 21:21:48 +02:00
										 |  |  | 
 | 
					
						
							|  |  |  | double CalcDb::powerFromdB(double powerdB) | 
					
						
							|  |  |  | { | 
					
						
							|  |  |  |     return pow(10.0, powerdB / 10.0); | 
					
						
							|  |  |  | } | 
					
						
							| 
									
										
										
										
											2021-06-27 01:52:42 +02:00
										 |  |  | 
 | 
					
						
							|  |  |  | double CalcDb::frexp10(double arg, int *exp) | 
					
						
							|  |  |  | { | 
					
						
							|  |  |  |    *exp = (arg == 0) ? 0 : 1 + (int)std::floor(std::log10(std::fabs(arg) ) ); | 
					
						
							|  |  |  |    return arg * std::pow(10 , -(*exp)); | 
					
						
							|  |  |  | } |