:
[0]PETSC ERROR: User provided function() line 0 in unknown directory unknown
file
application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0
[unset]: aborting job:
application called MPI_Abort(MPI_COMM_WORLD, 59) - process 0
R
On Jan 26, 2012, at 9:05 PM, Xuefei (Rebecca
that can get rid of those zeros in the matrix generated based
on the stencil information?
Thanks very much!
Best regards,
Rebecca
On Jan 20, 2012, at 5:23 PM, Matthew Knepley wrote:
On Fri, Jan 20, 2012 at 7:07 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Here is the output
Hi, Jed,
Thanks very much!
Problem solved.
Cheers,
Rebecca
On Jan 23, 2012, at 1:55 PM, Jed Brown wrote:
-dm_preallocate_only
http://www.mcs.anl.gov/petsc/petsc-dev/docs/manualpages/DM/DMSetMatrixPreallocateOnly.html
On Mon, Jan 23, 2012 at 15:52, Xuefei (Rebecca) Yuan xyuan
20, 2012 at 7:07 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Here is the output:
On Jan 20, 2012, at 5:01 PM, Matthew Knepley wrote:
On Fri, Jan 20, 2012 at 6:55 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Hello Matt,
I tried several times for 3.1-p8 and dev version
On Jan 21, 2012, at 10:27 AM, Matthew Knepley wrote:
On Sat, Jan 21, 2012 at 12:19 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Hello,
After I was not able to get rid of the zeros in the matrix by
ierr = MatSetOption(jacobian, MAT_IGNORE_ZERO_ENTRIES,
PETSC_TRUE);CHKERRQ(ierr
Hello all,
This is a test for np=1 case of the nonzero structure of the jacobian matrix.
The jacobian matrix is created via
ierr =
DMDACreate2d(comm,DMDA_BOUNDARY_NONE,DMDA_BOUNDARY_NONE,DMDA_STENCIL_BOX,
parameters.mxfield, parameters.myfield, PETSC_DECIDE, PETSC_DECIDE, 4, 2, 0, 0,
);CHKERRQ(ierr);
ierr = MatSetOption(jacobian, MAT_IGNORE_ZERO_ENTRIES,
PETSC_TRUE);CHKERRQ(ierr);
nothing changed.
The version is 3.1-p8.
Thanks very much!
Best regards,
Rebecca
On Jan 20, 2012, at 4:09 PM, Matthew Knepley wrote:
On Fri, Jan 20, 2012 at 6:02 PM, Xuefei (Rebecca) Yuan
);
None of those works. What is wrong here?
Thanks,
R
On Jan 20, 2012, at 4:32 PM, Matthew Knepley wrote:
On Fri, Jan 20, 2012 at 6:28 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Hello Matt,
I have changed the code as
ierr = MatSetOption(jacobian, MAT_IGNORE_ZERO_ENTRIES
On Fri, Jan 20, 2012 at 6:55 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Hello Matt,
I tried several times for 3.1-p8 and dev version by putting MatSetOption
Are you sure your entries are exactly 0.0?
Yes, because I have looked at the output, they are 0.00e+00.
R
Here is the output:
On Jan 20, 2012, at 5:01 PM, Matthew Knepley wrote:
On Fri, Jan 20, 2012 at 6:55 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Hello Matt,
I tried several times for 3.1-p8 and dev version by putting MatSetOption
Are you sure your entries are exactly 0.0
:
On Fri, Jan 20, 2012 at 7:07 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Here is the output:
On Jan 20, 2012, at 5:01 PM, Matthew Knepley wrote:
On Fri, Jan 20, 2012 at 6:55 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Hello Matt,
I tried several times for 3.1-p8
I have replace INSERT_VALUES by ADD_VALUES,
ierr = MatSetValuesStencil(jacobian, 1, row, 6, col, v,
ADD_VALUES);CHKERRQ(ierr);
but still cannot get rid of those zeros...
Cheers,
R
On Jan 20, 2012, at 5:28 PM, Xuefei (Rebecca) Yuan wrote:
I did use
ierr = MatSetValuesStencil(jacobian
Hello,
I use ksp_view to get the information about the lu pc as follow:
PC Object: 1 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
matrix ordering: nd
factor fill ratio given 0, needed 0
Factored matrix follows:
Matrix
Hello all,
I have a piece of code that need to store the RHS vector b of the ksp object
via vecview.
However, the routine is like:
ierr = DMMGSetSNESLocal(dmmg, FormFunctionLocal,
FormJacobianLocal,0,0);CHKERRQ(ierr);
ierr = DMMGSolve(dmmg);CHKERRQ(ierr);
The rhs evaluation
Hello Matt,
On Jan 17, 2012, at 12:53 PM, Matthew Knepley wrote:
On Tue, Jan 17, 2012 at 2:51 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Hello all,
I have a piece of code that need to store the RHS vector b of the ksp object
via vecview.
However, the routine is like
-ksp_view_binary?
This will save all the ksp content in binary form (more than one set of
jacobian and rhs), if I only want to save the first initial Jacobian and rhs,
will let -ksp_max_it 1 -snes_max_it 1 give me that?
Thanks,
Rebecca
Hi Satish,
Thanks very much!
Best regards,
Rebecca
On Jan 12, 2012, at 11:09 AM, Satish Balay wrote:
On Thu, 12 Jan 2012, Xuefei (Rebecca) Yuan wrote:
Hello all,
Is there any sample matrix and rhs data files for the test of ex10.c?
You can get some sample datafiles [mat,vec
. After the solver works well and you know the matrix structure
should not change,
then use SAME_NONZERO_PATTERN for optimization.
Hong
On Aug 31, 2011, at 2:42 PM, Barry Smith wrote:
On Aug 31, 2011, at 4:40 PM, Xuefei (Rebecca) Yuan wrote:
Hello all,
Is the default nonzero pattern
Hello all,
Is the default nonzero pattern for SuperLU_dist is SAME_NONZERO_PATTERN?
Thanks!
Best,
Rebecca
,SAME_NONZERO_PATTERN);CHKERRQ(ierr);
ierr = DMMGSolve(dmmg);CHKERRQ(ierr);
This will specify the nonzero pattern as the same? How could I check if I have
set this up right?
Thanks,
Rebecca
On Aug 31, 2011, at 2:42 PM, Barry Smith wrote:
On Aug 31, 2011, at 4:40 PM, Xuefei (Rebecca) Yuan
Hello all,
Could I convert a C code using PETSc library to a UPC code using PETSc?
Thanks,
Xuefei (Rebecca) Yuan
Postdoctoral Fellow
Lawrence Berkeley National Laboratory
Tel: 1-510-486-7031
Hello,
I would like to print out the nonzero structure of the analytic Jacobian
through PETSc, is there any example on that?
Thanks very much!
Xuefei (Rebecca) Yuan
Postdoctoral Fellow
Lawrence Berkeley National Laboratory
Tel: 1-510-486-7031
is wrong here?
Thanks very much!
Xuefei (Rebecca) Yuan
Postdoctoral Fellow
Lawrence Berkeley National Laboratory
Tel: 1-510-486-7031
(): memory corruption:
0x01e31280 ***
Any idea what is wrong here?
Thanks very much!
Xuefei (Rebecca) Yuan
Postdoctoral Fellow
Lawrence Berkeley National Laboratory
Tel: 1-510-486-7031
Hello,
I would like to use SuperLU as a direct solver for a large size problem running
on multiple processors.
I have write the analytical Jacobian via FormJacobianLocal(), however, I do not
know how to set the direct solver.
Any thoughts?
Thanks a lot!
Rebecca
Dear all,
I have assembled a 2d matrix through the FormJacobianLocal() routine.
I would like to store this sparse MXN matrix A in the following form:
line1: M N numberofnonzeros
line2: i j value
line3: i j
Hi Jed,
I need to have this form for the simple input for PDSLin.
What do you mean by COO?
Any examples to show how to pull that new array for reading from the AIJ
structure?
Thanks a lot!
R
On Jul 14, 2011, at 11:28 AM, Jed Brown wrote:
On Thu, Jul 14, 2011 at 13:18, Xuefei (Rebecca
Dear Matt and Jed,
Thanks very much! That helps!
Best regards,
Rebecca
On Jul 14, 2011, at 11:36 AM, Matthew Knepley wrote:
On Thu, Jul 14, 2011 at 6:30 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Hi Jed,
I need to have this form for the simple input for PDSLin.
What do you
/impls/aij/mpi/mumps/mumps.c
which convert petsc seqaij and mpiaij to triples: row[], col[], val[] .
You may take a look at these routines.
Hong
On Thu, Jul 14, 2011 at 1:37 PM, Xuefei (Rebecca) Yuan xyuan at lbl.gov
wrote:
Dear Matt and Jed,
Thanks very much! That helps!
Best regards
29 matches
Mail list logo