On Jul 28, 2011, at 10:57 PM, Debao Shao wrote:
> Dear all:
>
> I'm using PETSC iterative solver(PCILU & KSPGMRES) in OOQP, but when the to
> be solved matrix is ill-posed, say condition number ~ 1e+4, then, KSPGMRES
> will exceed its maximal iteration number(default as 1), while when I
>
6 15 6 15
>>>> 7 16 7 16
>>>> 8 17 8 17
>>>> 90 90
>>>> 101101
>>>> 112112
>>>> 123123
>>>> 134134
>>>> 145145
>>>> 156156
>>>> 167167
>>>> 178178
>>>>
>>>>
>>>> which is simply the initial ordering with a change in the order in which
>>>> processors handle nodes. Could it be that the partitioning is not unique
>>>> and each time the algorithm only tries to obtain the "best" possible
>>>> ordering depending on the initial distribution? If so, how should I know
>>>> what ordering to start with?
>>>>
>>>
>>> Yes, ParMetis does not provide a unique "best" ordering, which is at
>>> least NP-complete if not worse.
>>>
>>>Matt
>>>
>>>
>>>> I am really confused and would appreciate if someone could provide some
>>>> insights.
>>>>
>>>> Thanks,
>>>> Mohammad
>>>>
>>>
>>>
>>>
>>> --
>>> What most experimenters take for granted before they begin their
>>> experiments is infinitely more interesting than any results to which their
>>> experiments lead.
>>> -- Norbert Wiener
>>>
>>
>>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110728/c9fe82af/attachment-0001.htm>
rder in which
>> processors handle nodes. Could it be that the partitioning is not unique
>> and each time the algorithm only tries to obtain the "best" possible
>> ordering depending on the initial distribution? If so, how should I know
>> what ordering to start with?
>>
>
> Yes, ParMetis does not provide a unique "best" ordering, which is at least
> NP-complete if not worse.
>
>Matt
>
>
>> I am really confused and would appreciate if someone could provide some
>> insights.
>>
>> Thanks,
>> Mohammad
>>
>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which their
> experiments lead.
> -- Norbert Wiener
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110728/81e6f267/attachment-0001.htm>
Dear all:
I'm using PETSC iterative solver(PCILU & KSPGMRES) in OOQP, but when the to be
solved matrix is ill-posed, say condition number ~ 1e+4, then, KSPGMRES will
exceed its maximal iteration number(default as 1), while when I checked the
same data with cholmod direct solver, it can sol
Dear all:
I'm using PETSC iterative solver(PCILU & KSPGMRES) in OOQP, but when the to be
solved matrix is ill-posed, say condition number ~ 1e+4, then, KSPGMRES will
exceed its maximal iteration number(default as 1), while when I checked the
same data with cholmod direct solver, it can solv
n HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110728/f8715986/attachment-0001.htm>
Thanks for your reply, Shri. Responses below:
On 2011-07-28, at 1:03 PM, Shri wrote:
>
> - Original Message -
>> Hi Barry, Shri,
>>
>> I reworked my implementation to include only the time-independent
>> electrical problem, and to make use of SNESVISetVariableBounds().
>> However, it do
Dear Hong,
Yes, there is no memory corruption for small size problem test after I add
those options. I have not tested on big size problems yet.
For using superlu_dist as a direct solver, I would like to make sure that the
there is no memory corruption for large size problem, i.e., the memory
Rebecca:
Turn off orderings and some options, e.g.,
-mat_superlu_dist_equil NO -mat_superlu_dist_rowperm NATURAL
-mat_superlu_dist_colperm NATURAL
Do you still get memory corruption?
Hong
> Hello all,
>
> I tried to use superlu as a direct solver running on Hopper, but found that
> there are s
write-strings -Wno-strict-aliasing -Wno-unknown-pragmas -g3 -o ex29
>> >>> ex29.o ComputeRHS.o
>> >>> -L/home/zlwei/soft/mercurial/petsc-dev/arch-linux2-c-debug/lib
>> -lpetsc
>> >>>
>> -Wl,-rpath,/home/zlwei/soft/mercurial/petsc-dev/arch-linux2-c-debug/lib
>> >>> -lflapack -lfblas -lm -L/usr/lib/gcc/x86_64-redhat-linux/4.1.2 -ldl
>> -lmpich
>> >>> -lopa -lmpl -lrt -lpthread -lgcc_s -lmpichf90 -lgfortran -lm -lm -ldl
>> >>> -lmpich -lopa -lmpl -lrt -lpthread -lgcc_s -ldl
>> >>> gcc: ex29.o: No such file or directory
>> >>> make: [ex29] Error 1 (ignored)
>> >>>
>> >>> The program and modified makefile are attached. Could you please take
>> a
>> >>> look and give me some suggestions.
>> >>>
>> >>> best,
>> >>> Alan
>> >>
>> >
>> >
>>
>>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110728/3eab77dd/attachment.htm>
On Jul 28, 2011, at 12:44 PM, Adam Byrd wrote:
> Dear All,
>
> I'm attempting to find the most efficient/simple way to collectively retrieve
> just a handful of values (local and non local) from a parallel matrix. I'm
> filling two 4x4 arrays with 4 values from 8 different rows of the solution
- Original Message -
> Hi Barry, Shri,
>
> I reworked my implementation to include only the time-independent
> electrical problem, and to make use of SNESVISetVariableBounds().
> However, it does not work quite as I expected: Constraints on
> variables that are functions of other variable
anks.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110728/4327ab15/attachment.htm>
Hello all,
I tried to use superlu as a direct solver running on Hopper, but found that
there are some memory corruption errors:
x/xyuan> cd $PBS_O_WORKDIR
Directory:
/global/homes/x/xyuan/Workspace_Nersc/cartmhdpdslin/trunk/test_superlu_as_direct_solver/m256_p1024
test_superlu_as_direct_solver/
Ping :
I tested similar calling procedure using
petsc-dev/src/ksp/ksp/examples/tutorials/ex5.c
successfully.
Try following:
1) switch to petsc-dev. See
http://www.mcs.anl.gov/petsc/petsc-as/developers/index.html on how to
get it.
2) send us a short and stand-alone code that reproduce the error. We
-gnu-intel all
xxx=xxx
it is odd.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110728/464c53e4/attachment.htm>
mercurial/petsc-dev/arch-linux2-c-debug/lib
> >>> -lflapack -lfblas -lm -L/usr/lib/gcc/x86_64-redhat-linux/4.1.2 -ldl
> -lmpich
> >>> -lopa -lmpl -lrt -lpthread -lgcc_s -lmpichf90 -lgfortran -lm -lm -ldl
> >>> -lmpich -lopa -lmpl -lrt -lpthread -lgcc_s -ldl
> >>> gcc: ex29.o: No such file or directory
> >>> make: [ex29] Error 1 (ignored)
> >>>
> >>> The program and modified makefile are attached. Could you please take a
> >>> look and give me some suggestions.
> >>>
> >>> best,
> >>> Alan
> >>
> >
> >
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110728/4826a630/attachment.htm>
Hi Barry, Shri,
I reworked my implementation to include only the time-independent electrical
problem, and to make use of SNESVISetVariableBounds(). However, it does not
work quite as I expected: Constraints on variables that are functions of other
variables have no effect on the variables by wh
Hello Matt,
"Undefined reference" is usually a linker error. As you noticed in the
console normally a two step process is applied: First the .c file(s)
is(are) compiled into .o files. Then they are linked. As Berend pointed
out the error happened in the treatment of the .c file. So why is it
linki
Dear Matt,
I can't directly see anything wrong with your settings.
Have you tried compiling your example outside of eclipse with the PETSc
makefile system?
The previous pastebin output you emailed me does suggest that the error
is during the compile stage: it cannot find the correct headers
d
Dear Matt,
What do you have under
Properties->C/C++ Build -> Settings -> GCC Linker Libraries for both
"Libraries" and "Library Search Path"?
Regards,
Berend.
On 07/28/2011 02:04 AM, Matt Bockman wrote:
> Hi Berend,
>
> Under Properties->C/C++ Build -> Settings -> GCC C Compiler ->
> Directori
>>> -Wl,-rpath,/home/zlwei/soft/mercurial/petsc-dev/arch-linux2-c-debug/lib
>>> -lflapack -lfblas -lm -L/usr/lib/gcc/x86_64-redhat-linux/4.1.2 -ldl
-lmpich
>>> -lopa -lmpl -lrt -lpthread -lgcc_s -lmpichf90 -lgfortran -lm -lm -ldl
>>> -lmpich -lopa -lmpl -lrt -lpthread -lgcc_s -ldl
>>> gcc: ex29.o: No such file or directory
>>> make: [ex29] Error 1 (ignored)
>>>
>>> The program and modified makefile are attached. Could you please take a
>>> look and give me some suggestions.
>>>
>>> best,
>>> Alan
>>
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110728/98cc332a/attachment.htm>
Dear Matt,
The error is definitely in the compiler stage: it cannot find that
appropriate petsc.h include file.
1. Have you included petsc.h in the source file?
2. What do you have under
-> Settings -> Includes
?
Also, just to make sure, you can set the appropriate values for
PETSC_DIR and PET
Dear Matt,
Does it say this during the compiling? Or linking?
If it says this during the compiling, it means that eclipse cannot find
the PETSc header files. So, it must be the setting of the "Includes".
You might want to "hard-code" the directory, just to make sure.
It is indeed not completel
Dear Matt,
I use Eclipse and have eclipse make the makefiles.
It is just a matter of indicating to eclipse where the PETSc
headers/libraries are to be found, so if you have a C project which
needs PETSc headers and libraries:
To do this, click on your managed C project with the right sided mous
25 matches
Mail list logo