This is the mail archive of the gsl-discuss@sourceware.org mailing list for the GSL project.


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

- no subject -


I am a student who is trying to use the histogram methods that were
published , but my problem is that i have tried to get the logarithmic
values of these histograms using these 2 methods :

first method for normalization : 

ÂMake continuous double probability density functions by appropriate
normalization
Â*/
template <typename T, typename V>
gsl_histogram** CHistogram<T,
V>::getContinuousDoubleProbabilityDensities()
{
ÂÂÂint i;

ÂÂÂif ( this->doubleHistograms != NULL ) {
ÂÂÂ// Clone histograms and normalize
ÂÂÂif ( this->continuousDoubleProbabilityDensities == NULL ) {
ÂÂÂÂÂÂthis->continuousDoubleProbabilityDensities = new gsl_histogram*[
this->numberOfHistograms ];

ÂÂÂÂÂÂfor ( i = 0; i < this->numberOfHistograms; i++ ) {
ÂÂÂÂÂÂthis->continuousDoubleProbabilityDensities[ i ] =
gsl_histogram_clone( this->doubleHistograms[ i ] );
ÂÂÂÂÂÂgsl_histogram_scale( this->continuousDoubleProbabilityDensities[ i
],
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂ1.0 / ( gsl_histogram_sum(
this->continuousDoubleProbabilityDensities[ i ] )
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂ* this->binWidth ) );
ÂÂÂÂÂÂ}
ÂÂÂ}

ÂÂÂreturn this->continuousDoubleProbabilityDensities;
ÂÂÂ}
ÂÂÂelse {
ÂÂÂthis->printError( "getContinuousDoubleProbablilityDensities(),"
ÂÂÂÂÂÂÂÂÂÂ"No continuous double probability densities available!" );

ÂÂÂreturn NULL;
ÂÂÂ}
}

second method to get log : 

/*
Â* Get logarithmized continuous double probability density functions
Â*/
template <typename T, typename V>
gsl_histogram** CHistogram<T,
V>::getLogContinuousDoubleProbabilityDensities()
{
ÂÂÂint i, j;

ÂÂÂif ( this->doubleHistograms != NULL ) {
ÂÂÂ// Make continuous double probability density functions
ÂÂÂif ( this->continuousDoubleProbabilityDensities == NULL ) {
ÂÂÂÂÂÂgetContinuousDoubleProbabilityDensities();
ÂÂÂ}

ÂÂÂ// Create histograms and logarithmize content
ÂÂÂif ( this->logContinuousDoubleProbabilityDensities == NULL ) {
ÂÂÂÂÂÂthis->logContinuousDoubleProbabilityDensities = new
gsl_histogram*[ this->numberOfHistograms ];

ÂÂÂÂÂÂfor ( i = 0; i < this->numberOfHistograms; i++ ) {
ÂÂÂÂÂÂthis->logContinuousDoubleProbabilityDensities[ i ] =
gsl_histogram_alloc( this->numberOfBins );
ÂÂÂÂÂÂgsl_histogram_set_ranges_uniform(
this->logContinuousDoubleProbabilityDensities[ i ],
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂthis->minValue,
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂthis->maxValue );

ÂÂÂÂÂÂfor ( j = 0; j < this->numberOfBins; j++ ) {
ÂÂÂÂÂÂÂÂÂgsl_histogram_accumulate(
this->logContinuousDoubleProbabilityDensities[ i ],
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂthis->minValue + ( j + 0.5 ) * this->binWidth,
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂlog10( gsl_histogram_get(
this->continuousDoubleProbabilityDensities[ i ],
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂj ) ) );
ÂÂÂÂÂÂÂÂÂÂÂÂÂ//std::cout << "result newÂÂ "
<<logContinuousDoubleProbabilityDensities<<std::endl;



ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂ}
ÂÂÂÂÂÂ}
ÂÂÂ}

ÂÂÂreturn this->logContinuousDoubleProbabilityDensities;


ÂÂÂ}
ÂÂÂelse {
ÂÂÂthis->printError( "getLogContinuousDoubleProbabilityDensities()",
ÂÂÂÂÂÂÂÂÂÂ"No logarithmized continuous double probability densities
available!" );
ÂÂÂreturn NULL;
ÂÂÂ}
}

/*
Â* Get logarithmized continuous complex probability density functions
Â*/
template <typename T, typename V>
gsl_histogram2d** CHistogram<T,
V>::getLogContinuousComplexProbabilityDensities()
{
ÂÂÂint i, j, k;

ÂÂÂif ( this->complexHistograms != NULL ) {
ÂÂÂ// Make discrete complex probability density functions
ÂÂÂif ( this->discreteComplexProbabilityDensities == NULL ) {
ÂÂÂÂÂÂgetDiscreteComplexProbabilityDensities();
ÂÂÂ}

ÂÂÂ// Create histograms and logarithmize content
ÂÂÂif ( this->logDiscreteComplexProbabilityDensities == NULL ) {
ÂÂÂÂÂÂthis->logDiscreteComplexProbabilityDensities = new
gsl_histogram2d*[ this->numberOfHistograms ];

ÂÂÂÂÂÂfor ( i = 0; i < this->numberOfHistograms; i++ ) {
ÂÂÂÂÂÂthis->logDiscreteComplexProbabilityDensities[ i ] =
gsl_histogram2d_alloc( this->numberOfBins,
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂthis->numberOfBins );
ÂÂÂÂÂÂgsl_histogram2d_set_ranges_uniform(
this->logDiscreteComplexProbabilityDensities[ i ],
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂthis->minValue,
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂthis->maxValue,
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂthis->minValue,
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂthis->maxValue );

ÂÂÂÂÂÂfor ( j = 0; j < this->numberOfBins; j++ ) {
ÂÂÂÂÂÂÂÂÂfor ( k = 0; k < this->numberOfHistograms; k++ ) {
ÂÂÂÂÂÂÂÂÂgsl_histogram2d_accumulate(
this->logDiscreteComplexProbabilityDensities[ i ],
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂthis->minValue + ( j + 0.5 ) * this->binWidth,
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂthis->minValue + ( k + 0.5 ) * this->binWidth,
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂlog10( gsl_histogram2d_get(
this->discreteComplexProbabilityDensities[ i ],
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂj,
ÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂÂk ) ) );
ÂÂÂÂÂÂÂÂÂ}
ÂÂÂÂÂÂ}
ÂÂÂÂÂÂ}
ÂÂÂ}

ÂÂÂreturn this->logDiscreteComplexProbabilityDensities;
ÂÂÂ}
ÂÂÂelse {
ÂÂÂthis->printError( "getLogDiscreteComplexProbabilityDensities()",
ÂÂÂÂÂÂÂÂÂÂ"No logarithmized discrete complex probability densities
available!" );
ÂÂÂreturn NULL;
ÂÂÂ}
}


but then he results of the histogram are either negative infinty (which 
is ok ) or positive values which is not ok as my probabilities should 
range between negative infinity and zero after getting log ....so any 
idea about this ,please help 

Best Regards , 

Sherry 


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]