Compare commits
2 Commits
53cd72a4b0
...
c971550535
| Author | SHA1 | Date | |
|---|---|---|---|
| c971550535 | |||
| a42e5a7756 |
@@ -117,4 +117,33 @@ inline PIVector2D<double> abs(const PIVector2D<complexd> & v) {
|
|||||||
return result;
|
return result;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
||||||
|
/**
|
||||||
|
* @brief Inline funtion of compare with zero different types
|
||||||
|
*
|
||||||
|
* @param v is input parameter of type T
|
||||||
|
* @return true if zero, false if not zero
|
||||||
|
*/
|
||||||
|
template<typename T, typename std::enable_if<std::is_floating_point<T>::value, int>::type = 0>
|
||||||
|
inline bool PIMathFloatNullCompare(const T v) {
|
||||||
|
static_assert(std::is_floating_point<T>::value, "Type must be floating point");
|
||||||
|
return (piAbs(v) < T(1E-200));
|
||||||
|
}
|
||||||
|
|
||||||
|
/**
|
||||||
|
* @brief Inline funtion of compare with zero colmplexf type
|
||||||
|
*
|
||||||
|
* @param v is input parameter of type colmplexf
|
||||||
|
* @return true if zero, false if not zero
|
||||||
|
*/
|
||||||
|
template<typename T,typename std::enable_if<
|
||||||
|
std::is_floating_point<decltype(T::real)>::value ||
|
||||||
|
std::is_floating_point<decltype(T::imag)>::value
|
||||||
|
, int>::type = 0>
|
||||||
|
inline bool PIMathFloatNullCompare(const T v) {
|
||||||
|
static_assert(std::is_floating_point<decltype(v.real)>::value, "Type must be floating point");
|
||||||
|
static_assert(std::is_floating_point<decltype(v.imag)>::value, "Type must be floating point");
|
||||||
|
return (abs(v) < float(1E-200));
|
||||||
|
}
|
||||||
|
|
||||||
#endif // PIMATHCOMPLEX_H
|
#endif // PIMATHCOMPLEX_H
|
||||||
|
|||||||
@@ -28,40 +28,6 @@
|
|||||||
#include "pimathvector.h"
|
#include "pimathvector.h"
|
||||||
#include "pimathcomplex.h"
|
#include "pimathcomplex.h"
|
||||||
|
|
||||||
/**
|
|
||||||
* @brief Inline funtion of compare with zero different types
|
|
||||||
*
|
|
||||||
* @param v is input parameter of type T
|
|
||||||
* @return true if zero, false if not zero
|
|
||||||
*/
|
|
||||||
template<typename T>
|
|
||||||
inline bool _PIMathMatrixNullCompare(const T v) {
|
|
||||||
static_assert(std::is_floating_point<T>::value, "Type must be floating point");
|
|
||||||
return (piAbs(v) < T(1E-200));
|
|
||||||
}
|
|
||||||
|
|
||||||
/**
|
|
||||||
* @brief Inline funtion of compare with zero colmplexf type
|
|
||||||
*
|
|
||||||
* @param v is input parameter of type colmplexf
|
|
||||||
* @return true if zero, false if not zero
|
|
||||||
*/
|
|
||||||
template<>
|
|
||||||
inline bool _PIMathMatrixNullCompare<complexf>(const complexf v) {
|
|
||||||
return (abs(v) < float(1E-200));
|
|
||||||
}
|
|
||||||
|
|
||||||
/**
|
|
||||||
* @brief Inline funtion of compare with zero complexd type
|
|
||||||
*
|
|
||||||
* @param v is input parameter of type colmplexd
|
|
||||||
* @return true if zero, false if not zero
|
|
||||||
*/
|
|
||||||
template<>
|
|
||||||
inline bool _PIMathMatrixNullCompare<complexd>(const complexd v) {
|
|
||||||
return (abs(v) < double(1E-200));
|
|
||||||
}
|
|
||||||
|
|
||||||
|
|
||||||
/// Matrix templated
|
/// Matrix templated
|
||||||
|
|
||||||
@@ -1225,7 +1191,7 @@ public:
|
|||||||
for (uint k = i; k < _V2D::cols_; ++k) smat.element(k, j) -= mul * smat.element(k, i);
|
for (uint k = i; k < _V2D::cols_; ++k) smat.element(k, j) -= mul * smat.element(k, i);
|
||||||
}
|
}
|
||||||
if (i < _V2D::cols_ - 1) {
|
if (i < _V2D::cols_ - 1) {
|
||||||
if (_PIMathMatrixNullCompare(smat.element(i + 1, i + 1))) {
|
if (PIMathFloatNullCompare(smat.element(i + 1, i + 1))) {
|
||||||
if (ok != 0) *ok = false;
|
if (ok != 0) *ok = false;
|
||||||
return *this;
|
return *this;
|
||||||
}
|
}
|
||||||
@@ -1275,7 +1241,7 @@ public:
|
|||||||
if (sv != 0) (*sv)[j] -= mul * (*sv)[i];
|
if (sv != 0) (*sv)[j] -= mul * (*sv)[i];
|
||||||
}
|
}
|
||||||
if (i < _V2D::cols_ - 1) {
|
if (i < _V2D::cols_ - 1) {
|
||||||
if (_PIMathMatrixNullCompare(smat.element(i + 1, i + 1))) {
|
if (PIMathFloatNullCompare(smat.element(i + 1, i + 1))) {
|
||||||
if (ok != 0) *ok = false;
|
if (ok != 0) *ok = false;
|
||||||
return *this;
|
return *this;
|
||||||
}
|
}
|
||||||
|
|||||||
@@ -105,7 +105,7 @@ void PIMutex::init() {
|
|||||||
pthread_mutexattr_t attr;
|
pthread_mutexattr_t attr;
|
||||||
memset(&attr, 0, sizeof(attr));
|
memset(&attr, 0, sizeof(attr));
|
||||||
pthread_mutexattr_init(&attr);
|
pthread_mutexattr_init(&attr);
|
||||||
pthread_mutexattr_settype(&attr, PTHREAD_MUTEX_NORMAL);
|
pthread_mutexattr_settype(&attr, PTHREAD_MUTEX_RECURSIVE);
|
||||||
memset(&(PRIVATE->mutex), 0, sizeof(PRIVATE->mutex));
|
memset(&(PRIVATE->mutex), 0, sizeof(PRIVATE->mutex));
|
||||||
pthread_mutex_init(&(PRIVATE->mutex), &attr);
|
pthread_mutex_init(&(PRIVATE->mutex), &attr);
|
||||||
pthread_mutexattr_destroy(&attr);
|
pthread_mutexattr_destroy(&attr);
|
||||||
|
|||||||
Reference in New Issue
Block a user