[Trilinos-Users] Multidimensional Arrays

[John Mitchell]

Thanks Chris.  I believe that's what I was looking for.

Christopher Siefert wrote:
> John,
>
> I believe an Epetra_BlockMap is the item you want, with elements of size
> 3.
>
> Epetra_BlockMap (int NumGlobalElements, int ElementSize, int IndexBase,
> const Epetra_Comm &Comm)
>
> -Chris
>
> On Mon, 2010-02-01 at 14:27 -0700, John Mitchell wrote:
>   
>> I have a question on multidimensional arrays with Epetra_Maps and 
>> Epetra_Vectors.
>>
>> For example, suppose I have a set of n points.  Associated with each 
>> point I have coordinates (x,y,z) and a global Id that I use to create an 
>> Epetra_Map.  For a given point local index j, I can look up the 
>> coordinates vec[3*j+k], k=0,1,2. 
>>
>> I would like to create an Epetra_Map that understands this and stores 
>> the data so that I don't have to create 3 Epetra_Vectors (for X, Y, Z). 
>>
>> Is this possible?  If not, what is the recommended best practice?
>>
>> I see that there is Epetra_MultiVector with the following constructor, 
>> but that doesn't appear to do what I want (unless I do a copy of all 
>> incoming data into an array double *A that is not interlaced and then 
>> use the constructor below to copy one more time)
>>
>> Epetra_MultiVector(Epetra_DataAccess CV, const Epetra_BlockMap& Map, 
>> double *A, int MyLDA, int NumVectors);
>>
>> Thanks,
>> John
>>
>>
>>
>> _______________________________________________
>> Trilinos-Users mailing list
>> Trilinos-Users at software.sandia.gov
>> http://software.sandia.gov/mailman/listinfo/trilinos-users
>>