Skip to content

Bug?: S2Fun and densities #302

Closed
Closed
@kilir

Description

Hi Ralf,
using the mtex5.0 beta 2 version, I notice that for a close-to-uniform distribution on the sphere, I get density values of around 0.3 while with mtex4.5.2 they were approaching 1.

(For example plot(vector3d.rand(1e6),'contourf','halfwidth',30*degree,'minmax') or using calcDensity with 5.0 in contrast to kernelDensityEstimation in 4.5.2.)
The same applies to other distributions which by now have characteristically lower densities. However, pole figures plotted from a uniform ODF have a density of 1.

I'd assume that this is a bug or is there any reason for density estimates to come close to uniform to be at 0.3 or am I missing something right here?

The sample happens for calcPDF:

cs=crystalSymmetry('321')
odf =uniformODF(cs)
pf = calcPDF(odf,Miller(1,0,0,cs))

where the value is actually 0.28 and different form the 1 produced by plotPDF(odf,...).

Cheers,
Rüdiger

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions