fortran - Difference of finite difference solution between two timesteps -

-

i have solution discretized differential equation given by

f(i)

where i spatial index. how can find difference between solution @ each adjacent time step? more clear:

the solution defined array

real,dimension(0:10) :: f

i discretize differential equation , solve stepping forward in time. if time index k, portion of code 

do k=1,25    = 1,10       f(i) = f(i+1)+f(i-1)+f(i)    end end

i can print solution, f(i) @ each time step k following code

print*, "print f(i) k=, k print "(//(5(5x,e22.14)))", f

how can find difference between solution @ each adjacent time step? is, time steps k+1,k. store value in new array g, has dimension given by

real,dimension(0:10) :: g

so trying find 

!g(i)=abs(f(i;k+1)-f(i;k))...not correct code.

how can this? way implement code? not sure how using if /then statements or whatever code need needed this. thanks

typically, in explicit time integration methods or iterative methods, have save last time-step last solution, current time-step solution , possibly more.

so have 

  real,dimension(0:10) :: f0, f

where f0 previous value

you iterate jacobi or gauss-seidel discretization:

f = f0  k=1,25     = 1,9       f(i) = f(i+1)+f(i-1)+f(i)    end     max_diff = maxval(abs(f-f0))     if (diff small enough) exit     f0 = f  end

if have time-evolving problem heat equation:

f = f0  k=1,25     = 1,9       f(i) = f0(i) + dt * viscosity * (f0(i+1)+f0(i-1)+f0(i))    end     max_diff = maxval(abs(f-f0))     f0 = f  end

Comments

Post a Comment

Popular posts from this blog

ios - RestKit 0.20 — CoreData: error: Failed to call designated initializer on NSManagedObject class (again) -

-
this error has shown on many times, seems. prominent answer issue seems restkit 0.20 — coredata: error: failed call designated initializer on nsmanagedobject class however, solution posted doesn't work me. let's start beginning: i want upload file server. want parse response make sure uploaded. if not, try again later. model has boolean variable called nsnumber *senttoserver purpose (in coredata, bool saved nsnumber ) . because don't know whether sending going successful, save entity coredata model first before sending via restkit. when post request successful, server returns json string this: {"id":14,"created":true} notice that, of course, id not set on coredata model yet, both id , created new informations model. here's model: @interface recording : nsmanagedobject @property (nonatomic, retain) nsstring * audiofilepath; @property (nonatomic, retain) nsdate * createdat; @property (nonatomic, retain) nsstring * deviceid; @...
Read more

java - Digest auth with Spring Security using javaconfig -

-
so i'm trying create digest authentication spring following documentation trying translate xml "requirements" in java requirements. let's have xml in docs: <bean id="digestfilter" class= "org.springframework.security.web.authentication.www.digestauthenticationfilter"> <property name="userdetailsservice" ref="jdbcdaoimpl"/> <property name="authenticationentrypoint" ref="digestentrypoint"/> <property name="usercache" ref="usercache"/> </bean> <bean id="digestentrypoint" class= "org.springframework.security.web.authentication.www.digestauthenticationentrypoint"> <property name="realmname" value="contacts realm via digest authentication"/> <property name="key" value="acegi"/> <property name="noncevalidityseconds" value="10"/> </b...
Read more

laravel - PDOException in Connector.php line 55: SQLSTATE[HY000] [1045] Access denied for user 'root'@'localhost' (using password: YES) -

-
my .env: app_env=local app_debug=true app_key=myappkey app_url=http://localhost db_connection=mysql db_host=localhost db_port=3306 db_database=gestiontaller db_username=root db_password='root' cache_driver=file session_driver=file queue_driver=sync redis_host=127.0.0.1 redis_password=null redis_port=6379 mail_driver=smtp mail_host=mailtrap.io mail_port=2525 mail_username=null mail_password=null mail_encryption=null i've tried password ", ' , without, restarting server every time. i've tried changing localhost 172.0.0.1. yes, password correct, running mysql -u root -p introducing password root in prompt can go mysql. users table: +------------------+-----------+ | user | host | +------------------+-----------+ | debian-sys-maint | localhost | | mysql.sys | localhost | | root | localhost | | root1 | localhost | +------------------+-----------+ i've added homestead.yaml this: database_port: 3306 ...
Read more