On Thu, Sep 19, 2013 at 11:44 AM, Reza Yaghmaie wrote:
>
> Dear Sir/Madam,
>
> I am very new to PETSC and I have some issues running my code.
> I want to create a global M*1 matrix (not a vector since I want to
> multiply this matrix with an other 1*M matrix and get a M*M matrix and I
> know that
Thanks Matthew for the reply.
It would be time consuming for me to implement an other way for creating a
matrix and doing matmult and I am not concerned with how much time this
will take.
So is there any way which you could possibly guide me through using
matmpibaij?
Thanks alot,
-Reza
On Thu, Sep 19, 2013 at 12:01 PM, Reza Yaghmaie wrote:
>
> Thanks Matthew for the reply.
> It would be time consuming for me to implement an other way for creating a
> matrix and doing matmult and I am not concerned with how much time this
> will take.
> So is there any way which you could possibl
I'm not sure I understand your statement about matrix shells being time
consuming to implement. They are easy, see this example:
http://www.mcs.anl.gov/petsc/petsc-current/src/ksp/ksp/examples/tutorials/ex14f.F.html
Randy M.
On Sep 19, 2013, at 12:01 PM, Reza Yaghmaie wrote:
>
> Thanks Mat
Dear Sir/Madam,
I am very new to PETSC and I have some issues running my code.
I want to create a global M*1 matrix (not a vector since I want to multiply
this matrix with an other 1*M matrix and get a M*M matrix and I know that
transpose function on vectors does not work so this is why I need to
El 31/07/2013, a las 12:03, Heikki Virtanen escribió:
> Hi, Sorry to bother you, but I still have problems with MATMPIBAIJ matrices.
> It is said here
>
> http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatCreateBAIJ.html
>
> that these matrices are created using block compress
> El 29/07/2013, a las 14:54, Heikki Virtanen escribió:
>
> > Hi, I try to solve an eigenvalue problem (Ax = lambda B x) where A
> > and B are complex matrices. Unfortunately, my data structures
> > are not capable of handling directly complex numbers, yet. So,
> > I have to use
> >
> >
El 29/07/2013, a las 14:54, Heikki Virtanen escribió:
> Hi, I try to solve an eigenvalue problem (Ax = lambda B x) where A
> and B are complex matrices. Unfortunately, my data structures
> are not capable of handling directly complex numbers, yet. So,
> I have to use
>
> [re(A
Hi, I try to solve an eigenvalue problem (Ax = lambda B x) where A
and B are complex matrices. Unfortunately, my data structures
are not capable of handling directly complex numbers, yet. So,
I have to use
[re(A) -im(A)]
RealA = []
[imA) r
