fortran - MPI_ALLGATHER error parallelizing code -

-

i'm trying parallelize following code.

subroutine log_likelihood(y, theta, lli, ll)     doubleprecision, allocatable, intent(in)    :: y(:)      doubleprecision, intent(in)                 :: theta(2)     doubleprecision, allocatable, intent(out)   :: lli(:)     doubleprecision, intent(out)                :: ll     integer                                     :: allocate (lli(size(y))) lli = 0.0d0 ll = 0.0d0 = 1, size(y)     lli(i) = -log(sqrt(theta(2))) - 0.5*log(2.0d0*pi) &            - (1.0d0/(2.0d0*theta(2)))*((y(i)-theta(1))**2) end ll = sum(lli) end subroutine log_likelihood

to this, i'm trying use mpi_allgather. code wrote

subroutine log_likelihood(y, theta, lli, ll)     doubleprecision, allocatable, intent(in)    :: y(:)      doubleprecision, intent(in)                 :: theta(2)     doubleprecision, allocatable, intent(out)   :: lli(:)     doubleprecision, intent(out)                :: ll     integer                                     :: i, size_y, diff size_y=size(y) allocate (lli(size_y)) !broadcasting call mpi_bcast(theta, 1, mpi_double_precision, 0, mpi_comm_world, ierr) call mpi_bcast(y, 1, mpi_double_precision, 0, mpi_comm_world, ierr)  ! determine how many points handle each proc points_per_proc = (size_y + numprocs - 1)/numprocs ! determine start , end index proc's points istart = proc_num * points_per_proc + 1 iend = min((proc_num + 1)*points_per_proc, size_y) diff = iend-istart+1 allocate(proc_contrib(istart:iend)) = istart, iend     proc_contrib(i) = -log(sqrt(theta(2))) - 0.5*log(2.0d0*pi) &                     - (1.0d0/(2.0d0*theta(2)))*((y(i)-theta(1))**2) end call mpi_allgather(proc_contrib, diff, mpi_double_precision, &                    lli, diff, mpi_double_precision, &                    mpi_comm_world, ierr) ll = sum(lli) end subroutine log_likelihood

when try run program, following error.

$ mpiexec -n 2 ./mle.x  fatal error in pmpi_allgather: internal mpi error!, error stack: pmpi_allgather(961)......: mpi_allgather(sbuf=0x7ff2f251b860, scount=1500000, mpi_double_precision, rbuf=0x7ff2f2ad5650, rcount=3000000, mpi_double_precision, mpi_comm_world) failed mpir_allgather_impl(807).:  mpir_allgather(766)......:  mpir_allgather_intra(560):  mpir_localcopy(357)......: memcpy arguments alias each other, dst=0x7ff2f2ad5650 src=0x7ff2f251b860 len=12000000  =================================================================================== =   bad termination of 1 of application processes =   exit code: 1 =   cleaning remaining processes =   can ignore below cleanup messages ===================================================================================

can please explain me i'm doing wrong?

thanks!

i able solve problem. serial , parallel versions of code available @ https://bitbucket.org/ignacio82/bhhh

subroutine log_likelihood(y, n, theta, lli, ll)     integer, intent(in)                         :: n     doubleprecision, intent(in)                 :: y(n)      doubleprecision, intent(in)                 :: theta(2)     doubleprecision, intent(out)                :: lli(n)     doubleprecision, intent(out)                :: ll     integer                                     ::  = istart, iend     proc_contrib(i-istart+1) = -log(sqrt(theta(2))) - 0.5*log(2.0d0*pi) &                     - (1.0d0/(2.0d0*theta(2)))*((y(i)-theta(1))**2) end  if ( mod(n,numprocs)==0 )     call mpi_allgather(proc_contrib, points_per_proc, mpi_double_precision, &                        lli, points_per_proc, mpi_double_precision, &                        mpi_comm_world, ierr)     else if (numprocs-1 == proc_num )      recvcounts(numprocs) = iend-istart+1     call mpi_allgatherv(proc_contrib, points_per_proc, mpi_double_precision, &                         lli, recvcounts, displs, mpi_double_precision, &                         mpi_comm_world, ierr)  end if  ll = sum(lli) end subroutine log_likelihood

Comments

Post a Comment

Popular posts from this blog

c++ - Creating new partition disk winapi -

-
i'am trying create new partition , mount volume new, thought createfile let me this, code : lpctstr lpfilename=l"\\\\.\\device\\harddisk0\\partition3"; handle handl=createfile( lpfilename, generic_read | generic_write, file_share_read | file_share_write, null, create_always, file_attribute_normal, null ); if (handl==invalid_handle_value) { qdebug()<<"handl invalid"<<" error"<<getlasterror();} bool success = definedosdevice(ddd_raw_target_path,l"i:",l"\\device\\harddisk0\\partition3"); if(!success) qdebug()<<" definedosdevice failed "<<getlasterror(); bflag = getvolumenameforvolumemountpoint( l"i:\\", // input volume mount point or directory /** u in dir
Read more

Android Prevent Bluetooth Pairing Dialog -

-
i'm developing internal application uses bluetooth printing. want bluetooth pairing occur without user input. have managed working trapping android.bluetooth.device.action.pairing_request broadcast. in broadcast receiver call setpin method, , pairing works ok, bluetoothpairingdialog displayed second or two, disappears - see link below. https://github.com/android/platform_packages_apps_settings/blob/master/src/com/android/settings/bluetooth/bluetoothpairingdialog.java since broadcast non-ordered, can't call abortbroadcast() , , wondering if there other way prevent pairing dialog appearing. can hook window manager in way? i haven't been able come way without modifying sdk. if you're oem, it's easy (i'm on 4.3): in packages/apps/settings/androidmanifest.xml, comment intent filter pairing dialog: <activity android:name=".bluetooth.bluetoothpairingdialog" android:label="@string/bluetooth_pairing_request"
Read more

VBA function to include CDATA -

-
i have following function include cdata: function cdatasection() dim objdom domdocument dim objkmlrootelement ikmldomelement dim objkmlelement ikmldomelement dim cdata ikmldomcdatasection set objdom = new domdocument set objkmlrootelement = objdom.createelement("balloonstyle") objdom.appendchild objkmlrootelement set objkmlelement = objdom.createelement("text") objkmlrootelement.appendchild objkmlelement set cdata = objdom.createcdatasection("text") cdata.data = "<![cdata[<b>latitude = $[latitude]</b>?]]>;" end function when run above, getting error "user defined data type not found" function. yes ripster says find library reference. vba window tools>references for list of relevant references see... http://msdn.microsoft.com/en-us/library/windows/desktop/ms763701%28v=vs.85%29.aspx however, if have hard time tracking down reference, replace dim sta
Read more