[aspect-devel] Querying Aspect plugins directly from C++

[Jonathan Perry-Houts]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Ah, the shear_thinning test is an excellent starting point! That's
really close to what I want to do. I like the idea of including the
postprocessor in the same file as the material model. That way I can
have a sort of built-in unit test and keep a simple test.prm file
around to run it.

Thanks for the tips! I'll let you know if I come up with anything
general enough to be worth including with Aspect.

Cheers,
Jonathan

On 11/18/2013 02:24 PM, Wolfgang Bangerth wrote:
> On 11/18/2013 10:48 AM, Timo Heister wrote:
>> Dear Jonathan,
>> 
>> a problem is that a material model can potentially depend on
>> many things inside aspect (geometry, etc.) which will be
>> difficult for you to reproduce in a separate program.
>> 
>> I would create a .prm that sets up a tiny test problem, does 2
>> time steps, and writes out the results. Then you can compare
>> things like the statistics file. If you want to test something
>> specific, you can always write a postprocessor that generates a
>> file with the information you want (say look up viscosity using
>> representative_point and a range of temperatures).
>> 
>> Note that we have a lot of unittests in the tests/ directory that
>> do something similar.
> 
> As a starting point, take a look at tests/shear_thinning.prm 
> tests/shear_thinning.cc that implements a model where viscosity = 1
> / (1 + norm(eps(u))) We then simply solve this on a very coarse
> mesh (=very fast) and integrate this viscosity over the entire
> domain in a postprocessor.
> 
> If you use cmake to install aspect, just say cmake ...your usual
> command line... ctest to run all of the tests. The shear thinning
> test will be one of them. Most tests will fail because the output
> is slightly different from the one on our testing machine, but you
> get the idea.
> 
> Best W.
> 
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.14 (GNU/Linux)
Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/

iQEcBAEBAgAGBQJSiqkjAAoJEGe6xJ1FYRpRmbAH/R0XQHJq07B4+RW5FRWKlXkq
O95+CuoYjG9yrlBaYxtBd2RZL2tMh5Z/UDP8BcA32BeHl8EilrzZ1j5HN+pEeoz/
Ym+mKmkEvrx+apaZMJLiGibLWMoKlHBXR6F2oMPO2BR827xzsDj4FhxpNudrv6kQ
OAhJ2TOa99M+dTbn60WSQOx4Bpl3+8ckph2J2PDklvkZJ4g9qSkAacmdeBcRe2qs
gCWd7jBqVbpGcVv+4Dio7atLktW3EplLO8KvORCCCs9qIH9KlGzBOFodZUYBziQu
4PgLZOMMTrSw9iwEbixNnvN8fszhDElQc70kYq8BIZNk2xh/tORyZCR3B2PjjJg=
=K3ij
-----END PGP SIGNATURE-----