parallel processing - Reshape an outer layer of a 3D array to 1D array -

-

i have 3d array, msh(0:m+1,0:n+1,0:l+1) representative of cartesian mesh grid. reshape outer layer 1d array. best way (efficient memory point of view)?

currently reshape plane plane starting -z,+z,:

array(1:m*n)=reshape(msh(1:m,1:n,0),m*n)    array(m*n+1:2*m*n)=reshape(msh(1:m,1:n,l+1),m*n)

my first question weather reshaping being done in memory efficient way! or there other way using loop faster , more memory efficient (the size of arrays huge).

secondly, if have several of these arrays in different processors best way gather these outer layers in master node (e.g. rank=0) using mpi.

it's o(2*m*n+2*l*m+2*l*n) problem way it, one option series of loops:

allocate(array(2*m*n+2*l*m+2*l*n)) h = 1 j=1,n    i=1,m       array(h) = msh(i,j,0)       array(h+m*n+1) = msh(i,j,l+1)       h = h+1    enddo enddo k=1,l    i=1,m       array(h) = msh(i,0,k)       array(h+m*l+1) = msh(i,n+1,k)       h = h + 1    enddo enddo k=1,l    j=1,n       array(h) = msh(0,j,k)       array(h+l*n+1) = msh(m+1,j,k)       h = h + 1    enddo enddo

then use normal mpi subroutines pass array master.


Comments

Post a Comment

Popular posts from this blog

python - Subclassed QStyledItemDelegate ignores Stylesheet -

-
i'm trying subclass qstyleditemdelegate remove focus rectangle in qcombobox . even though i'm calling base implementation of paint function , nothing else, result different. looks if parts of style-sheet influence bounding box of item taken consideration. class pstyleditemdelegate(qstyleditemdelegate): def __init__(self, *args, **kwds): super(pstyleditemdelegate, self).__init__(*args, **kwds) def paint(self, *args, **kwargs): qstyleditemdelegate.paint(*args, **kwargs) what have paint unmodified qstyleditemdelegate ? as suggested tried replacing pyside pyqt4 , works, appears bug. update pyside 1.1.2 1.2.1 result same. unfortunately switch breaks other parts of code, if there no other suggestions i'm going accept answer. edit bug tracked here
Read more

java - HttpClient 3.1 Connection pooling vs HttpClient 4.3.2 -

-
i have used httpclient3.1(java) connection pooling below : import org.apache.commons.httpclient.httpclient; import org.apache.commons.httpclient.multithreadedhttpconnectionmanager; import org.apache.commons.httpclient.params.httpclientparams; import org.apache.commons.httpclient.params.httpmethodparams; import org.apache.log4j.logger; public class httpclientmanager { private static logger logger = logger.getlogger(httpclientmanager.class); private static multithreadedhttpconnectionmanager connectionmanager = new multithreadedhttpconnectionmanager(); private static httpclientparams clientparams = new httpclientparams(); static { try { clientparams.makelenient(); clientparams.setauthenticationpreemptive(false); clientparams.setparameter(httpmethodparams.so_timeout, new integer(30000)); // set so_timeout clientparams.setparameter(httpmethodparams.head_body_check_timeout, new integer(3000)); // head responses
Read more

SQL: Divide the sum of values in one table with the count of rows in another -

-
is possible to, in single query, sum bunch of values in 1 table, , divide count of rows in another? there common id/key in both tables... fyi, win points least info provided in question. answer it's possible...here psuedo code since i'm guessing @ tables , columns select summing / counting whynot (select id, sum(whatever) summing where_ever ) inner join (select id, count(1) counting what_ever)b on a.id = b.id question = 'vague' group a.id
Read more