Hi Manav,
On Fri, Sep 20, 2013 at 10:22 PM, Manav Bhatia <[email protected]>wrote:
> Did you try this in 2D or 3D?
>
2D
> I do not get this error in 2D, but it shows up in 3D.
>
Hmm, OK.
>
> The following code is able to reproduce the issue at my end.
>
> int main (int argc, char* const argv[])
> {
> // Initialize libMesh.
> LibMeshInit init (argc, argv);
>
>
> {
> ParallelMesh m(init.comm());
> MeshTools::Generation::build_cube(m, 5,5,5, 0., 1., 0., 1., 0.,
> 1., HEX8);
> EquationSystems es(m);
> ExplicitSystem& s = es.add_system<ExplicitSystem>("mys");
> FEType fe;
> s.add_variable("v", fe);
> es.init();
> class Func: public FEMFunctionBase<Real> {
> public:
> virtual AutoPtr<FEMFunctionBase<Real> > clone () const {
> return AutoPtr<FEMFunctionBase<Real> > (new Func); };
> virtual Real operator() (const FEMContext&, const Point& p,
> const Real time = 0.) { return 0.; };
> virtual void operator() (const FEMContext&, const Point& p,
> const Real time,
> DenseVector<Real>& output) {
> output.zero(); };
>
>
> };
> Func ff;
> s.project_solution(&ff);
> }
> return 0;
> }
>
>
> Please let me know if you see the same behavior.
>
Thanks for the test case! I'll play with this in a little while and see if
I can reproduce.
Best,
Paul
------------------------------------------------------------------------------
LIMITED TIME SALE - Full Year of Microsoft Training For Just $49.99!
1,500+ hours of tutorials including VisualStudio 2012, Windows 8, SharePoint
2013, SQL 2012, MVC 4, more. BEST VALUE: New Multi-Library Power Pack includes
Mobile, Cloud, Java, and UX Design. Lowest price ever! Ends 9/22/13.
http://pubads.g.doubleclick.net/gampad/clk?id=64545871&iu=/4140/ostg.clktrk
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
