[petsc-users] [RFC] Docs: TeX -> HTML

2016-07-23 Thread Marco Zocca
Dear all, following the discussion at PETSc'16, I have tried to render the TeX-based manual into HTML with latex2html [1] and pandoc [2] . Neither attempt was successful, because of the presence of certain external TeX packages used for rendering various custom aspects of the manual. There is

Re: [petsc-users] HDF5 and PETSc

2016-07-11 Thread Marco Zocca
ed me initially. So my question becomes: please synchronize the output of ./configure --help with manpages and pdf manual :) Thanks again, Marco On 11 July 2016 at 09:57, Marco Zocca wrote: > Good morning, > >Does the HDF5 functionality need to be explicitly requested at > confi

[petsc-users] HDF5 and PETSc

2016-07-11 Thread Marco Zocca
Good morning, Does the HDF5 functionality need to be explicitly requested at configure time? I just noticed that my default configuration on a single-node machine does not compile any relevant symbol. I do not have HDF5 installed on my system yet, but I assumed PETSc includes it by default, or

[petsc-users] Re. PETSc user meeting 2016

2016-07-02 Thread Marco Zocca
poster presenters, but first and foremost the PETSc team, the sponsors and our excellent host Karl. Hoping to meet you soon again, perhaps at PETSc'17, Kind regards, Marco Zocca https://github.com/ocramz/petsc-hs https://github.com/ocramz/petsc-hs-docker

[petsc-users] re. AIJ preallocation man page

2016-05-23 Thread Marco Zocca
Dear all, I'm struggling with this statement: "It is recommended that you call both of the above preallocation routines [MatSeqAIJSetPreallocation and MatMPIAIJSetPreallocation] for simplicity." (source: https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MATAIJ.html ) I saw that,

[petsc-users] petsc 3.7 and slepc 3.6.1 incompatibility

2016-04-27 Thread Marco Zocca
estimate for a compatible release of SLEPc ? Thank you and kind regards, Marco Zocca

[petsc-users] erratic bug with nested PETSc and SLEPc

2016-01-25 Thread Marco Zocca
Dear all, I have a simple code with a matrix filling, assembly and output to stdout. This in turn is wrapped within a SLEPc bracket, and that in turn in a PETSc bracket, both called with default options (`XInitializeNoArguments`). I run this on my laptop, MPI comm size == 1. Issue: _sometimes_ t

Re: [petsc-users] installation on cloud platform

2016-01-04 Thread Marco Zocca
Hi Satish, thank you for the input; I was really looking for something that lets one abstract out "where" the code lives, so as to possibly work both in a single-node and cluster setting. This is why a "container" approach sounds meaningful. Configure once, run many. For message-passing codes s

Re: [petsc-users] installation on cloud platform

2016-01-04 Thread Marco Zocca
--download-mpich > --download-fblaslapack --with-metis=1 --download-metis=1 --with-debugging=0 > $ export PETSC_DIR=/home/ubuntu/petsc-3.6.3 > $ export PETSC_ARCH=arch-linux2-c-opt > $ make all > $ export PATH=$PATH:$PETSC_DIR/$PETSC_ARCH/bin > > Tabrez > > > On 01/03/2016

[petsc-users] installation on cloud platform

2016-01-03 Thread Marco Zocca
Dear all, has anyone here tried/managed to install PETSc on e.g. Amazon AWS or the Google Compute Engine? I believe some extra components are needed for coordination, e.g. Kubernetes or Mesos (in turn requiring that the library be compiled within some sort of container, e.g. Docker), but I'm a

[petsc-users] LP solvers based on PETSc ?

2015-12-29 Thread Marco Zocca
Dear all, Is anyone here aware of linear program solvers, of any form, being implemented on top of PETSc ? Thank you in advance, Marco

Re: [petsc-users] Eigenvalue problem in batch petsc-users Digest, Vol 84, Issue 59

2015-12-18 Thread Marco Zocca
Hi Venkatesh, it sounds like a job for the batch scheduler; again, a "bash" problem, with some additional variables: On a RHEL cluster, I was using Torque/PBS, which gives the `qsub`, `qstat` commands. I would guess that Cray has analogous ones . PETSc and SLEPc are instead toolkits for splitti

[petsc-users] shared components for real- and complex-scalar builds

2015-11-16 Thread Marco Zocca
Currently, I'm re-building fblaslapack and mpich along with PETSc for both cases, but would it be possible to share either between the two builds? Thank you, Marco

[petsc-users] modeling a dual mesh

2015-11-14 Thread Marco Zocca
What construct can I use to build and keep track (i.e. map over) a staggered index set, i.e. the dual (vertices, elements) meshes for a FEM application. I am looking for the equivalent, on a regular mesh, of the `sieve` construction ( http://www.mcs.anl.gov/petsc/documentation/tutorials/sieve.pdf

Re: [petsc-users] on SNESSetJacobian semantics

2015-10-22 Thread Marco Zocca
Thank you, but this is not really helping; why are there two slots to supply the Jacobian of the same function, i.e. Amat and J ? On 22 October 2015 at 16:30, Matthew Knepley wrote: > On Thu, Oct 22, 2015 at 9:20 AM, Marco Zocca > wrote: > >> The signature of SNE

[petsc-users] on SNESSetJacobian semantics

2015-10-22 Thread Marco Zocca
The signature of SNESSetJacobian is PetscErrorCode SNESSetJacobian(SNES snes,Mat Amat,Mat Pmat,PetscErrorCode (*J)(SNES,Vec,Mat,Mat,void*),void *ctx) however it would seem redundant to supply both a constant matrix and a function to compute the Jacobian. The manual says of J "(if NULL then SNES

[petsc-users] [ANN] PETSc Haskell bindings

2015-10-04 Thread Marco Zocca
forward to all feedback; Kind regards, Marco Zocca

Re: [petsc-users] human-readable examples and PETSc idioms

2015-07-10 Thread Marco Zocca
Thank you for the reply; let me improve my question: in e.g. a FEM code, we need: *) a mesh table :: [element -> [face -> [edge -> [node] ] ] ] *) one or more reference elements with associated basis over elements and/or faces and Jacobians, or a quadrature rule to represent integration over

[petsc-users] human-readable examples and PETSc idioms

2015-07-09 Thread Marco Zocca
Hello, I'm looking for some non-optimized PETSc code; namely, I struggle a bit with generalizing the provided examples. On one hand inlined elementary operations help keep track of the FLOPs but make for hard-to read code. For instance, KSP example 2, ( http://www.mcs.anl.gov/petsc/petsc-curren