Domain decomposition

[mojtabaAmiraslanpour]

Hello,

I am trying to use the new functionality for domain decomposition for the pipflow case but I do not see my desired decomposition. This is what I am setting in the config.xml file:

<domain>
    ...
    <mABx>15</mABx>
    <mABy>8</mABy>
    <mABz>2</mABz>
</domain>

What I get is (10x6x4):

[mojtabaAmiraslanpour]

Any thoughts on this?

[gzavo]

Hi @mojtabaAmiraslanpour,

I'm not sure why you get this decomposition. This is the relevant code:

HemoCell/core/hemoCell.cpp

Lines 447 to 460 in ccccc63

	SparseBlockStructure3D sb = createRegularDistribution3D(management.getBoundingBox(),
	(*cfg)["domain"]["mABx"].read<int>(),
	(*cfg)["domain"]["mABy"].read<int>(),
	(*cfg)["domain"]["mABz"].read<int>());
	hlog << "(HemoCell) Domain management overwritten from config file." << endl;
	ExplicitThreadAttribution * eta = new ExplicitThreadAttribution(BlockToMpi);
	domain_lattice_management = new MultiBlockManagement3D(sb,eta,management.getEnvelopeWidth(),management.getRefinementLevel());
	lattice = new MultiBlockLattice3D<T,DESCRIPTOR>(*domain_lattice_management,
	defaultMultiBlockPolicy3D().getBlockCommunicator(),
	defaultMultiBlockPolicy3D().getCombinedStatistics(),
	defaultMultiBlockPolicy3D().getMultiCellAccess<T, DESCRIPTOR>(),
	new GuoExternalForceBGKdynamics<T, DESCRIPTOR>(1.0/param::tau));

It can be that this does not play nice together with the voxelizer. Maybe try another example where no voxelization is used, just to see if this is the case.

[mojtabaAmiraslanpour]

I managed to make it work. I was assuming it is already added to the available examples. I used initializeLattice and everything worked just fine.