Changeset 27376


Ignore:
Timestamp:
11/10/22 09:31:23 (2 years ago)
Author:
Mathieu Morlighem
Message:

BUG: fixing size allocation

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/classes/FemModel.cpp

    r27375 r27376  
    15951595
    15961596        /*Now send and receive vector for nodes on partition edge*/
    1597         /*Now send and receive vector for vertices on partition edge*/
    15981597        IssmDouble       **send_buffers  = xNew<IssmDouble*>(num_procs);
    15991598        ISSM_MPI_Request  *send_requests = xNew<ISSM_MPI_Request>(num_procs);
    16001599        for (int rank = 0;rank<num_procs;rank++){
    16011600                #ifdef _HAVE_AD_
    1602                 send_buffers[rank] = xNew<IssmDouble>(this->vertices->Size(),"t"); //only one alloc, "t" is required by adolc
     1601                send_buffers[rank] = xNew<IssmDouble>(this->nodes->Size(),"t"); //only one alloc, "t" is required by adolc
    16031602                #else
    1604                 send_buffers[rank] = xNew<IssmDouble>(this->vertices->Size());
     1603                send_buffers[rank] = xNew<IssmDouble>(this->nodes->Size());
    16051604                #endif
    16061605                send_requests[rank] = ISSM_MPI_REQUEST_NULL;
    16071606        }
    1608         IssmDouble* recv_buffer = xNew<IssmDouble>(this->vertices->Size());
     1607        IssmDouble* recv_buffer = xNew<IssmDouble>(this->nodes->Size());
    16091608
    16101609        for(int rank=0;rank<num_procs;rank++){
Note: See TracChangeset for help on using the changeset viewer.