[petsc-users] vector of struct

2014-03-24 Thread Dinesh Kumar
Hi, I am trying to implement a 3-D Surface Registration code using PetSc. I want to create a vector of structures i.e. struct Point { double x, y, z; }; Then create a PetSc vector that stores array of "Points types". Can someone point me to the right direction. regards --dinesh

Re: [petsc-users] Building PETSc with Intel mpi

2014-03-24 Thread Qin Lu
It turns out the machine's Linux/glibc versions are too old. The configure passed the error after I switched to a new machine.   Thanks a lot, Qin From: Barry Smith To: Qin Lu Cc: petsc-users Sent: Saturday, March 22, 2014 5:30 PM Subject: Re: [petsc-users

Re: [petsc-users] hypre-boomeramg precondition

2014-03-24 Thread Barry Smith
On Mar 24, 2014, at 2:54 PM, Qiyue Lu wrote: > Dear All: > > I am solving a series of matrix which are from the same FEA model but have > different number of meshes. In petsc, I use CR as the solver and > hypre-boomeramg as the preconditioner. It works for a small system with 2 > million DOF

[petsc-users] hypre-boomeramg precondition

2014-03-24 Thread Qiyue Lu
Dear All: I am solving a series of matrix which are from the same FEA model but have different number of meshes. In petsc, I use CR as the solver and hypre-boomeramg as the preconditioner. It works for a small system with 2 million DOFs. Also works for a system with 6 million DOFs. But always fail

[petsc-users] Local vecs of a DMComposite() in F90 code

2014-03-24 Thread Gautam Bisht
Hi, I'm trying to write a fortran code that uses TS with DMComposite for a multiphysics problem. How can I access local vectors created using a DMCreateGlobalVector() from a Fortran code? I tried DMCompositeGetLocalVectorsF90(), but I believe there isn't such a function. Also, what is the differe

Re: [petsc-users] gamg failure with petsc-dev

2014-03-24 Thread Stephan Kramer
On 21/03/14 11:34, Stephan Kramer wrote: On 21/03/14 04:24, Jed Brown wrote: Stephan Kramer writes: We have been having some problems with GAMG on petsc-dev (master) for cases that worked fine on petsc 3.4. We're solving a Stokes equation (just the velocity block) for a simple convection in a

Re: [petsc-users] reg: Eigenvectors of general complex eigen value problem

2014-03-24 Thread venkatesh g
Thanks a lot. I did what you told. It is working now like a charm. However, I have another query, which I will post separately. On Mon, Mar 24, 2014 at 3:21 PM, Jose E. Roman wrote: > > El 24/03/2014, a las 10:28, venkatesh g escribió: > > > Ok. I am doing the following > > > > 1. Running ".

Re: [petsc-users] reg: Eigenvectors of general complex eigen value problem

2014-03-24 Thread Jose E. Roman
El 24/03/2014, a las 10:28, venkatesh g escribió: > Ok. I am doing the following > > 1. Running "./ex7 -f1 A2 -f2 B2 -st_type sinvert -evecs VC" > 2. In MATLAB I do "vc=PetscBinaryRead('VC');" > > The output is one Eigenvector in VC which is 800x1. > > And I did the normalization vc=vc/max(vc

Re: [petsc-users] reg: Eigenvectors of general complex eigen value problem

2014-03-24 Thread venkatesh g
Ok. I am doing the following 1. Running "./ex7 -f1 A2 -f2 B2 -st_type sinvert -evecs VC" 2. In MATLAB I do "vc=PetscBinaryRead('VC');" The output is one Eigenvector in VC which is 800x1. And I did the normalization vc=vc/max(vc) and compared with the normalized original vector.. Also I must tell

Re: [petsc-users] reg: Eigenvectors of general complex eigen value problem

2014-03-24 Thread Jose E. Roman
El 24/03/2014, a las 05:59, venkatesh g escribió: > Ok. I used -st_type sinvert only, > > I am comparing it eigs after normalization. Even now the eigenvectors differ. > > Whether I should use a different eps solver ? > > Pls let me know. > It works for me, so you should provide more det