[petsc-dev] Issues with DMComplex Cone Orientations in Fortran

[Matthew Knepley]

On Tue, Aug 28, 2012 at 4:58 PM, Chris Eldred <chris.eldred at gmail.com>wrote:

> Building on the previous example code (for the Doublet mesh), I am
> getting an error when trying to set Cone orientations:
>
> [0]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [0]PETSC ERROR: Argument out of range!
> [0]PETSC ERROR: Cone orientation 0 is not in the valid range [-1. 0)!
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Development HG revision:
> b6fe82991deee4a0d3f9a20654bc7750a6b1fe0f  HG Date: Mon Aug 27 13:17:10
> 2012 -0500
> [0]PETSC ERROR: See docs/changes/index.html for recent updates.
> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [0]PETSC ERROR: See docs/index.html for manual pages.
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: ./test on a arch-linu named Puget-101334 by user Tue
> Aug 28 15:54:36 2012
> [0]PETSC ERROR: Libraries linked from
> /home/user/Desktop/LIBRARIES/petsc-dev/arch-linux2-cxx-debug/lib
> [0]PETSC ERROR: Configure run at Mon Aug 27 12:55:14 2012
> [0]PETSC ERROR: Configure options --download-boost --download-chaco
> --download-ctetgen --download-f-blas-lapack --download-fiat
> --download-generator --download-metis --download-ml --download-mpich
> --download-parmetis --download-scientificpython --download-triangle
> --with-clanguage=cxx --with-dynamic-loading --with-shared-libraries
> --with-sieve PETSC_ARCH=arch-linux2-cxx-debug
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
> [0]PETSC ERROR: DMComplexSetConeOrientation() line 1207 in
> /home/user/Desktop/LIBRARIES/petsc-dev/src/dm/impls/complex/complex.c
> application called MPI_Abort(MPI_COMM_WORLD, 63) - process 0
> [unset]: aborting job:
> application called MPI_Abort(MPI_COMM_WORLD, 63) - process 0
>
> The code is:
>
> Orientations(1) = 0
> Orientations(2) = 1
> pOrientations => Orientations
> call DMComplexSetConeOrientation(dm,4,pOrientations,ierr)
> CHKERRQ(ierr)
> call DMComplexGetConeOrientation(dm,4,pOrientations,ierr)
> CHKERRQ(ierr)
> write(*,*) pOrientations
>
> Mesh point 4 is an edge (so it has a cone of size 2)- It is unclear to
> me why the cone orientation range is [-1,0). My understanding of cone
>

This means that the number of dof on that point is 0.

   Matt


> orientations is that they specified the order in which cone's were
> traversed for a point- ie for mesh point 4 (an edge) I could give it a
> cone Orientation of (1,0) which would mean look at point 1 of the cone
> first, then point 0.
>
> --
> Chris Eldred
> DOE Computational Science Graduate Fellow
> Graduate Student, Atmospheric Science, Colorado State University
> B.S. Applied Computational Physics, Carnegie Mellon University, 2009
> chris.eldred at gmail.com
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120828/735affb1/attachment.html>