[CIG-SHORT] writing to output_points

[Brad Aagaard]

Rowena,

In addition to a small bug in how we searched for points (see below), 
the origin of the problem was that due to rounding errors PyLith did not 
find the points (-40,10,0) and (-40,0,0) in the mesh. Moving these 
points epsilon below the ground surface solves this problem.

I have fixed the search algorithm so that it skips cohesive cells (which 
have zero volume) and will give the proper error message associated with 
a point not being found within the domain. This fix will appear in the 
next bugfix release.

Regards,
Brad


On 6/29/12 12:14 PM, Rowena Lohman wrote:
> Good afternoon,
>
> I'm trying to expand on the greens function example in the 2D/strike
> slip example directory, and have been successful except for one issue
> - writing output.  I switched to 3D for my example and believe I have
> made the correct changes to the rest of the cfg and fault impulse
> files since I get no errors until the end.  Note that the grid I'm
> using is pretty coarse at this point - I'm just trying to get it to
> run to completion.
>
> Here's what happens at the writing stage:  If I use a SMALL number of
> output points (15 or less), it works fine.  However, if I use any
> larger number I get the following error:
>
> mpinemesis(21991) malloc: *** error for object 0x1054a5c34: pointer
> being freed was not allocated
> *** set a breakpoint in malloc_error_break to debug
> [0]0:Return code = 0, signaled with Abort trap: 6
>
> I THINK I've included everything such that pylith
> --problem=pylith.problems.GreenFns should work.  Do I just need to
> figure out how to adjust the input cfg files to go back with just
> writing to one of the domains?  I can do the interpolation myself
> later, but thought it would be nice to use the pylith output points
> version if possible.
>
> Thanks!
> R
>
>
>
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>