Re: [Libmesh-users] Stiffness assembling is very slow at the first time

2008-03-22 Thread Shengli Xu
Hi Ben and Roy In libMesh 0.5.0, The assemble function runs very quickly. But in libMesh0.6.2, It is very slowly. The memory used in libMesh0.5.0 is much larger than in libMesh0.6.2 before call assemble function. The program is the same. Why? I debugged my programe. The function of apply_periodic_b

Re: [Libmesh-users] Stiffness assembling is very slow at the first time

2008-03-20 Thread Benjamin Kirk
> You must be adding the constraint rows after the matrix has been > preallocated... Regardless of whether PeriodicBoundary calls would > fix your problem, we ought to make sure that users are able and > encouraged to do efficient preallocation for arbitrary user > constraints, too. The DOF conne

Re: [Libmesh-users] Stiffness assembling is very slow at the first time

2008-03-19 Thread Shengli Xu
2008/3/19, Roy Stogner <[EMAIL PROTECTED]>: > > > On Wed, 19 Mar 2008, Shengli Xu wrote: > > > > In the system, there are peroidic boundary conditions applied through " > > dof_map.constrain_element_matrix_and_vector(Ke,Fe,dof_indices);" before > > system.matrix->add_matrix(Ke,dof_indices); > > sys

Re: [Libmesh-users] Stiffness assembling is very slow at the first time

2008-03-19 Thread Roy Stogner
On Wed, 19 Mar 2008, Shengli Xu wrote: > In the system, there are peroidic boundary conditions applied through " > dof_map.constrain_element_matrix_and_vector(Ke,Fe,dof_indices);" before > system.matrix->add_matrix(Ke,dof_indices); > system.rhs->add_vector(Fe,dof_indices); > during elements asse

Re: [Libmesh-users] Stiffness assembling is very slow at the first time

2008-03-19 Thread Shengli Xu
Hi Ben and Roy, The version of libMesh is 0.6.2. I modified src/mesh/unstructured_mesh.C, add "(name.rfind(".unv") < name.size())" on line 385 to skip renumber nodes and elements for .unv input file. I use PETSc2.3.3. Running the code with -info, I find that Number of mallocs during MatSetValues()

Re: [Libmesh-users] Stiffness assembling is very slow at the first time

2008-03-19 Thread Roy Stogner
On Wed, 19 Mar 2008, Shengli Xu wrote: > micro_fluid_loadcase11_assemble is the performance log of calling the > function of assembling stiffness matrix. > | micro_fluid_loadcase11_assemble1 5515.6511 5515.651089 You might want to do some more fine-grained performance logging.

Re: [Libmesh-users] Stiffness assembling is very slow at the first time

2008-03-19 Thread Benjamin Kirk
> Hi, Libmesh Users > I find the stiffness matrix assembling is very slow at the first time. But > it is very quickly at the second time and laters for the same system. > It needs 2 hours to assemble the stiffness matrix of the 3d stokes system > with 15,000 dofs. How to speed up the process of ass

Re: [Libmesh-users] Stiffness assembling is very slow at the first time

2008-03-19 Thread Shengli Xu
Hi John The program runns in opt mode. The information of system is as follows: Mesh Information: mesh_dimension()=3 spatial_dimension()=3 n_nodes()=1331 n_elem()=1000 n_local_elem()=1000 n_active_elem()=1000 n_subdomains()=1 n_processors()=1 processor_id()=0 EquationSystems

Re: [Libmesh-users] Stiffness assembling is very slow at the first time

2008-03-18 Thread John Peterson
It should not take 2 hours to assemble a system with 15,000 dofs. Something is seriously wrong, perhaps you are in debug mode? -J On Tue, Mar 18, 2008 at 9:02 PM, Shengli Xu <[EMAIL PROTECTED]> wrote: > Hi, Libmesh Users > I find the stiffness matrix assembling is very slow at the first time. Bu

[Libmesh-users] Stiffness assembling is very slow at the first time

2008-03-18 Thread Shengli Xu
Hi, Libmesh Users I find the stiffness matrix assembling is very slow at the first time. But it is very quickly at the second time and laters for the same system. It needs 2 hours to assemble the stiffness matrix of the 3d stokes system with 15,000 dofs. How to speed up the process of assembling th