parallel processing - Calling a gfortran subroutine under OpenMP -


i have problem getting subroutine work under openmp. variable declarations (private, shared, firstprivate) seem okay, code works when subroutine statements written in line code. created toy problem demonstrates problem, proga.f works, while progb.f not. common region both programs is:

    program proga/b   use modprog     integer kmp     double precision dgam,gamma(10),phi(10),pphi(10),nderiv(10,10)     open(20,file='a/b.dat')     nn=10     dgam=0.1d-6     10 i=1,nn 10   gamma(i)=dfloat(i)     call gamset(gamma)     call errvec(phi)

the parallel region of (which works) is

!$omp parallel shared(nn,dgam,phi,nderiv)  !$omp& private(pphi,cm) !$omp& firstprivate(gamma)      100 kmp=1,nn     gamma(kmp)=gamma(kmp)+dgam     15 i=1,nn 15   cm(i)=gamma(i)     20 i=1,nn 20   pphi(i)=cm(i)**2      40 j=1,nn 40   nderiv(j,kmp)=(pphi(j)-phi(j))/dgam     gamma(kmp)=gamma(kmp)-dgam 100   continue  !$omp end parallel

while parallel region of progb (which doesn't work) 

!$omp parallel shared(nn,dgam,phi,nderiv)  !$omp& private(pphi,cm) !$omp& firstprivate(gamma)      100 kmp=1,nn     gamma(kmp)=gamma(kmp)+dgam     call gamset(gamma)     call errvec(pphi)     40 j=1,nn 40   nderiv(j,kmp)=(pphi(j)-phi(j))/dgam     gamma(kmp)=gamma(kmp)-dgam 100   continue  !$omp end parallel

the common subroutines , module 2 programs 

    subroutine gamset(gamma)     use modprog     double precision gamma(10)     15 i=1,nn 15   cm(i)=gamma(i)     return     end  c             subroutine errvec(phi)     use modprog     double precision phi(10)     20 i=1,nn 20   phi(i)=cm(i)**2      return     end  c     module modprog     integer nn     double precision cm(10)     end module


Comments

Post a Comment

Popular posts from this blog

javascript - Create a stacked percentage column -

i have collection of data : var data = [ {"p301a":"10","p301b":"7","p301c":"7","p301d":"3","p301e":"8","p301f":"8","p301g":"8","p301h":"8","p301i":"8","p301j":"8","p301k":"8","p301l":"8","p301m":"8","p301n":"8","p301o":"8","age":"31-40 years","profesion":"2","position":"2"}, {"p301a":"5","p301b":"4","p301c":"4","p301d":"4","p301e":"4","p301f":"4","p301g":"4","p301h":"4","p301i":"4","p301j":"4","p301k":"4...
Read more

Optimising Firebase database by automatically overwriting data -

in firebase app have following structure: posts latest all users can write posts > latest . want limit entries 20 posts save space in database. , posts beyond 20 unnecessary never shown on homepage. how can limit posts 20, when user writes posts > latest posts @ end drop off (be deleted) automatically? this sample should you. you need : const max_log_count = 20; exports.removeold = functions.database.ref('/posts/latest/{postid}').oncreate(event => { const parentref = event.data.ref.parent; return parentref.once('value').then(snapshot => { if (snapshot.numchildren() >= max_log_count) { let childcount = 0; const updates = {}; snapshot.foreach(function(child) { if (++childcount <= snapshot.numchildren() - max_log_count) { updates[child.key] = null; } }); // update parent. removes children. ...
Read more

javascript - Angular UI-Grid customTemplate directive causing rows to load slowly/? -

Image
each grid cell can have more 1 customtemplate loaded. have directive holds watcher on model_col_field listen cell updates on scroll when dom rows updated , recompile celltemplate . without watcher , $compile rows repeat on scroll, showing same set of rows. is there more angular way of approaching or way improve performance? relevant controller code: create gridoptions.data & gridoptions.columndefs : var customtemp = '<div ng-repeat="opt in grid.appscope.colvalueoptionsmodel" ng-init="tempstring=model_col_field.clazz+opt.key;"><comparison-directive dbo="model_col_field" attobj="opt" fff="grid.appscope.getattributetemplate(tempstring)" /></div>'; $scope.gridoptions.columndefs = [ { name:'yindex',pinnedleft:true,celltemplate: customtemp, displayname:'index'}, { name:'col1',celltemplate: customtemp}, { name:'col2',celltemplate: c...
Read more