Lid driven cavity working

src/FiniteDifference/eq_unit_test.hpp

@@ -110,12 +110,12 @@ BOOST_AUTO_TEST_SUITE( eq_test_suite )
 BOOST_AUTO_TEST_CASE(lid_driven_cavity)
 {
 	// Domain
-	Box<2,float> domain({0.0,0.0},{1.0,1.0});
+	Box<2,float> domain({0.0,0.0},{3.0,1.0});
 
 	// Ghost
 	Ghost<2,float> g(0.01);
 
-	long int sz[] = {8,8};
+	long int sz[] = {256,64};
 	size_t szu[2];
 	szu[0] = (size_t)sz[0];
 	szu[1] = (size_t)sz[1];

src/Vector/Vector_eigen.hpp

@@ -123,12 +123,18 @@ class Vector<T,Eigen::Matrix<T, Eigen::Dynamic, 1>>
 		// set the staggered position of the property
 		openfpm::vector<comb<scheme::Sys_eqs_typ::dims>> stag_pos[sizeof...(pos)];
 
+		// interpolation points for each property
+		openfpm::vector<std::vector<comb<scheme::Sys_eqs_typ::dims>>> interp_pos[sizeof...(pos)];
+
 		// fill the staggered position vector for each property
 		stag_set_position<scheme::Sys_eqs_typ::dims,typename Grid_dst::value_type::type> ssp(stag_pos);
 
 		typedef boost::mpl::vector_c<unsigned int,pos ... > v_pos_type;
 		boost::mpl::for_each_ref<v_pos_type>(ssp);
 
+		interp_points<scheme::Sys_eqs_typ::dims,v_pos_type,typename Grid_dst::value_type::type> itp(interp_pos,stag_pos);
+		boost::mpl::for_each_ref<v_pos_type>(itp);
+
 		// shift the start and stop by the padding
 		grid_key_dx<scheme::Sys_eqs_typ::dims> start_k = grid_key_dx<scheme::Sys_eqs_typ::dims>(start) + pd.getKP1();
 		grid_key_dx<scheme::Sys_eqs_typ::dims> stop_k = grid_key_dx<scheme::Sys_eqs_typ::dims>(stop) + pd.getKP1();
@@ -154,7 +160,7 @@ class Vector<T,Eigen::Matrix<T, Eigen::Dynamic, 1>>
 
 			// Transform this id into an id for the Eigen vector
 
-			interp_ele<scheme,Grid_dst,Eigen::Matrix<T, Eigen::Dynamic, 1>,T,sizeof...(pos)> cp(key_dst,g_dst,v,key_src,g_map,stag_pos);
+			interp_ele<scheme,Grid_dst,Eigen::Matrix<T, Eigen::Dynamic, 1>,T,sizeof...(pos)> cp(key_dst,g_dst,v,key_src,g_map,interp_pos);
 
 			boost::mpl::for_each_ref<boost::mpl::range_c<int,0,v_size::value>>(cp);