python - Pyomo DAE toolbox cannot discretize the domain appropriately -


i'm using pyomo dae toolbox solve partial differential equation , domain length of 1 independent variable (say m.x) small (like 1e-6).

it seems me dae toolbox cannot discretize domain appropriately because discretized m.x had 2 points (i used continuousset.get_finite_elements() method check that) specified 10 finite element points in option (see code below).

m.x = continuousset(bounds=(0, 1e-6))  disc=transformationfactory('dae.finite_difference') disc.apply_to(m, nfe=10, wrt=m.x, scheme='backward')

i increased upper bound of m.x as

m.x = continuousset(bounds=(0, 1e-5))

and tried discretize again. time got 10 points in x domain.

however, when increased number of finite element points as

m.x = continuousset(bounds=(0, 1e-5)) disc.apply_to(m, nfe=100, wrt=m.x, scheme='backward')

the returned discretized domain x still had 10 points.

i want know if there 1e-6 tolerance or limitation dae toolbox. if there is, how can customize value?

thank in advance!

yes, there limit on domain length. discretization transformation rounds new discretization points 6 decimal points. because seeing issues same discretization point being added multiple times because of slight differences in floating point round off. @qtothec mentioned, best solution modify scaling bounds on continuousset around (0,1).

the number of decimal points include not option can modified user can go source code , change number here , here.


Comments

Post a Comment

Popular posts from this blog

What is happening when Matlab is starting a "parallel pool"? -

running parallel cpu processes in matlab starts command parpool() according documentation , function: [creates] special job on pool of workers, , [connects] matlab client parallel pool. this function takes bit of time execute, on order of 30 seconds. in other multi-cpu paradigms openmp, parallel execution seems totally transparent -- i've never noticed behavior analogous matlab (granted i'm not experienced parallel programming). so, happening between time parpool() called , when finishes executing? takes long? parallel computing toolbox enables run matlab code in parallel using several different paradigms (e.g. jobs , tasks, parfor , spmd , parfeval , batch processing), , run either locally (parallelised across cores in local machine) or remotely (parallelised across machines in cluster - either 1 own, or 1 in cloud). in of these cases, code run on matlab workers , copies of matlab without interactive desktop. if you're intending run on remo...
Read more

angular - DownloadURL return null in below code -

when change file updated in firebase storage when try downloadurl show null. filechange(event: any) { const imagefolder: string = this.employee.id; const filelist: filelist = event.target.files; const file: file = filelist[0]; const storageref = firebase.storage().ref().child(`${imagefolder}/profile.jpg`).put(file); console.log(storageref.snapshot.downloadurl); } the put() asynchronous. if want url after file uploaded have this: firebase.storage().ref().child(`${imagefolder}/profile.jpg`).put(file).then((snapshot) => { storageref = snapshot.downloadurl: console.log(snapshot.downloadurl); });
Read more

php - Cannot override Laravel Spark authentication with own implementation -

during login process of laravel spark need request token of remote api. api stores exact same users , passwords laravel spark application. therefore need username , non hashed password of user during authentication process. i thought overriding authenticated method solution problem. in routes/web.php override post /login endpoint , pointing endpoint own logincontroller: <?php namespace app\http\controllers\auth; use illuminate\support\facades\log; use illuminate\http\request; use laravel\spark\http\controllers\auth\logincontroller sparklogincontroller; class logincontroller extends sparklogincontroller { /** * create new login controller instance. * * @return void */ public function __construct() { parent::__construct(); } /** * handle successful authentication attempt. * * @param request $request * @param \illuminate\contracts\auth\authenticatable $user * @return response */ public fu...
Read more