Opened 12 years ago

Closed 12 years ago

Last modified 12 years ago

#368 closed defect (worksforme)

ANUGA parallel leaking along domain portion edges

Reported by: martins Owned by: steve
Priority: normal Milestone:
Component: Architecture and API Version:
Severity: major Keywords: Parallel
Cc:

Description

When I run a simulation using multiple processors I get this effect along the domain portion edges that quickly spreads ruining the simulation.

Attachments (5)

ThreeSimulations.zip (60.0 MB) - added by martins 12 years ago.
The same simulation run on one, two and three Cpus
TwoSimulationsNoCheck.zip (94.2 MB) - added by martins 12 years ago.
Two 6 Cpu Simulations at different mesh resolutions
ForSteve.zip (40.5 MB) - added by martins 12 years ago.
run_model_parallel.py (13.8 KB) - added by steve 12 years ago.
Script with elevation set for the sequential domain
Busselton.sww (14.3 MB) - added by martins 12 years ago.
Some offshore effects along domain segment edges, real or viewer artefact?

Change History (12)

Changed 12 years ago by martins

The same simulation run on one, two and three Cpus

comment:1 Changed 12 years ago by steve

Tried to run your simulation but got the error

Traceback (most recent call last):

File "run_model_parallel.py", line 27, in <module>

from setup_model import project, trigs_min

ImportError?: No module named setup_model

Can you upload the extra modules needed to run your problem.

Wondering if the check pointing may be causing the problem. Without check pointing do you get the same error?

Changed 12 years ago by martins

Two 6 Cpu Simulations at different mesh resolutions

Changed 12 years ago by martins

comment:2 Changed 12 years ago by steve

  • Owner changed from ole to steve

comment:3 Changed 12 years ago by steve

  • Status changed from new to assigned

Changed 12 years ago by steve

Script with elevation set for the sequential domain

comment:4 Changed 12 years ago by steve

The problem was created by the setting of the elevation quantity independently on each of the parallel domains.

A proper fix would involve setting up a parallel set_quantity method which communicates the vertex values from the processors owning the vertices to the neighbouring processors.

But a quick fix is to ensure the set up the elevation is in the sequential part of the code, i.e. where the original domain is created on processor 0. See attached script.

This problem only occurs for the elevation, as the other important quantities (stage, x and y momentum) are updated in the evolve loop, but elevation is assumed to be continuous and is not updated.

There should be a method for testing the continuity of elevation in the parallel case.

comment:5 Changed 12 years ago by steve

  • Component changed from Functionality and features to Architecture and API
  • Priority changed from high to normal

Changed 12 years ago by martins

Some offshore effects along domain segment edges, real or viewer artefact?

comment:6 Changed 12 years ago by martins

  • Resolution set to worksforme
  • Status changed from assigned to closed

comment:7 Changed 12 years ago by martins

With the latest version of the anuga_parallel code, there is a member function for the domain which will merge the parallel sww file.

Just put

domain.sww_merge()

after the evolve loop. (look at run_parallel_sw_merimbula.py in the anuga_parallel directory).

Can also use

domain.sww_merge(delete_old=True)

if you are brave and want to delete the individual sww files.

Look in anuga_parallel/test_parallel_frac_op.py

There is a parallel version of the Inlet_operator which must be imported from anuga_parallel.parallel_operator_factory

ie

from anuga_parallel.parallel_operator_factory import Inlet_operator

Note: See TracTickets for help on using tickets.