Re: [Libmesh-users] Exporting assembled matrix into Petsc binary

2011-07-18 Thread Omar Al-Abbasi
Thanks a lot. -- AppSumo Presents a FREE Video for the SourceForge Community by Eric Ries, the creator of the Lean Startup Methodology on "Lean Startup Secrets Revealed." This video shows you how to validate your ideas,

Re: [Libmesh-users] Exporting assembled matrix into Petsc binary

2011-07-17 Thread Omar Al-Abbasi
Thanks Jed for the reply, I tried what you suggested but since the matrix that I have is called matrix_A, I tried the following: Mat A = matrix_A.mat(); matrix_A is created by SparseMatrix& matrix_A = *eigen_system.matrix_A; I got the following error error: ‘struct libMesh::SparseMatrix’ has n

[Libmesh-users] Exporting assembled matrix into Petsc binary

2011-07-16 Thread Omar Al-Abbasi
Hi everyone, Is it possible to export the stiffness matrix once it has been assembled into Petsc binary format? Regards, Omar -- AppSumo Presents a FREE Video for the SourceForge Community by Eric Ries, the creator of t

Re: [Libmesh-users] Setting a target eigenvalue rather than selecting the position of spectrum

2011-06-29 Thread Omar Al-Abbasi
Thanks for the suggestions. John are you suggesting that within my Libmesh code and prior to specifying the position of the spectrum I shall write ... (-st_shift .75) and then I specify the spectrum to LARGEST_MAGNITUDE ?? I would appreciate your help. Omar ---

[Libmesh-users] Setting a target eigenvalue rather than selecting the position of spectrum

2011-06-29 Thread Omar Al-Abbasi
Hi, I am using Libmesh with Slepc for solving the eigen system that I have. I know that in Slepc there is an option to search for eigenvalues around specific target (say eigenvalues around 0.75), so is there a way in Libmesh to specify this target rather than setting the position of the spectrum t

Re: [Libmesh-users] Defining the Dirichlet boundary conditions

2011-05-11 Thread Omar Al-Abbasi
e for Petsc 3.1. If you are using a different > version there will most likely be a slightly different way to tell it to > keep the nonzero pattern. I recommend setting this each time before you > start filling the matrix up. Petsc seems to like to forget about this > option. > >

Re: [Libmesh-users] Defining the Dirichlet boundary conditions

2011-05-10 Thread Omar Al-Abbasi
Thank you all for all the suggestions you made. In fact I tried doing two seperate sweeps as Boyce suggested, but it seems that I am doing a mistake which leads to no change in mass matrix. This how I am doing that after the first sweep for assembling the matrices: . . . for ( ; el != end_el; ++el

[Libmesh-users] Defining the Dirichlet boundary conditions

2011-05-10 Thread Omar Al-Abbasi
Hi everyone, I am trying to define the boundary conditions in a way which I am more familiar with. I have only one type of BC which is Drichilet and I am doing that simply by making indices that are on the outer surface of the domain to be zero in the stiffness (k) and mass (m) matrix and 1 for th

[Libmesh-users] Exporting mass and stiffness matrices?

2011-04-15 Thread Omar Al-Abbasi
Hi, Is it possible to export the mass and stiffness matrices that is created by Libmesh to any data file? Regards, Omar -- Benefiting from Server Virtualization: Beyond Initial Workload Consolidation -- Increasing the us

Re: [Libmesh-users] How to search for smallest eigenvalues?

2011-04-15 Thread Omar Al-Abbasi
About searching for the smallest eigenvalues. Ok I figured it out. In fact that option was there in example 17 but in a commented line. :-) On Thu, Apr 14, 2011 at 12:32 PM, Omar Al-Abbasi wrote: > Hi everyone, > > I am solving a generalized eigenvalue problem and I followed e

[Libmesh-users] How to search for smallest eigenvalues?

2011-04-14 Thread Omar Al-Abbasi
Hi everyone, I am solving a generalized eigenvalue problem and I followed example 17 successfully until now. My question is how to specify the option for SLEPc so it searches for the smallest eigenvalues. I know that this option is there but unfortunately I don't know how to implement it. Many th

Re: [Libmesh-users] Getting PETSc error after modifying the stiffness matrix

2011-04-13 Thread Omar Al-Abbasi
Thank you all. So it seems that the problem is memory allocation. Ok is there any suggestion to solve this problem. On Wed, Apr 13, 2011 at 5:42 PM, Derek Gaston wrote: > On Apr 13, 2011, at 3:32 PM, Roy Stogner wrote: > > > No; you're getting it because you're trying to write to unallocated >

[Libmesh-users] Getting PETSc error after modifying the stiffness matrix

2011-04-13 Thread Omar Al-Abbasi
Hi everyone, I am trying to add a term to a the stiffness matrix which is V and it appears as follows: for (unsigned int i=0; ihttp://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#Signal[0]PETSCERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory c