duoram-utils.h 17 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727
  1. enum step
  2. {
  3. write_out = 0,
  4. num_steps
  5. };
  6. size_t duoram_progress[step::num_steps] = { 0 };
  7. // const size_t expo = 5;
  8. // const size_t db_nitems = 1ULL << expo;
  9. const size_t number_of_writes = 2;
  10. // const size_t number_of_ind_reads = 5;
  11. typedef int64_t DB_t;
  12. struct DuORAM_Write
  13. {
  14. size_t shift;
  15. DB_t CW;
  16. };
  17. size_t target_index = 2;
  18. size_t bytes_written2 =0;
  19. DB_t final_cw;
  20. DB_t * DB;
  21. DB_t * updated_DB;
  22. DB_t * blinds;
  23. DB_t * updated_blinds;
  24. DB_t * blinded_DB;
  25. DB_t * blinded_DB_recv;
  26. DB_t * updated_blinded_DB_recv;
  27. DB_t M;
  28. DB_t * reading_temp;
  29. DB_t distinguised_value[number_of_writes];
  30. DB_t * b;
  31. DB_t * c;
  32. DB_t * d;
  33. DB_t one_shareA, one_shareB, one_shareC;
  34. int8_t * reading_b;
  35. int8_t * reading_c;
  36. int8_t * reading_d;
  37. int8_t * writing_b;
  38. int8_t * writing_c;
  39. int8_t * writing_d;
  40. void generate_flags_and_standard_basis_vectors2(bool party, size_t db_nitems)
  41. {
  42. for(size_t j = 0; j < db_nitems; ++j)
  43. {
  44. reading_b[j] = 0;
  45. reading_c[j] = 0;
  46. reading_d[j] = 0;
  47. if(j > 5 && j < 42)
  48. {
  49. reading_b[j] = 1;
  50. reading_c[j] = 1;
  51. reading_d[j] = 1;
  52. }
  53. if(j > 75 && (j % 3) == 0)
  54. {
  55. reading_b[j] = 1;
  56. reading_c[j] = 1;
  57. reading_d[j] = 1;
  58. }
  59. }
  60. if(party)
  61. {
  62. reading_b[target_index] = 1;
  63. reading_c[target_index] = 1;
  64. reading_d[target_index] = 1;
  65. }
  66. if(!party)
  67. {
  68. for(size_t j = 0; j < db_nitems; ++j)
  69. {
  70. reading_b[j] = -reading_b[j];
  71. reading_c[j] = -reading_c[j];
  72. reading_d[j] = -reading_d[j];
  73. }
  74. }
  75. for(size_t j = 0; j < db_nitems; ++j)
  76. {
  77. b[j] = j;
  78. c[j] = j;
  79. d[j] = j;
  80. writing_b[j] = 0;
  81. writing_c[j] = 0;
  82. writing_d[j] = 0;
  83. if(j > 5 && j < 42)
  84. {
  85. writing_b[j] = 1;
  86. writing_c[j] = 1;
  87. writing_d[j] = 1;
  88. }
  89. if(j > 75 && (j % 3) == 0)
  90. {
  91. writing_b[j] = 1;
  92. writing_c[j] = 1;
  93. writing_d[j] = 1;
  94. }
  95. }
  96. if(party)
  97. {
  98. writing_b[target_index] = 1;
  99. writing_c[target_index] = 1;
  100. writing_d[target_index] = 1;
  101. }
  102. if(party)
  103. {
  104. // printf("PARTY 1\n");
  105. // for(size_t j = 0; j < db_nitems; ++j)
  106. // {
  107. // if(writing_b[j]) pm0 += 1;
  108. // if(writing_c[j]) pm1 += 1;
  109. // if(writing_d[j]) pm2 += 1;
  110. // }
  111. }
  112. else
  113. {
  114. // printf("PARTY 0\n");
  115. // for(size_t j = 0; j < db_nitems; ++j)
  116. // {
  117. // if(writing_b[j]) pm0 -= 1;
  118. // if(writing_c[j]) pm1 -= 1;
  119. // if(writing_d[j]) pm2 -= 1;
  120. // }
  121. for(size_t j = 0; j < db_nitems; ++j)
  122. {
  123. writing_b[j] = -writing_b[j];
  124. writing_c[j] = -writing_c[j];
  125. writing_d[j] = -writing_d[j];
  126. b[j] = -b[j];
  127. c[j] = -c[j];
  128. d[j] = -d[j];
  129. }
  130. }
  131. if(party)
  132. {
  133. b[target_index] = -47920;
  134. c[target_index] = -47920;
  135. d[target_index] = -47920;
  136. }
  137. else
  138. {
  139. b[target_index] = 0;
  140. c[target_index] = 0;
  141. d[target_index] = 0;
  142. }
  143. }
  144. void setup(DB_t * DB, DB_t * updated_DB, DB_t * blinded_DB_recv, DB_t * blinds, DB_t * updated_blinds, DB_t * updated_blinded_DB_recv, size_t db_nitems, bool party)
  145. {
  146. for(size_t j = 0; j < db_nitems; ++j)
  147. {
  148. if(party)
  149. {
  150. DB[j] = 0;
  151. updated_DB[j] = 0;
  152. blinded_DB_recv[j] = 0;
  153. }
  154. else
  155. {
  156. DB[j] = 0;
  157. updated_DB[j] = 0;
  158. blinded_DB_recv[j] = 0;
  159. }
  160. }
  161. for(size_t j = 0; j < db_nitems; ++j)
  162. {
  163. blinds[j] = 0;
  164. updated_blinds[j] = blinds[j];
  165. updated_blinded_DB_recv[j] = blinded_DB_recv[j];
  166. }
  167. }
  168. void debug_(tcp::socket& in2, tcp::socket& sb, size_t db_nitems)
  169. {
  170. for(size_t j = 0; j < db_nitems; ++j)
  171. {
  172. DB_t debug_blinds2;
  173. boost::asio::read(in2, boost::asio::buffer(&debug_blinds2, sizeof(debug_blinds2)));
  174. assert(blinds[j] == debug_blinds2);
  175. }
  176. for(size_t jj = 0; jj < db_nitems; ++jj)
  177. {
  178. DB_t debug_refresh;
  179. boost::asio::write(sb, boost::asio::buffer(&updated_blinded_DB_recv[jj], sizeof(updated_blinded_DB_recv[jj])));
  180. boost::asio::read(sb, boost::asio::buffer(&debug_refresh, sizeof(debug_refresh)));
  181. assert(debug_refresh == DB[jj] + blinds[jj]);
  182. }
  183. DB_t DB_out;
  184. boost::asio::write(sb, boost::asio::buffer(&DB[target_index], sizeof(DB[target_index])));
  185. boost::asio::read(sb, boost::asio::buffer(&DB_out, sizeof(DB_out)));
  186. DB_out = DB_out + DB[target_index];
  187. std::cout << "DB_out = " << DB_out << std::endl;
  188. }
  189. DB_t print_reconstruction(tcp::socket& sb, DB_t output)
  190. {
  191. DB_t out_reconstruction;
  192. boost::asio::write(sb, boost::asio::buffer(&output, sizeof(output)));
  193. boost::asio::read(sb, boost::asio::buffer(&out_reconstruction, sizeof(out_reconstruction)));
  194. out_reconstruction = out_reconstruction + output;
  195. return out_reconstruction;
  196. }
  197. int read_database_shares(bool party, size_t db_nitems)
  198. {
  199. if(party)
  200. {
  201. int const in { open( "DB1", O_RDONLY ) };
  202. size_t r = read(in, DB, db_nitems * sizeof(DB_t));
  203. if(r < 0) {
  204. perror("Read error");
  205. close(in);
  206. return 1;
  207. }
  208. }
  209. if(!party)
  210. {
  211. int const in { open( "DB0", O_RDONLY ) };
  212. size_t r = read(in, DB, db_nitems * sizeof(DB_t));
  213. if(r < 0) {
  214. perror("Read error");
  215. close(in);
  216. return 1;
  217. }
  218. }
  219. return 0;
  220. }
  221. int read_flags(bool party, size_t db_nitems)
  222. {
  223. if(!party)
  224. {
  225. int const in0 { open( "party0_read_flags_b", O_RDONLY ) };
  226. size_t r = read(in0, reading_b, db_nitems * sizeof(reading_b[0]));
  227. if(r < 0) {
  228. perror("Read error");
  229. close(in0);
  230. return 1;
  231. }
  232. int const in1 { open( "party0_read_flags_c", O_RDONLY ) };
  233. r = read(in1, reading_c, db_nitems * sizeof(reading_c[0]));
  234. if(r < 0) {
  235. perror("Read error");
  236. close(in1);
  237. return 1;
  238. }
  239. int const in2 { open( "party0_read_flags_d", O_RDONLY ) };
  240. r = read(in2, reading_d, db_nitems * sizeof(reading_d[0]));
  241. if(r < 0) {
  242. perror("Read error");
  243. close(in2);
  244. return 1;
  245. }
  246. int const in0_w { open( "party0_write_flags_b", O_RDONLY ) };
  247. r = read(in0_w, writing_b, db_nitems * sizeof(writing_b[0]));
  248. if(r < 0) {
  249. perror("Read error");
  250. close(in0_w);
  251. return 1;
  252. }
  253. int const in1_w { open( "party0_write_flags_c", O_RDONLY ) };
  254. r = read(in1_w, writing_c, db_nitems * sizeof(writing_c[0]));
  255. if(r < 0) {
  256. perror("Read error");
  257. close(in1_w);
  258. return 1;
  259. }
  260. int const in2_w { open( "party0_write_flags_d", O_RDONLY ) };
  261. r = read(in2_w, writing_d, db_nitems * sizeof(writing_d[0]));
  262. if(r < 0) {
  263. perror("Read error");
  264. close(in2_w);
  265. return 1;
  266. }
  267. int const in0_w_ { open( "party0_write_b", O_RDONLY ) };
  268. r = read(in0_w_, b, db_nitems * sizeof(b[0]));
  269. if(r < 0) {
  270. perror("Read error");
  271. close(in0_w_);
  272. return 1;
  273. }
  274. int const in1_w_ { open( "party0_write_c", O_RDONLY ) };
  275. r = read(in1_w_, c, db_nitems * sizeof(c[0]));
  276. if(r < 0) {
  277. perror("Read error");
  278. close(in1_w_);
  279. return 1;
  280. }
  281. int const in2_w_ { open( "party0_write_d", O_RDONLY ) };
  282. r = read(in2_w_, d, db_nitems * sizeof(d[0]));
  283. if(r < 0) {
  284. perror("Read error");
  285. close(in2_w_);
  286. return 1;
  287. }
  288. }
  289. if(party)
  290. {
  291. int const in0 { open( "party1_read_flags_b", O_RDONLY ) };
  292. size_t r = read(in0, reading_b, db_nitems * sizeof(reading_b[0]));
  293. if(r < 0) {
  294. perror("Read error");
  295. close(in0);
  296. return 1;
  297. }
  298. int const in1 { open( "party1_read_flags_c", O_RDONLY ) };
  299. r = read(in1, reading_c, db_nitems * sizeof(reading_c[0]));
  300. if(r < 0) {
  301. perror("Read error");
  302. close(in1);
  303. return 1;
  304. }
  305. int const in2 { open( "party1_read_flags_d", O_RDONLY ) };
  306. r = read(in2, reading_d, db_nitems * sizeof(reading_d[0]));
  307. if(r < 0) {
  308. perror("Read error");
  309. close(in2);
  310. return 1;
  311. }
  312. int const in0_w { open( "party1_write_flags_b", O_RDONLY ) };
  313. r = read(in0_w, writing_b, db_nitems * sizeof(writing_b[0]));
  314. if(r < 0) {
  315. perror("Read error");
  316. close(in0_w);
  317. return 1;
  318. }
  319. int const in1_w { open( "party1_write_flags_c", O_RDONLY ) };
  320. r = read(in1_w, writing_c, db_nitems * sizeof(writing_c[0]));
  321. if(r < 0) {
  322. perror("Read error");
  323. close(in1_w);
  324. return 1;
  325. }
  326. int const in2_w { open( "party1_write_flags_d", O_RDONLY ) };
  327. r = read(in2_w, writing_d, db_nitems * sizeof(writing_d[0]));
  328. if(r < 0) {
  329. perror("Read error");
  330. close(in2_w);
  331. return 1;
  332. }
  333. int const in0_w_ { open( "party1_write_b", O_RDONLY ) };
  334. r = read(in0_w_, b, db_nitems * sizeof(b[0]));
  335. if(r < 0) {
  336. perror("Read error");
  337. close(in0_w);
  338. return 1;
  339. }
  340. int const in1_w_ { open( "party1_write_c", O_RDONLY ) };
  341. r = read(in1_w_, c, db_nitems * sizeof(c[0]));
  342. if(r < 0) {
  343. perror("Read error");
  344. close(in1_w);
  345. return 1;
  346. }
  347. int const in2_w_ { open( "party1_write_d", O_RDONLY ) };
  348. r = read(in2_w_, d, db_nitems * sizeof(d[0]));
  349. if(r < 0) {
  350. perror("Read error");
  351. close(in2_w);
  352. return 1;
  353. }
  354. }
  355. return 0;
  356. }
  357. int read_flags_for_reading(bool party, size_t db_nitems)
  358. {
  359. if(!party)
  360. {
  361. int const in0 { open( "party0_read_flags_b", O_RDONLY ) };
  362. size_t r = read(in0, reading_b, db_nitems * sizeof(reading_b[0]));
  363. if(r < 0) {
  364. perror("Read error");
  365. close(in0);
  366. return 1;
  367. }
  368. int const in1 { open( "party0_read_flags_c", O_RDONLY ) };
  369. r = read(in1, reading_c, db_nitems * sizeof(reading_c[0]));
  370. if(r < 0) {
  371. perror("Read error");
  372. close(in1);
  373. return 1;
  374. }
  375. int const in2 { open( "party0_read_flags_d", O_RDONLY ) };
  376. r = read(in2, reading_d, db_nitems * sizeof(reading_d[0]));
  377. if(r < 0) {
  378. perror("Read error");
  379. close(in2);
  380. return 1;
  381. }
  382. }
  383. if(party)
  384. {
  385. int const in0 { open( "party1_read_flags_b", O_RDONLY ) };
  386. size_t r = read(in0, reading_b, db_nitems * sizeof(reading_b[0]));
  387. if(r < 0) {
  388. perror("Read error");
  389. close(in0);
  390. return 1;
  391. }
  392. int const in1 { open( "party1_read_flags_c", O_RDONLY ) };
  393. r = read(in1, reading_c, db_nitems * sizeof(reading_c[0]));
  394. if(r < 0) {
  395. perror("Read error");
  396. close(in1);
  397. return 1;
  398. }
  399. int const in2 { open( "party1_read_flags_d", O_RDONLY ) };
  400. r = read(in2, reading_d, db_nitems * sizeof(reading_d[0]));
  401. if(r < 0) {
  402. perror("Read error");
  403. close(in2);
  404. return 1;
  405. }
  406. }
  407. return 0;
  408. }
  409. int read_flags_for_writing(bool party, size_t db_nitems)
  410. {
  411. if(!party)
  412. {
  413. int const in0_w { open( "party0_write_flags_b", O_RDONLY ) };
  414. size_t r = read(in0_w, writing_b, db_nitems * sizeof(writing_b[0]));
  415. if(r < 0) {
  416. perror("Read error");
  417. close(in0_w);
  418. return 1;
  419. }
  420. int const in1_w { open( "party0_write_flags_c", O_RDONLY ) };
  421. r = read(in1_w, writing_c, db_nitems * sizeof(writing_c[0]));
  422. if(r < 0) {
  423. perror("Read error");
  424. close(in1_w);
  425. return 1;
  426. }
  427. int const in2_w { open( "party0_write_flags_d", O_RDONLY ) };
  428. r = read(in2_w, writing_d, db_nitems * sizeof(writing_d[0]));
  429. if(r < 0) {
  430. perror("Read error");
  431. close(in2_w);
  432. return 1;
  433. }
  434. }
  435. if(party)
  436. {
  437. int const in0_w { open( "party1_write_flags_b", O_RDONLY ) };
  438. size_t r = read(in0_w, writing_b, db_nitems * sizeof(writing_b[0]));
  439. if(r < 0) {
  440. perror("Read error");
  441. close(in0_w);
  442. return 1;
  443. }
  444. int const in1_w { open( "party1_write_flags_c", O_RDONLY ) };
  445. r = read(in1_w, writing_c, db_nitems * sizeof(writing_c[0]));
  446. if(r < 0) {
  447. perror("Read error");
  448. close(in1_w);
  449. return 1;
  450. }
  451. int const in2_w { open( "party1_write_flags_d", O_RDONLY ) };
  452. r = read(in2_w, writing_d, db_nitems * sizeof(writing_d[0]));
  453. if(r < 0) {
  454. perror("Read error");
  455. close(in2_w);
  456. return 1;
  457. }
  458. }
  459. return 0;
  460. }
  461. int read_flags_for_updating(bool party, size_t db_nitems)
  462. {
  463. if(!party)
  464. {
  465. int const in0_w_ { open( "party0_write_b", O_RDONLY ) };
  466. size_t r = read(in0_w_, b, db_nitems * sizeof(b[0]));
  467. if(r < 0) {
  468. perror("Read error");
  469. close(in0_w_);
  470. return 1;
  471. }
  472. int const in1_w_ { open( "party0_write_c", O_RDONLY ) };
  473. r = read(in1_w_, c, db_nitems * sizeof(c[0]));
  474. if(r < 0) {
  475. perror("Read error");
  476. close(in1_w_);
  477. return 1;
  478. }
  479. int const in2_w_ { open( "party0_write_d", O_RDONLY ) };
  480. r = read(in2_w_, d, db_nitems * sizeof(d[0]));
  481. if(r < 0) {
  482. perror("Read error");
  483. close(in2_w_);
  484. return 1;
  485. }
  486. }
  487. if(party)
  488. {
  489. int const in0_w_ { open( "party1_write_b", O_RDONLY ) };
  490. size_t r = read(in0_w_, b, db_nitems * sizeof(b[0]));
  491. if(r < 0) {
  492. perror("Read error");
  493. close(in0_w_);
  494. return 1;
  495. }
  496. int const in1_w_ { open( "party1_write_c", O_RDONLY ) };
  497. r = read(in1_w_, c, db_nitems * sizeof(c[0]));
  498. if(r < 0) {
  499. perror("Read error");
  500. close(in1_w_);
  501. return 1;
  502. }
  503. int const in2_w_ { open( "party1_write_d", O_RDONLY ) };
  504. r = read(in2_w_, d, db_nitems * sizeof(d[0]));
  505. if(r < 0) {
  506. perror("Read error");
  507. close(in2_w_);
  508. return 1;
  509. }
  510. }
  511. return 0;
  512. }
  513. int read_flags_for_refreshing(size_t db_nitems)
  514. {
  515. int const in1_w { open( "party0_write_flags_c", O_RDONLY ) };
  516. size_t r = read(in1_w, writing_c, db_nitems * sizeof(writing_c[0]));
  517. if(r < 0) {
  518. perror("Read error");
  519. close(in1_w);
  520. return 1;
  521. }
  522. int const in2_w { open( "party1_write_flags_d", O_RDONLY ) };
  523. r = read(in2_w, writing_d, db_nitems * sizeof(writing_d[0]));
  524. if(r < 0) {
  525. perror("Read error");
  526. close(in2_w);
  527. return 1;
  528. }
  529. int const in1_w_ { open( "party0_write_c", O_RDONLY ) };
  530. r = read(in1_w_, c, db_nitems * sizeof(c[0]));
  531. if(r < 0) {
  532. perror("Read error");
  533. close(in1_w_);
  534. return 1;
  535. }
  536. int const in2_w_ { open( "party1_write_d", O_RDONLY ) };
  537. r = read(in2_w_, d, db_nitems * sizeof(d[0]));
  538. if(r < 0) {
  539. perror("Read error");
  540. close(in2_w_);
  541. return 1;
  542. }
  543. return 0;
  544. }
  545. int read_flags_for_generating_cancellation_terms(size_t db_nitems)
  546. {
  547. int const in2 { open( "party1_read_flags_d", O_RDONLY ) };
  548. size_t r = read(in2, reading_d, db_nitems * sizeof(reading_d[0]));
  549. if(r < 0) {
  550. perror("Read error");
  551. close(in2);
  552. return 1;
  553. }
  554. int const in2_ { open( "party0_read_flags_c", O_RDONLY ) };
  555. r = read(in2_, reading_c, db_nitems * sizeof(reading_c[0]));
  556. if(r < 0) {
  557. perror("Read error");
  558. close(in2_);
  559. return 1;
  560. }
  561. return 0;
  562. }
  563. void generate_random_distinguished_points(bool party)
  564. {
  565. if(party)
  566. {
  567. for(size_t j = 0; j < number_of_writes; ++j)
  568. {
  569. distinguised_value[j] = j + 2;
  570. }
  571. }
  572. if(!party)
  573. {
  574. for(size_t j = 0; j < number_of_writes; ++j)
  575. {
  576. distinguised_value[j] = j + 2;
  577. }
  578. }
  579. }
  580. // template<typename T1>
  581. // void rotate_all( T1 * orginal_vector1, T1 * rotated_vector1,
  582. // T1 * orginal_vector3, T1 * rotated_vector3,
  583. // T1 * orginal_vector5, T1 * rotated_vector5,
  584. // size_t rotate_by, size_t db_nitems)
  585. // {
  586. // for(size_t j = 0; j < db_nitems; ++j)
  587. // {
  588. // rotated_vector1[(j + rotate_by) % db_nitems] = orginal_vector1[j];
  589. // rotated_vector3[(j + rotate_by) % db_nitems] = orginal_vector3[j];
  590. // rotated_vector5[(j + rotate_by) % db_nitems] = orginal_vector5[j];
  591. // }
  592. // }
  593. DB_t dot_product_with_bool(DB_t D[], int8_t flags[], size_t db_nitems, size_t rotate_by = 0)
  594. {
  595. DB_t result = 0;
  596. for(size_t j = 0; j < db_nitems; ++j)
  597. {
  598. result = result + (D[(j + rotate_by) % db_nitems] * flags[j]);
  599. }
  600. return result;
  601. }