matrix - Parallelization Fortran 90 -

-

i started learn how parallelize programs in fortran , trying include in program.

the problem using lapack (to calculate eigenvalues , eigenvectors) , apparently creating conflict parallelization routine, because depending on matrix size calculates eigenvalues (wrong values small sizes, , doesn't compute them @ matrices equal or bigger 10x10).

when run program 1 thread works correctly. 

program nuevo  use omp_lib      implicit none   integer ::i,j,k,l,kk,jj,red,ok,info,lwork,stat1    integer :: m,n,q,rea,i0,j0,chunk,chunksize   !!!! rea= number of realizations >1   parameter(m=10,n=10,rea=12,lwork=2*m*n,ok=m*n,chunksize=rea/2)     complex*8 t,aa,tt,ww !! caution dimension of variables   integer, dimension(1)::  seed, old   integer, parameter ::  ikind=selected_real_kind(p=15)   real (kind=ikind) ::  pi   real ( kind = 8 ) wtime   real*8 d, dd(m,n),zz, x(2*m*n*rea), xx(m*n,2*rea), ss, hh, harvest(2*n*m*rea), w   complex*16 a(m*n,m*n), work(3*n-2,3*n-2)   real*8 rwork(2*m*n), b(m*n), f   !!!!!!!!!!!!random numbers!!!!!!!!!!!!!!!!!!!      seed(1) = 67937     call random_seed     call random_seed(size=k)     call random_seed(get=old(1:k))     call random_number(harvest)         x=harvest   !!!!!!!!!!!random matrix!!!!!!!!!!!!!!!!!!!!!!        xx=reshape( x, (/ m*n, 2*rea /) ) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! wtime = omp_get_wtime ( ) !call omp_set_num_threads(1)    t=1.0d0                  !!  hoppings    chunk=chunksize  !$omp parallel !$omp schedule(static,chunk) ordered     jj=0,2*(rea-1),2  !!for realizations   w=0.1d0              !! desorder ratio   ww=n   q=m*n   aa=(0.0d0,1.0d0) !! complex number    a=0d0    b=0d0    f=2/ww     j=1,(m*n)     if (xx(j,jj+1)<0.5)        ss=1.0d0      else        ss=-1.0d0      end if     a(j,j)=ss*0.5d0*w*xx(j,jj+2)        !!  diagonal terms slighty random modifications      end      !!!!!!!!!!!!!!!!!!!!!!!!!!!     !!!!!!  matrix    !!!!!!!!!     !!!!!!!!!!!!!!!!!!!!!!!!!!!    !! in x direction  (i odd)          i=1,m,2             j=1,n-2,2                a((j*m)+i,((j+1)*m)+i)=t             end          end   !!  in x direction  (i even)          i=2,m,2             j=1,n-1,2                a(((j-1)*m)+i,(j*m)+i)=t             end          end    !!  in y direction             j=1,n             i=1,m-1             a(((j-1)*m)+i,((j-1)*m)+i+1)=t*exp(pi*j*f*aa) !!  phase factor due magnetic field              end          end !! periodic boundaries in x direction            i=1,m,2            a(i,(i+m*n-m))=t          end  !!    periodic boundaries in y direction              j=1,n                a(((j-1)*m)+1,j*m)=t*(exp(-pi*j*f*aa))              end  !!=========conjugate===============           i0=1,m*n             j0=1,m*n                a(j0,i0)=dconjg(a(i0,j0))             end          end  call zheev('v', 'u', ok, a, ok, b, work, lwork, rwork, info) !!lapack routine     i=1,m*n         write(*,*) b(i)         write(*,*) ' '     end  end     !! realizations   !$omp end !$omp end parallel  end program nuevo


Comments

Post a Comment

Popular posts from this blog

Load Balancing in Bluemix using custom domain and DNS SRV records -

-
i'm trying load balancing of application in bluemix using custom domain , srv record in dns manager. created own doamin ccbluemix.dynu.com , added custom domain bluemix. created application in region , created route custom domain myapp.ccbluemix.dynu.com deployed same application in uk region , created route custom domain myapp.ccbluemix.dynu.com there well. now whenever user hits route want load balancing between , uk region alternate requests served , uk regions. for created srv records like _http._tcp.ccbluemix.dynu.com srv eu-gb.mybluemix.net [priority: 0] [weight: 50] 10 _http._tcp.ccbluemix.dynu.com srv mybluemix.net [priority: 0] [weight: 50] 10 but nothing happening , i'm not able access application using custom route @ all. if create 'a' record cc-load-bal-app.ccbluemix.dynu.com 75.126.81.68 i can access app using custom route ends hitting server. how can achieve load balancing using srv records? wrong creating srv records here? i read ne...
Read more

oracle - pls-00402 alias required in select list of cursor to avoid duplicate column names -

-
my cursor : select xzis.msisdn msisdn, xzis.in_account_number in_account_number xxobw_zte_in_sync xzis not exists (select 1 hz_contact_points hcp, hz_cust_accounts hca hcp.owner_table_id = hca.party_id , hcp.owner_table_name = 'hz_parties' , hcp.phone_number = substr (xzis.msisdn, -8) , nvl (hca.attribute12, 'nil') = xzis.in_account_number); i have given alias columns, still facing issue. suggestions? it's cursor. 1 has no such problems. problem is, when you're joining 2 tables, or apply wrong aliases in subselect, cursor may @ end return 2 columns same name/alias. the problem not need alias, no 2 columns can have same alias (the name counting alias, when no alias given). apparently has happened, far can tell, it's ...
Read more

python - Consider setting $PYTHONHOME to <prefix>[:<exec_prefix>] error -

-
Image
i getting same error when python manage.py runserver. there no such error when have ubuntu 15.10. got started when upgraded ubuntu 16.04. question might duplicate have tried solution provided question have applied command dpkg --configure -a, have done apt-get update,upgrade,dist-upgrade, clean, -f install. have reinstalled python2.7 , python3 no success though. here screenshot while trying run django application i got after re-installing python please me? bashrc file # ~/.bashrc: executed bash(1) non-login shells. # see /usr/share/doc/bash/examples/startup-files (in package bash-doc) # examples # if not running interactively, don't case $- in *i*) ;; *) return;; esac export pythonhome="/usr/bin/python" echo $virtualenvwrapper_python echo 'tushant' # don't put duplicate lines or lines starting space in history. # see bash(1) more options histcontrol=ignoreboth # append history file, don't overwrite shopt -s histappend # setting h...
Read more