next up previous contents
Next: ext : Spurious Up: AUTO Demos : Previous: AUTO Demos :

pvl : Use of the Subroutine PVLS.

 

Consider Bratu's equation

equation2866

with boundary conditions tex2html_wrap_inline10274 , tex2html_wrap_inline9798 As in demo exp, a solution curve requires one free parameter; here tex2html_wrap_inline9556 .

Note that additional parameters are specified in the user-supplied subroutine PVLS in file pvls.f, namely, tex2html_wrap_inline9564 (the tex2html_wrap_inline9458 -norm of tex2html_wrap_inline9574 ), tex2html_wrap_inline9562 (the minimum of tex2html_wrap_inline9576 on the space-interval [0,1] ), tex2html_wrap_inline9706 (the boundary value tex2html_wrap_inline10294  ). These additional parameters should be considered as ``solution measures'' for output purposes; they should not be treated as true continuation parameters.

Note also that four free parameters are specified in the AUTO-constants file r.pvl.1, namely, tex2html_wrap_inline9556 , tex2html_wrap_inline9564 , tex2html_wrap_inline9562 , and tex2html_wrap_inline9706 . The first one in this list, tex2html_wrap_inline9556 , is the true continuation parameter. The parameters tex2html_wrap_inline9564 , tex2html_wrap_inline9562 , and tex2html_wrap_inline9706 are overspecified so that their values will appear in the output. However, it is essential that the true continuation parameter appear first. For example, it would be an error to specify the parameters in the following order : tex2html_wrap_inline9564 , tex2html_wrap_inline9556 , tex2html_wrap_inline9562 , tex2html_wrap_inline9706 .

In general, true continuation parameters must appear first in the parameter-specification in the AUTO constants-file. Overspecified parameters will be printed, and can be defined in PVLS, but they are not part of the intrinsic continuation procedure.

As this demo also illustrates (see the UZR values in r.pvl.1), labeled solutions can also be output at selected values of the overspecified parameters.

 

COMMAND ACTION
mkdir pvl create an empty work directory
cd pvl change directory
@dm pvl copy the demo files to the work directory
cp r.pvl.1 r.pvl get the constants-file
@r pvl compute a solution branch
@sv pvl save output-files as p.pvl, q.pvl, d.pvl
Table 15.1: Commands for running demo pvl.

 


next up previous contents
Next: ext : Spurious Up: AUTO Demos : Previous: AUTO Demos :

Gabriel Lord
Tue Jan 25 09:54:03 GMT 2000