fortran - why move_alloc is not working in gfortran (4.6.3), but it is in ifort? -


i've read here move_alloc works since gfortran 4.2. i'm gfortran 4.6 installed in ubuntu 12.04 move_alloc not working! move_alloc used 5 times inside loop runs 10 times. after compiled (without error or warning) gfrotran, program runs 1 step (some prints verify mistake) of loop , shows "segmentation fault(kernel image recorded)'. however, when ifort used, program runs , works correctly. have tried use gfortran 4.4.6 in centos. both computers x86_64.

other important information: piece of code in subroutine, inside module, once don't know size of vector allocated move_alloc. these vector attribute intent (out) in subroutine. xray_all, yray_all , elem_all double precision , other integer. main , module in different files. here piece of code use move_alloc:

program main double precision,allocatable,dimension(:)::xrayall,yrayall (...)other allocatable variables call yyyy(....,ray_indent,xray_all,...) end program main   module xxxx subroutine yyyy j=1,10   <lots of calculation>       allocate( vec_aux( 1:(i+size(ray_indent) ) ) )      vec_aux(1:size(ray_indent))=ray_indent      vec_aux(size(ray_indent)+1:)=j         call move_alloc(vec_aux,ray_indent)       allocate( vec_auxreal( 1:(i+size(xray_all) ) ) )     vec_auxreal(1:size(xray_all))=xray_all     vec_auxreal(size(xray_all)+1:)=xray     call move_alloc(vec_auxreal,xray_all)                allocate( vec_auxreal( 1:(i+size(yray_all) ) ) )     vec_auxreal(1:size(yray_all))=yray_all     vec_auxreal(size(yray_all)+1:)=yray     call move_alloc(vec_auxreal,yray_all)                   elemsize=count(icol/=0);       allocate( vec_auxreal( 1:(elemsize+size(elem_all) ) ) )      vec_auxreal(1:size(elem_all))=elem_all      vec_auxreal(size(elem_all)+1:)=elem(1:elemsize)      call move_alloc(vec_auxreal,elem_all)                allocate( vec_aux( 1:(elemsize+size(icol_all) ) ) )      vec_aux(1:size(icol_all))=icol_all      vec_aux(size(icol_all)+1:)=icol(1:elemsize)      call move_alloc(vec_aux,icol_all)      allocate( vec_aux( 1:(elemsize+size(irow_all) ) ) )     vec_aux(1:size(irow_all))=irow_all     vec_aux(size(irow_all)+1:)=j+control           !     call move_alloc(vec_aux,irow_all) end end module xxxx end subroutine yyyy

i have found solution! in gfortran necessary add following if statement in 5 expressions:

 allocate( vec_auxreal( 1:(elemsize+size(elem_all) ) ) )  if (j/=1) vec_auxreal(1:size(elem_all))=elem_all   vec_auxreal(size(elem_all)+1:)=elem(1:elemsize)   call move_alloc(vec_auxreal,elem_all)

this happens because, in gfortran if vector still empty, not recognized nothing added. in ifort (tested in version 12.0), if statement not necessary program work.


Comments

Post a Comment

Popular posts from this blog

c# - DetailsView in ASP.Net - How to add another column on the side/add a control in each row? -

i have detailsview control on page used edit various fields of record, works in respect. i looking way add 1 column (and if works, why not more) right, absolutely read-only, show same fields of record comparison purposes. i aware there no obvious way such thing out of box detailsview . have looked other controls (transposed gridview , recommended formview , listview ), nothing satisfies. have special data binding setup using detailsview , can't out of without losing features. anyone on how "hack in" additional columns (for display only) on detailsview ? the solution have now, use second detailsview , visible set false in aspx. in code, make sure databind hidden detailsview hosts data third column first, initial detailsview named itemdetails . and in item created event, pass third column html rendering of hidden controls (in last code block) : protected void itemdetails_itemcreated(object sender, eventargs e) { if (dataite
Read more

javascript - firefox memory leak -

i have following script (i use raphaeljs also) draws circle, waits 1 second, draws circle again, etc (you can copy paste - functional, need include raphael raphaeljs.com): <div id="holder"> </div> <script src="raphael.js" type="text/javascript"></script> <script> var paper = raphael(10, 50, 320, 200); var rad = math.pi / 180; var finangle = 0; function sector(cx, cy, r, startangle, endangle, params) { var x1 = cx + r * math.cos(-startangle * rad), x2 = cx + r * math.cos(-endangle * rad), y1 = cy + r * math.sin(-startangle * rad), y2 = cy + r * math.sin(-endangle * rad); return ["m", cx, cy, "l", x1, y1, "a", r, r, 0, +(endangle - startangle > 180), 0, x2, y2]; } var path = sector(100,100,50,0,180); var arc = paper.path(path); function draw() { path = sector(100,100,50,0,finangle);
Read more

Trying to import CSV file to a SQL Server database using asp.net and c# - can't find what I'm missing -

i student , quite new programming, , given task using asp.net , c#, without being taught either. plan learn teach ourselves. the task stuck on making website import csv file sql server database of last year's room bookings around campus. the table has following columns: request_id; priority; module_id; day; start_time; length; park; students; room_code; status; semester_id; linked_request; week_1; week_2; week_3; week_4; week_5; week_6; week_7; week_8; week_9; week_10; week_11; week_12; week_13; week_14; week_15; ) in code, try make work 3 columns (request_id,priority,module) starters. when press import button, want read .csv file (with fixed name @ fixed directory) eg "1,2,3" - comma separating fields [delimiter] - , import sql server datab
Read more