Create vectors

[Matthew Knepley]

Yes, Vec is just a regular type.

  Matt

On Mon, Nov 9, 2009 at 7:29 AM, <jarunan at ascomp.ch> wrote:

>
> I am solving multi-level grid (similar to multi grid but not the same). In
> each iteration, each level is solved separately but solutions are mapped to
> eachother. Each level has different size of matrix and vector. And each test
> case has different numbers of grid level.
>
> I have a difficulty to create vectors and matrices for each level,
> preparing for the computation, as I do not want to create and destroy them
> in every iteration. I am thinking of something similar to array pointer (the
> code is in fortran) e.g.,
>
> Type(real_array), Dimension(:), allocatable:: pointername
> Allocate(pointername(level_numbers))
>
> do i = 1, level_numbers
>   allocate(pointername(i)%p(size))
> enddo
>
>
> Is it possible to create pointer to vectors?
>
>
> Thank you
> Jarunan
>
>
> Quoting jarunan at ascomp.ch:
>
>
>> Hello,
>>
>> Is there an equivalent way to allocating array pointer for creating
>> vectors (or vector pointers)?
>>
>>
>> Jarunan
>>
>
>
>
> --
> Jarunan Panyasantisuk
> Development Engineer
> ASCOMP GmbH, Technoparkstr. 1
> CH-8005 Zurich, Switzerland
> Phone : +41 44 445 4072
> Fax   : +41 44 445 4075
> E-mail: jarunan at ascomp.ch
> www.ascomp.ch
>



-- 
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/20091109/3ffd9151/attachment.htm>