, a las 16:35, Barry Smith escribió:
> >
> >
> >
> >> On Aug 30, 2020, at 2:57 AM, Eda Oktay wrote:
> >>
> >> Dear Matt,
> >>
> >> First of all I figured out that I asked wrongly. It's not ParMETIS
> giving the same result. It is CHA
Okay thank you so much!
Mark Adams , 30 Ağu 2020 Paz, 15:34 tarihinde şunu yazdı:
>
>>
>>
>> So, if ParMETIS gives different edge cut as it is expected,
>> MatPartitioningGetUseEdgeWeights and MatPartitioningSetUseEdgeWeights works
>> correctly. Why can't CHACO?
>>
>>>
>>>
> Chaco does not
ereas when the option changes to
chaco, index set remains the same.
Thanks!
Eda
Eda Oktay , 30 Ağu 2020 Paz, 11:05 tarihinde şunu
yazdı:
> And is edge weights being less than 10 still valid? Although I am not
> getting any errors, when the elements in my matrix are larger than 10, even
> Pa
And is edge weights being less than 10 still valid? Although I am not
getting any errors, when the elements in my matrix are larger than 10, even
ParMETIS doesn't give different results.
Eda Oktay , 30 Ağu 2020 Paz, 10:57 tarihinde şunu
yazdı:
> Dear Matt,
>
> First of all I figured o
, at 2:23 PM, Eda Oktay wrote:
> >
> > Hi all,
> >
> > I am trying to partition a sparse matrix by using ParMETIS. I am
> converting my matrix to adjacency type and then applying partitioning.
>
> You don't need to do this. Just pass your original matri
Hi all,
I am trying to partition a sparse matrix by using ParMETIS. I am converting
my matrix to adjacency type and then applying partitioning. Default, I
understood that partitioning doesn't use edge-weights. However, when I used
the following codes I saw from ex15 and used
onstrates the problem?
>
> Barry
>
>
> > On Aug 14, 2020, at 8:53 AM, Eda Oktay wrote:
> >
> > Hi all,
> >
> > I am trying to try something. I am using the same MatPartitioning codes
> for both CHACO and ParMETIS:
> >
> > ierr =
> MatC
Hi all,
I am trying to try something. I am using the same MatPartitioning codes for
both CHACO and ParMETIS:
ierr = MatConvert(SymmA,MATMPIADJ,MAT_INITIAL_MATRIX,);CHKERRQ(ierr);
ierr = MatPartitioningCreate(MPI_COMM_WORLD,);CHKERRQ(ierr);
ierr =
>
> Or you can ask sysadmin on this machine to install python3-distutils or
> equivalent package.
>
> Satish
>
> On Fri, 31 Jul 2020, Eda Oktay wrote:
>
> > Hi all,
> >
> > I configured petsc-3.13.2 with following options:
> >
> > ./configure --wi
ed, 29 Jul 2020, Satish Balay via petsc-users wrote:
>
> > I can reproduce this - by uninstalling system make.
> >
> > Will check.
> >
> > Satish
> >
> > On Wed, 29 Jul 2020, Eda Oktay wrote:
> >
> > > Dear Barry,
> > >
> >
Dear Satish,
Thank you so much for your replies. I installed git and I was able to
configure PETSc.
Thanks again!
Eda
Satish Balay , 30 Tem 2020 Per, 18:32 tarihinde şunu yazdı:
>
> On Thu, 30 Jul 2020, Eda Oktay wrote:
>
> > Dear Satish,
> >
> > I am sorry but
Dear Satish,
I am sorry but I couldn't understand what should I do. Am I getting this
error because of case insensitivity?
If I am, why now? I mean I configured Petsc with exact same options before
in 2 machines. What can be the reason?
Thanks!
Eda
On Thu, Jul 30, 2020, 6:21 PM Satish Balay
en/maint should resolve
> the problem. It will force the cmake build to use the make obtained with
> --download-make
>
>Barry
>
>
> On Jul 28, 2020, at 9:48 AM, Matthew Knepley wrote:
>
> On Tue, Jul 28, 2020 at 10:45 AM Eda Oktay wrote:
>>
>> I am t
y@pj01 petsc]$ file x.tar.gz
> x.tar.gz: ASCII text, with no line terminators
> [balay@pj01 petsc]$
>
> On Mon, 27 Jul 2020, Satish Balay via petsc-users wrote:
>
> > Looks like python urlretrive() is failing on glitab tarball urls [but its
>
21 Jul 2020, at 11:06, Eda Oktay wrote:
> >
> > Dear Lawrence,
> >
> > I am using MPICC but not Mac, Fedora 25. If it will still work, I will
> try that.
> >
> > Thanks!
>
> It might be the case. When you observe the error, does "nslookup
> localhost" take a long time?
>
> Lawrence
Dear Lawrence,
I am using MPICC but not Mac, Fedora 25. If it will still work, I will try that.
Thanks!
Eda
Lawrence Mitchell , 21 Tem 2020 Sal, 12:58 tarihinde
şunu yazdı:
>
>
>
> > On 21 Jul 2020, at 10:49, Eda Oktay wrote:
> >
> > Hi all,
> >
>
Hi all,
I am using the following libraries and for some reason, I figured out
that if ı am disconnected to internet, my program is not working:
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
ose
>
>
> > El 19 jul 2020, a las 19:49, Matthew Knepley escribió:
> >
> > On Sun, Jul 19, 2020 at 11:01 AM Eda Oktay wrote:
> > Hi all,
> >
> > I am trying to install petsc 3.13.2 to my computer (fedora 25) with
> > slepc and some other librarie
Okay thanks again!
On Fri, Jul 10, 2020, 7:52 PM Jose E. Roman wrote:
> I don't know. This question belongs to the application part, not the
> solver part.
>
> > El 10 jul 2020, a las 18:45, Eda Oktay escribió:
> >
> > I looked at a2965*2965 sized matrix. I wa
hods work, iteration gets in
> action when the matrix size is larger than the subspace size (ncv). The
> tolerance is relevant when the solver iterates, not when it computes the
> solution with LAPACK.
>
> Jose
>
>
> > El 10 jul 2020, a las 18:32, Eda Oktay escribió:
> &
No, it is of size 9
Jose E. Roman , 10 Tem 2020 Cum, 19:31 tarihinde şunu yazdı:
>
> [Please respond to the list.]
>
> Is your matrix of size 8? This would explain the residuals.
>
> Jose
>
> > El 10 jul 2020, a las 17:10, Eda Oktay escribió:
> >
Hello,
I am currently using EPS to find 2 smallest eigenvalues and
corresponding eigenvectors of a matrix in order to be used in spectral
graph partitioning algorithms. I know that default tolerance is E-8
and I also tried the solver with tolerances E-2, E-4 and E-6 expecting
that in every matrix
;
> If inside the problem you can use -start_in_debugger to get the true
> location of the crash, if in a bash script you can use
> -on_error_attach_debugger
>
> Barry
>
>
> > On Jul 6, 2020, at 1:14 PM, Eda Oktay wrote:
> >
> > Hello,
> >
> > I
Hello,
I am trying to run my codes 100 times and then take the average
elapsed time. Although my program runs correctly for single time, when
I put them into loop, I am getting the following error:
*** Error in `./avg_coloredge_deneme_new_vecscatter_arastep':
munmap_chunk(): invalid pointer:
Thanks!
Jose E. Roman , 6 Tem 2020 Pzt, 12:10 tarihinde şunu yazdı:
>
> The index i should be a value between 0 and nconv-1
> https://slepc.upv.es/documentation/current/docs/manualpages/EPS/EPSGetEigenpair.html
> Jose
>
> > El 6 jul 2020, a las 11:06, Eda Oktay escribió:
>
Hello,
I am using EPSGetEigenpair to find some matrices' eigenpairs:
ierr = EPSGetType(eps,);CHKERRQ(ierr)
ierr = EPSGetDimensions(eps,,NULL,NULL);CHKERRQ(ierr);
ierr = MatCreateVecs(L,,NULL);CHKERRQ(ierr);
Vec *V;
VecDuplicateVecs(vr,nev,);
for (i=0;
.)
>
> On Jun 10, 2020, at 8:08 AM, Eda Oktay wrote:
>
> Dear Matt,
>
> Matthew Knepley , 10 Haz 2020 Çar, 16:03 tarihinde
> şunu yazdı:
>
>
> On Wed, Jun 10, 2020 at 8:56 AM Eda Oktay wrote:
>
>
> Hi all,
>
> I am trying to get all the rows of a parallel m
I got it, thanks so much!
Eda
On Wed, Jul 1, 2020, 2:42 PM Matthew Knepley wrote:
> On Wed, Jul 1, 2020 at 7:32 AM Eda Oktay wrote:
>
>> I'm sorry but I still can't understand how to put an input as these
>> parameters. It is the first time I saw 'const' for IS. If i
type? It is MPI
Dense.
Thanks!
Eda
Matthew Knepley , 1 Tem 2020 Çar, 13:43 tarihinde şunu yazdı:
>
> On Wed, Jul 1, 2020 at 6:34 AM Eda Oktay wrote:
>>
>> Der Barry,
>>
>> ı am trying to use your way but I couldn't understand how I create sub
>> matric
tors, they can't share the space with the matrix.
>
> Barry
>
>
>
> On Jun 10, 2020, at 11:36 AM, Matthew Knepley wrote:
>
> On Wed, Jun 10, 2020 at 12:26 PM Eda Oktay wrote:
>>
>> Matthew Knepley , 10 Haz 2020 Çar, 19:13 tarihinde
>> şun
derstanding what you are aiming to do, or if I am
> misunderstanding something.
>
> It sounds like you would benefit from having some routines in PETSc to do
> k-means (or other) clustering, by the way?
>
> Best regards,
> Richard
>
> On 4/29/20 3:47 AM, Eda Oktay
Oh okay, I misunderstood before. I will try that. Thanks!
Matthew Knepley , 22 Haz 2020 Pzt, 19:10 tarihinde
şunu yazdı:
>
> On Mon, Jun 22, 2020 at 11:39 AM Eda Oktay wrote:
>>
>> Matthew Knepley , 22 Haz 2020 Pzt, 14:43 tarihinde
>> şunu yazdı:
>> >
>>
Matthew Knepley , 22 Haz 2020 Pzt, 14:43 tarihinde
şunu yazdı:
>
> On Mon, Jun 22, 2020 at 7:14 AM Eda Oktay wrote:
>>
>> Hello everyone,
>>
>> I am trying to find elements in off diagonal blocks of a parallel
>> sparse matrix. That is why, I want to use MatGet
Hello everyone,
I am trying to find elements in off diagonal blocks of a parallel
sparse matrix. That is why, I want to use MatGetDiagonalBlock and from
the matrix I obtain, I want to obtain off diagonal elements by
subtracting it from my original matrix by using MatAXPY.
However, since
t communicators. For example, one is created with MPI_COMM_WOLD and
> the other with PETSC_COMM_WORLD.
>
> On Fri, Jun 19, 2020 at 12:48 PM Eda Oktay wrote:
>>
>> Dear Mark,
>>
>> I updated Petsc and again, here's the error:
>>
şunu yazdı:
>
> I don't know what is going on here. There was an update to this function
> about a year ago, so that might fix your problem.
>
> We would need you to test with a current version.
>
> Mark
>
> On Fri, Jun 19, 2020 at 11:23 AM Eda Oktay wrote:
>>
>
Hi all,
I am trying to find off block diagonal entries of a matrix and I am
trying to use MatFindOffBlockDiagonalEntries. However, although my
matrix is not NULL, I am getting an error message as follows:
1]PETSC ERROR: - Error Message
Hello everyone,
I am solving a graph partitioning problem. I found an unnormalized
Laplacian matrix of a graph, then since I have 4 processes and I am
using spectral partitioning algorithm, I calculated 4 eigenvectors
corresponding to 4 smallest eigenvalues of this Laplacian matrix.
Then, since
Matthew Knepley , 10 Haz 2020 Çar, 19:13 tarihinde
şunu yazdı:
>
> On Wed, Jun 10, 2020 at 12:07 PM Eda Oktay wrote:
>>
>> Der Matt,
>>
>> When I looked at the results, I found that there are some problems I
>> couldn't understand.
>>
>> First of
2020 Çar, 18:11 tarihinde
şunu yazdı:
>
> On Wed, Jun 10, 2020 at 10:09 AM Eda Oktay wrote:
>>
>> Dear Matt,
>>
>> I have one last question I believe. Up to creating a dense matrix I
>> did what you've suggested. Thank you so much for that.
>>
>> I cr
as rows?
Thanks!
Eda
Matthew Knepley , 10 Haz 2020 Çar, 16:16 tarihinde
şunu yazdı:
>
> On Wed, Jun 10, 2020 at 9:08 AM Eda Oktay wrote:
>>
>> Dear Matt,
>>
>> Matthew Knepley , 10 Haz 2020 Çar, 16:03 tarihinde
>> şunu yazdı:
>> >
>> > On Wed, Ju
Dear Matt,
Matthew Knepley , 10 Haz 2020 Çar, 16:03 tarihinde
şunu yazdı:
>
> On Wed, Jun 10, 2020 at 8:56 AM Eda Oktay wrote:
>>
>> Hi all,
>>
>> I am trying to get all the rows of a parallel matrix as individual
>> vectors. For instance, if I have 72*4 m
Hi all,
I am trying to get all the rows of a parallel matrix as individual
vectors. For instance, if I have 72*4 matrix, I want to get 72
different vectors having size 4.
As far as I understood, MatGetRow is only for local rows, so
MatGetOwnershipRange is used, however, when I tried this one, I
MATLAB k-means function, calling it like
>
> idx = kmeans(X,k)
>
> will give you the index set, but if you do
>
> [idx,C] = kmeans(X,k)
>
> then you will also get a matrix C which contains the cluster centroids. Is
> this not what you need?
>
> --Richard
>
>
I am sorry, I used VecDuplictaeVecs not MatDuplicateVecs
Eda Oktay , 22 May 2020 Cum, 20:31 tarihinde şunu
yazdı:
> Dear Richard,
>
> Thank you for your email. From MATLAB's kmeans() function I believe I got
> the final clustering index set, not centroids. What I am trying to do is
misunderstanding something.
>
> It sounds like you would benefit from having some routines in PETSc to do
> k-means (or other) clustering, by the way?
>
> Best regards,
> Richard
>
> On 4/29/20 3:47 AM, Eda Oktay wrote:
>
> Dear Richard,
>
> I am tryi
PETSc. I was thinking that I should implement some of these
> routines using PETSc. I can't get to this immediately, but I'm wondering
> if you might care to tell me a bit more about the clustering problems
> you need to solve and how having some support for this in PETSc might
> (or m
Hi all,
I created a parallel vector UV, by using VecDuplicateVecs since I need row
vectors of a matrix. However, I need the whole vector be in all processors,
which means I need to gather all and broadcast them to all processors. To
gather, I tried to use VecStrideGatherAll:
Vec UVG;
Hello all,
I am trying to form a matrix whose columns are eigenvectors I have
calculated before U = [v1,v2,v3]. Is there any easy way of forming this
matrix? My matrix should be parallel and I have created vectors as below,
where nev i s the number of requested eigenvalues. So each V[i]
Hi all,
I solved the problem by using different routine.
Thanks anyway!
Eda
Eda Oktay , 18 Ara 2019 Çar, 13:21 tarihinde şunu
yazdı:
> Hi all,
>
> I am trying to partition a matrix in multiple processes. However, at some
> point I get the following eror:
>
> [ 0] ***ASSERT
Hi all,
I solved the problem by using different routine.
Thanks anyway!
Eda
Eda Oktay , 18 Ara 2019 Çar, 10:24 tarihinde şunu
yazdı:
> Hello everyone,
>
> In a part of my code, I am trying to understand whether a matrix is equal
> to its transpose, that is why I am using MatEqua
Hi all,
I am trying to partition a matrix in multiple processes. However, at some
point I get the following eror:
[ 0] ***ASSERTION failed on line 176 of file
/home/edaoktay/petsc-3.11.1/arch-linux2-c-debug/externalpackages/git.parmetis/libparmetis/comm.c:
j == nnbrs
[ 2] [ 3] ***ASSERTION
Hello everyone,
In a part of my code, I am trying to understand whether a matrix is equal
to its transpose, that is why I am using MatEqual(). However, for most of
the matrices I took from University of Florida State Sparse Matrix Library,
I get the following error for more than 2 processes,
Hello all,
I am trying to find edge cut information of ParMETIS and CHACO. When I use
ParMETIS,
MatPartitioningViewImbalance(part,partitioning)
works and it gives also number of cuts.
However, when I used CHACO, it only gives imbalance information, not edge
cut. I have index sets but I couldn't
Hello all,
I am partitioning some graphs and I need to visualize them. I was wondering
if I can do it in PETSc. I know that there is some Draw routines but most
of them are for line graphs. ı didn't use PETSc for drawing before so I
don't know how and what to do.
Thanks!
Eda
Yes I use -mat_partitioning chaco.
Hello everyone,
I am trying to partition a matrix by using chaco. By using ex15 from mat
examples, I wrote the code below:
MatPartitioningCreate(MPI_COMM_WORLD, );
MatPartitioningSetAdjacency(part, A);
MatPartitioningSetFromOptions(part);
MatPartitioningApply(part, );
ISView(is,
Dear Matt,
You were right, I confused concepts. By setting vertex weights, I am able
to get different number of edge cuts.
Thanks,
Eda
Matthew Knepley , 1 Tem 2019 Pzt, 13:10 tarihinde şunu
yazdı:
> On Mon, Jul 1, 2019 at 3:00 AM Eda Oktay wrote:
>
>> Dear Matt,
>>
>>
, Jun 28, 2019, 1:49 PM Matthew Knepley wrote:
> On Fri, Jun 28, 2019 at 5:38 AM Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
>
>> Hello everyone,
>>
>> I am trying to change load balance in ParMETIS in order to give different
>> num
Hello everyone,
I am trying to change load balance in ParMETIS in order to give different
number of elements to different processors. In the manual, it states that
changing "ubvec" parameter should work. So, instead of 1.05, I gave 1.2,1.5
and so on but load balance didn't change. I always get
Hello everyone,
I am trying to partition a matrix into unequal parts: for example I am
trying to partition a 10*10 matrix into 4*4 and 6*6 submatrices. However,
ParMETIS is always dividing the matrix into 2 equal parts (5*5).
The reason why I am trying to do this is that I am using spectral
h process the
> "diagonal" block and the "off-diagonal" block.
>
>
>
> > On Jun 17, 2019, at 7:44 AM, Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
> >
> > Hello everyone,
> >
> > I am trying to find the off diagonal
I did but ı didn't send you that part:
PetscInt kk;
for (i=0;i, 17 Haz 2019 Pzt, 16:23 tarihinde şunu yazdı:
> You also do not seem to be initializing kk.
>
> On Mon, Jun 17, 2019 at 4:00 AM Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
>
>> He
Hello everyone,
I am trying to find the off diagonal blocks of a matrix. For example, I
partitioned a 10*10 matrix into 4*4 and 6*6 block sub matrices. I need to
find the other blocks in order to find the number of nonzero elements in
them.
I tried to use MatGetLocalSubMatrix but I couldn't
Hello everyone again,
I was making an index mistake in for loop. I corrected it and my problem
solved.
Thank you,
Eda
Eda Oktay , 17 Haz 2019 Pzt, 10:28 tarihinde şunu
yazdı:
> Hello everyone,
>
> I am trying to permute a matrix by using a sorted eigenvector in order to
> partitio
Hello everyone,
I am trying to permute a matrix by using a sorted eigenvector in order to
partition the matrix. I sorted the vector and got the index set idx. I am
using 2 processors and I have to divide the vector into 2 according to the
sign of elements (first negative signed ones, then
2019 Çar, 14:51 tarihinde şunu
yazdı:
> On Wed, May 15, 2019 at 7:35 AM Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
>
>> Hello,
>>
>> I am trying to divide a matrix into unequal sized parts into different
>> processors (for example I want
Hello,
I am trying to divide a matrix into unequal sized parts into different
processors (for example I want to divide 10*10 matrix into 4*4 and 6*6
submatrix in two processors). When my program reads a matrix from file, it
automatically divides it into equal parts and then I can't change local
ttps://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind
> >
> >
> >
> >> On Apr 22, 2019, at 7:31 AM, Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
> >>
> >> Hello,
> >>
> >> I am trying to partit
Hello everyone,
I am trying to run my program in petsc 3.11.0. It is currently working in
petsc 3.10.3.
I am using slepc in the program, so when I used 3.10.3 the slepc library is
in petsc-3.10.3/arch-linux2-c-debug. However, in 3.11.0, it is only in
ov> wrote:
>
>> Dump the index sets and the matrix in binary and send them
>>
>
> First, reduce your problem size to about 10.
>
> Matt
>
>
>> Il giorno mer 24 apr 2019 alle ore 13:21 Eda Oktay
>> ha scritto:
>>
>>> Since I am using squa
/petsc/petsc/src/dba4d5898a4236b42fbe8dff209f1243d2f2582a/src/mat/impls/aij/mpi/mpiaij.c#lines-1652
> The code assumes the number of leaves is A->rmap->n (the row local size of
> the matrix). Shouldn't this be the local size of the IS (rowp)? Barry?
>
> Il giorno mer 24 apr 2019
giorno lun 22 apr 2019 alle ore 11:28 Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> ha scritto:
>
>> Thank you for your answers. I figured out that for a matrix of size
>> 4253*4253, local size of Mat is 2127 and 2127 for 2 processors. However, I
>> wrote the p
Hello,
I am trying to partition an odd-numbered sized (for example 4253*4253),
square permutation matrix by using 2 processors with ParMETIS. The
permutation matrix is obtained by permuting the matrix by an index set "is"
(MatPermute(A,is,is,)). I checked the index set, it gives a permutation
and
ion and 2) that the size of the IS on
> each process matches the number of rows on that process.
>
> Note from the manual page: The index sets should be on the same
> communicator as Mat and have the same local sizes.
>
> Barry
>
>
> > On Apr 8, 2019, at 3:21 AM, Ed
tire
error message to petsc-ma...@mcs.anl.gov--
Thanks!
Eda
Eda Oktay , 25 Mar 2019 Pzt, 13:53 tarihinde şunu
yazdı:
> I attached whole program I wrote where the problem is in line 285. One of
> the matrices I used was airfoil1_binary, included in the folder. Also, I
> included makefil
gt;
>
> > On Apr 5, 2019, at 7:20 AM, Matthew Knepley via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
> >
> > On Fri, Apr 5, 2019 at 3:20 AM Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
> > Hello,
> >
> > I am trying to ca
Hello,
I am trying to calculate unweighted Laplacian of a matrix by using 2 cores.
If the size of matrix is in even number then my program works. However,
when I try to use a matrix having odd number for size, I guess since size
of the matrix cannot be divided into processors correctly, I get the
Hi Mark,
siz%size = 0 in my case because siz=362 and I use at most 2 processors.
Mark Adams , 2 Nis 2019 Sal, 15:12 tarihinde şunu yazdı:
>
>
> On Tue, Apr 2, 2019 at 6:55 AM Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
>
>> Hi Barry,
>>
>
unu
yazdı:
>
>
> > On Mar 27, 2019, at 7:33 AM, Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
> >
> > Hello,
> >
> > I am trying to permute a matrix A(of size 2n*2n) according to a sorted
> eigenvector vr (of size 2n) in paralle
Hello,
I am trying to permute a matrix A(of size 2n*2n) according to a sorted
eigenvector vr (of size 2n) in parallel using 2 processors (processor
number can change). However, I get an error in MatPermute line stating that
arguments are out of range and a new nonzero caused a malloc even if I
end in a small example that we can debug.
>
> Thanks,
>
> Matt
>
> On Mon, Mar 25, 2019 at 12:38 AM Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
>
>> Hello,
>>
>> I am trying to permute a vector A using following lines:
>
Hello,
I am trying to permute a vector A using following lines:
ierr =
ISCreateGeneral(PETSC_COMM_SELF,siz,idx,PETSC_COPY_VALUES,);CHKERRQ(ierr);
ierr = ISSetPermutation(is);CHKERRQ(ierr);
ierr = ISDuplicate(is,);CHKERRQ(ierr);
ierr =
INAL.c'
>
> > El 20 mar 2019, a las 8:25, Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> escribió:
> >
> > Hello,
> >
> > I am trying to compile a parallel program DENEME_TEMIZ_ENYENI_FINAL.c in
> PETSc. I wrote the following makefile but it
Hello,
I wrote a code computing element wise absolute value of a matrix. When I
run the code sequentially, it works. However, when I try to use it in
parallel with the same matrix, I get the following error:
[1]PETSC ERROR: Argument out of range
[1]PETSC ERROR: Only local rows
The absolute
Hello,
I am trying to compile a parallel program DENEME_TEMIZ_ENYENI_FINAL.c in
PETSc. I wrote the following makefile but it says that there is nothing to
do with the program:
export CLINKER = gcc
DENEME_TEMIZ_ENYENI_FINAL: DENEME_TEMIZ_ENYENI_FINAL.o chkopts
-${CLINKER} -o
Dear Professor Roman,
I decided to install petsc-3.10.3 with slepc since I did it once
succesfully. Thanks for your help.
Eda
Jose E. Roman , 19 Mar 2019 Sal, 15:15 tarihinde şunu
yazdı:
> What is the output of 'make check' in $PETSC_DIR ?
>
> > El 19 mar 2019, a las 13:02, Eda Okt
This is slepc.log:
Checking environment... done
Checking PETSc installation...
ERROR: Unable to link with PETSc
ERROR: See "arch-linux2-c-debug/lib/slepc/conf/configure.log" file for
details
Matthew Knepley , 19 Mar 2019 Sal, 13:36 tarihinde şunu
yazdı:
> On Tue, Mar 19, 2019 a
Hello,
I am trying to install PETSc with following configure options:
./configure --download-openmpi --download-openblas --download-slepc
--download-cmake --download-metis --download-parmetis
Compilation is done but after the following command, I got an error:
make
Hello,
I have a Laplacian matrix PL of matrix A and I try to partition A using
PARMETIS. Since PL is sequential and not adjacency matrix, I converted PL
to AL, then write the following code:
ierr = MatConvert(PL,MATMPIADJ,MAT_INITIAL_MATRIX,);CHKERRQ(ierr);
ierr =
in
/home/edaoktay/petsc-3.10.3/arch-linux2-c-debug/externalpackages/git.slepc/src/eps/interface/epssolve.c
I understood from this error that the matrix does not converge. However, I
tested in MATLAB and it converges. The matrix is true now, so this cannot
be from the matrix.
Thanks,
Eda
Eda Oktay , 11
Dear Matt,
I printed in wrong state, ncols gives right solution.
But I still can't understand the first problem.
Eda
Eda Oktay , 11 Mar 2019 Pzt, 16:05 tarihinde şunu
yazdı:
> Dear Matt,
>
> Thank you for answering. First of all, sizeof(vals) returns to number of
> entrie
Hello,
I have a following part of a code which tries to change the nonzero values
of matrix L with -1. However in MatSetValues line, something happens and
some of the values in matrix turns into 1.99665e-314 instead of -1. Type of
arr is defined as PetscScalar and arr is produced correctly. What
,
Eda
Matthew Knepley , 8 Mar 2019 Cum, 16:41 tarihinde şunu
yazdı:
> On Fri, Mar 8, 2019 at 8:38 AM Eda Oktay wrote:
>
>> Dear Matt,
>>
>> Thank you for your answer but I get an error even before compiling ex6.
>> The error is:
>>
>> error: too many arg
Thu, Mar 7, 2019 at 6:45 AM Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> wrote:
>
>> Hello,
>>
>> I have a code finding Laplacian of a matrix and then I need to find this
>> Laplacian's second smallest eigenpair. I used EPSGetEigenpair code but I
>&
Hello,
I have a code finding Laplacian of a matrix and then I need to find this
Laplacian's second smallest eigenpair. I used EPSGetEigenpair code but I
get the values like "0." or "4." (I could not get decimals) when I used
PetscPrintf(PETSC_COMM_WORLD," The second smallest eigenvalue: %g\n",kr)
; easy way is to read them into Matlab with mmread.m and write with PETSc's
> PetscBinaryWrite.m. This can be done similarly in python.
> Jose
>
>
> > El 12 feb 2019, a las 8:50, Eda Oktay via petsc-users <
> petsc-users@mcs.anl.gov> escribió:
> >
> > Hello,
&g
Hello,
I am trying to use MatTranspose function as in below (last row) but I get
an error like this:
[0]PETSC ERROR: - Error Message
--
[0]PETSC ERROR: Object is in wrong state
[0]PETSC ERROR: Not for unassembled
ite? Does this for loop give me these 3 arrays (arrays consist
of row indices,column indices and values)?
Eda
Matthew Knepley , 25 Şub 2019 Pzt, 18:37 tarihinde şunu
yazdı:
> On Mon, Feb 25, 2019 at 10:20 AM Eda Oktay wrote:
>
>> Thank you for your answer Matt, but I have some
scViewerDestroy
> <https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Viewer/PetscViewerDestroy.html#PetscViewerDestroy>.
> Do I write these functions into for loop or not?
>
> Eda
>
> Matthew Knepley , 25 Şub 2019 Pzt, 17:31 tarihinde
> şunu yazdı:
>
>&
1 - 100 of 108 matches
Mail list logo