Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Data projection bug #149

Open
PrasadBabarendaGamage opened this issue Jun 2, 2018 · 0 comments
Open

Data projection bug #149

PrasadBabarendaGamage opened this issue Jun 2, 2018 · 0 comments
Labels

Comments

@PrasadBabarendaGamage
Copy link
Contributor

A bug seems to have slipped through after the recent merge.

    # Create the data points
    dataPoints = iron.DataPoints()
    dataPointUserNumber = 1
    numberOfDataPoints = 1
    numberOfDimensions = 3
    dataPoints.CreateStart(dataPointUserNumber, region, numberOfDataPoints)
    dataPointLocations = np.zeros((numberOfDataPoints, numberOfDimensions))
    print("Number of data points: " + str(numberOfDataPoints))
    dataPointLocations[0,:] = np.array([1,1,1])
    dataList = [1]
    for dataPointIdx, dataPoint in enumerate(dataList):
        dataPoints.PositionSet(dataPoint,dataPointLocations[dataPointIdx,:])
    dataPoints.CreateFinish()

    # Set up data projection
    dataProjectionUserNumber = 1
    dataProjection = iron.DataProjection()
    dataProjection.CreateStart(dataProjectionUserNumber, dataPoints,
                               geometricField, iron.FieldVariableTypes.U)
    dataProjection.AbsoluteToleranceSet(1.0e-14)
    dataProjection.RelativeToleranceSet(1.0e-14)
    dataProjection.MaximumNumberOfIterationsSet(int(1e9))
    dataProjection.ProjectionTypeSet(
        iron.DataProjectionProjectionTypes.ALL_ELEMENTS)
    dataProjection.CreateFinish()

    # Evaluate data projection based on geometric field
    dataProjection.DataPointsProjectionEvaluate(
        iron.FieldParameterSetTypes.VALUES)
At line 324 of file /opencmiss/src/iron/src/data_projection_routines.f90
Fortran runtime error: Index '0' of dimension 1 of array 'closestdistances' below lower bound of 1

This error does not occur on Alys's mac build but it occurs on the linux GNU build. Likely because something is not initialised correctly.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

1 participant