fortran - Reading multiple files in parallel and storing them in shared variables -

-

i have program reads number of large files (~100s files, 120mb each) generated mpi program, can take time. each file contains variables in corresponding subdomain. want read in variables files , store them specific slice of 4 dimensional array. since takes considerable amount of time, parallelize piece of code openmp:

  6     subroutine read_old_restart   7         integer :: ii   8         integer :: thread_id   9         integer :: omp_get_thread_num  10         character(len=21) :: file_name  11   12         !$omp parallel private(ii,file_name)  13         ii=0,nproc_old-1  14             if(ii < 10)  15                 write(file_name,401) "input/restart_00", ii, ".out"  16             else if(ii < 100)  17                 write(file_name,402) "input/restart_0" , ii, ".out"  18             else  19                 write(file_name,403) "input/restart_"  , ii, ".out"  20             end if  21             print*, "thread = ", omp_get_thread_num(), "reading ", file_name  22 401 format(a16,i1,a4)  23 402 format(a15,i2,a4)  24 403 format(a14,i3,a4)  25             open (unit=321, file=trim(file_name), status="old", form="unformatted")  26             read(321) t                       ,                         &   27                       old_u         (:,:,:,ii),                         &   28                       old_v         (:,:,:,ii),                         &   29                       old_w         (:,:,:,ii),                         &  30                       old_p         (:,:,:,ii),                         &  31                       old_uc        (:,:,:,ii),                         &  32                       old_vc        (:,:,:,ii),                         &  33                       old_wc        (:,:,:,ii),                         &  34                       old_un2       (:,:,:,ii),                         &  35                       old_vn2       (:,:,:,ii),                         &  36                       old_wn2       (:,:,:,ii),                         &  37                       old_un1       (:,:,:,ii),                         &  38                       old_vn1       (:,:,:,ii),                         &  39                       old_wn1       (:,:,:,ii),                         &  40                       old_p1        (:,:,:,ii),                         &  41                       old_viscu     (:,:,:,ii),                         &  42                       old_viscv     (:,:,:,ii),                         &  43                       old_viscw     (:,:,:,ii),                         &  44                       old_convu     (:,:,:,ii),                         &  45                       old_convv     (:,:,:,ii),                         &  46                       old_convw     (:,:,:,ii),                         &  47                       statindex               ,                         &  48                       old_umn       (:,:,:,ii),                         &  49                       old_uumn      (:,:,:,ii),                         &  50                       old_urms      (:,:,:,ii),                         &  51                       old_mass_frac (:,:,:,:,:,ii), &                     52                       old_entht     (:,:,:,:,ii)  53             close (321)  54         end  55         !$omp end parallel  56     end subroutine read_old_restart

the code compiles , runs fine first loop of each thread. here output:

 thread =            3 reading input/restart_030.out  thread =            7 reading input/restart_067.out  thread =            2 reading input/restart_020.out  thread =            6 reading input/restart_058.out  thread =            9 reading input/restart_085.out  thread =            8 reading input/restart_076.out  thread =            5 reading input/restart_049.out  thread =            4 reading input/restart_040.out  thread =           11 reading input/restart_103.out  thread =            0 reading input/restart_000.out  thread =            1 reading input/restart_010.out  thread =           10 reading input/restart_094.out

the code appears running , gets stuck on above output. when running top, cannot cpu usage. idea why not working expected?

you should use private integer variable unit number , set different value each thread. using same file unit differently different threads recipe trouble. quite surprised not crash.


Comments

Post a Comment

Popular posts from this blog

SVG stroke-linecap doesn't work for circles in Firefox? -

-
i've got problem svg stroke-linecap attribute. i've got circular progress bar in angularjs , set outer circle (blue one) have rounded "ends". @ fiddle . <svg ... height="130" width="130"> <!-- ngif: background --> <circle ... ng-if="background" fill="#fff" class="ng-scope" stroke-width="13" stroke="#cc3399" r="57.5" cy="65" cx="65" stroke-linecap="round" /> <!-- end ngif: background --> <circle ... fill="none" stroke-dashoffset="36.12831551628261" stroke-dasharray="361.28315516282623" stroke-width="13" stroke="#432db3" stroke-linecap="round" r="57.5" cy="65" cx="65" transform="rotate(-89.9, 65, 65)" /> </svg> how can that? ...
Read more

routes - Laravel 4 Wildcard Routing to Different Controllers -

-
realizing i'm going take beating this, have several questions "best practices" when constructing laravel route. do laravel routes run "first come first serve" top bottom? if route found, pass control controller (or perform action) , terminate? i have part of application acts url shortner. so, route looks this: route::get('/{id}', array('as' => 'home.text', 'uses' => 'homecontroller@text'), function($id){ //make sure id exists before pushing controller. if (storage::where('shortcode','=', $id)->count() > 0) { // route found } }); how can improved? one, i'm not sure how pass processing on controller conditionally.. sort of controller::make() ..? 3) there better way in laravel check if record exists versus doing mode::where('column','=','value')->count() 4) have legacy urls http://url.com/?...
Read more

cross browser - XSLT namespace-alias Not Working in Firefox or Chrome -

-
i transforming xhtml xhtml, need xslt stylesheet part of resulting document (the stylesheet contained in <script type"text/template"> element. using the xsl:namespace-alias instruction, works fine in ie, fails in both chrome , firefox. here relevant code: <xsl:output doctype-system="about:legacy-compat" omit-xml-declaration="yes" indent="yes" method="html" media-type="application/xhml" encoding="utf-8"/> <xsl:namespace-alias stylesheet-prefix="wxsl" result-prefix="xsl" /> <xsl:template match="head"> <!-- code omitted clarity --> <script type="text/template"> <wxsl:stylesheet version="1.0" xmlns:xsl="http://www.w3.org/1999/xsl/transform" xmlns="http://www.w3.org/1999/xhtml" xmlns:xhtml="http://www.w3.org/1999/xhtml" xmlns:eli="local"> ...
Read more